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Quantum cluster approaches offer new perspectives to study the complexities of macroscopic cor- 
related fermion systems. These approaches can be understood as generalized mean-field theories. 
Quantum cluster approaches are non-perturbative and are always in the thermodynamic limit. 
Their quality can be systematically improved, and they provide complementary information to 
finite size simulations. They have been studied intensively in recent years and are now well es- 
tablished. After a brief historical review, this article comparatively discusses the nature and 
advantages of these cluster techniques. Applications to common models of correlated electron 
systems are reviewed. ^ 
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I. INTRODUCTION 
A. Brief history 



The theoretical description of interacting many- 
particle systems remains one of the grand challenges 
in condensed matter physics. Especially the field of 
strongly correlated electron systems has regained theo- 
retical and experimental interest with the discovery of 
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heavy Fermion compounds and high-temperature super- 
conductors. In this class of systems the strength of the 
interactions between particles is comparable to or larger 
than their kinetic energy, i.e. any theory based on a pcr- 
turbative expansion around the non-interacting limit is 
at least questionable. Theoretical tools to describe these 
systems are therefore faced with extreme difficulties, due 
to the non-perturbative nature of the problem. A large 
body of work has been devoted to a direct (numerically) 
exact solution of finite size systems using exact diago- 
nalization or Quantum Monte Carlo methods. Exact di- 
agonalization however is severely limited by the expo- 
nential growth of computational effort with system size, 
while Quantum Monte Carlo methods suffer from the sign 
problem at low temperatures. Another difficulty of these 
methods arises from their strong finite size effects, of- 
ten ruling out the reliable extraction of low energy scales 
that are important to capture the competition between 
different ground states often present in strongly corre- 
lated systems. 

Mean-field theories are defined in the thermodynamic 
limit and therefore do not face the finite size problems. 
With applications to a wide variety of extended systems 
from spin models to models of correlated electrons and/or 
bosons, mean-field theories are extremely popular and 
ubiquitous throughout science. The first mean-field the- 
ory which gained wide acceptance w as developed by P. 
Weiss for spin systems l|Weissl . Il907|l . The Curie- Weiss 
mean-field theory reduces the complexity of the thermo- 
dynamic lattice spin problem by mapping it onto that 
of a magnetic impurity embedded in a self-consistently 
determined mean magnetic field. 

Generally, mean-field theories divide the infinite num- 
ber of degrees of freedom into two sets. A small set of 
degrees of freedom is treated explicitly, while the effects 
of the remaining degrees of freedom are summarized as 
a mean-field acting on the first set. Here, by mean-field 
theory, we refer to the class of approximations which ac- 
count for the correlations between spatially localized de- 
grees of freedom explicitly, while treating those at longer 
length scales with an effective medium. Such local ap- 
proximations become exact in the limit of infinite coor- 
dination number or e quivalcntly infinite dimensions D 
l|Itzvkson and Drouffel il989) ; however non-local correc- 
tions become important in finite dimensions. The pur- 
pose of this review is to discuss methods for incorporating 
non-local corrections to local approximations. 

Many different local approximations have been de- 
veloped for systems with itinerant degrees of freedom. 
Early attempts focused on disordered syste ms, and in- 
cluded the virtual cr ystal approxim ation ijNordheimL 
Il981albl: iParmenteil Il955t (Schoen", '1969^ and the 
average-T matrix approxiin ation (,,Beebv an d EdwardJ 
Il962t TSchwartz et alV Il97l|) . However, the most suc- 
cessful local approximations for disordered systems is 
th e Coherent Po tential Appr oximation (CPA) develope d 
bv lSovenI lll967D and others lIShibaL Il97lt iTavlorl Il967t) . 
This method is distinguished from the others in that it 



becomes exact in both the limit of dilute and concen- 
trated disordered impurity systems, as well as the limit 
of infinite dimensions. 

There have been many attempts to extend the CPA 
formalism to correlated systems, starting with the Dy - 
namical CPA rPCPAl ofKakehashil ll2002li:ISumJ lll974l) . 
the XNCA of 'Kim et al. fl990V 'Kuramoto ^ 1985(l and 
the LNCA ofiGrewc C1987}; Grewc et al. (1988*). A great 
breakthrough was achieved with the formulation of the 
Dynamical Mean-Field Theory (DMFT) (for a review see 
1^0 rgcs et al, 1996; Pruschke et al, 199 5) in t he lim it of 
infinite dimension s bvlMetzner and VoUhardd |)1989 ') and 
iMuller-HartmannI l)l989b|) . The DCPA and the DMFT 
have been the most successful approaches and employ the 
same mapping between the cluster and the lattice prob- 
lems. They differ mostly in their starting philosophy. 
The DCPA employs the CPA equations to relate the im- 
purity solution to the lattice whereas in the DMFT the 
irreducible quantities calculated on the impurity are used 
to construct the lattice quantities. 

Despite the success of these mean-field approaches, 
they share the critical fiaw of neglecting the effects of 
non-local fluctuations. Thus they are unable to capture 
the physics of, e.g. spin waves in spin systems, local- 
ization in disordered systems, or spin-liquid physics in 
correlated electronic systems. Non-local corrections are 
required to treat even the initial effects of these phenom- 
ena and to describe phase transitions to states described 
by a non-local order parameter. 

The first attempt to a dd non- l ocal co rrections to mean- 
field theories was due to iBethd |I193,'tD by adding correc- 
tions to the Curie- Weiss mean-field theory. This was 
achieved by mapping the lattice problem onto a self- 
consistently embedded finite-size spin cluster composed 
of a central site and z nearest neighbors embedded in a 
mean- field. For small z, the resulting theory provides a 
remarkably large and accurate correction to t he transi- 
tion temperature llKikuchl IT951I: ISuzukiL [1981 . 

Many attempts have been ma de to apply si milar ideas 
to disordered electronic systems (iGonisl [l992() . Most ap- 
proaches were hampered by the difficulty of constructing 
a fully causal theory, with positive spectral functions. 
Several causal theories w ere developed including the em- 
bedded cluster meth od llGonisl . Il99^ and the molecu- 
lar CPA (MCPA ) bv lTsukadal l|1969|) (for a review see 
lDucastelleLll97l . These methods generally are obtained 
from the local approximation by replacing the impurity 
by a finite size cluster in real space. As a result these ap- 
proaches suffer from the lack of translational invariance, 
since the cluster has open boundary conditions and only 
the surface sites couple to the mean-field. 

Similar effort has been expended to find cluster ex- 
tensions to the DMFT, including most n otably the Dy- 
namical C luster Approximation (DCA) ijHettler et all 
I2000lll998() and the Cellular Dynamical Mean-Field The- 
ory (CDMFT) ijKotliar all 1200 1|) . Both cluster ap- 
proaches reduce the complexity of the lattice problem by 
mapping it to a finite size cluster self-consistently em- 
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bedded in a mean-field. As in the classical case, the 
self-consistency condition reflects the translationally in- 
variant nature of the original lattice problem. The main 
difference with their classical counterparts arises from the 
presence of quantum fluctuations. Mean-field theories for 
quantum systems with itinerant degrees of freedom cut 
off spatial fluctuations but take full account of tempo- 
ral fluctuations. As a result the mean-fleld is a time- 
or respectively frequency dependent quantity. Even an 
effective cluster problem consisting of only a single site 
(DMFT) is hence a highly non-trivial many-body prob- 
lem. CDMFT and DCA mainly differ in the nature of the 
effective cluster problem. The CDMFT shares an identi- 
cal mapping of the lattice to the cluster problem with the 
MCPA, and hence also violates translational symmetries 
on the cluster. The DCA maps the lattice to a periodic 
and therefore translationally invariant cluster. 

A numerically more tractable cluster approxima- 
tion to the ther r nodyn amic limit was developed by 
iGros and Valentil l)l994|) . In this formalism the self- 
consistent coupling to a mean-field is neglected. This 
leads to a theory in which the self-energy of an isolated 
finite size cluster is us ed to approxim a te the lattice prop- 
agator. As shown bv lSenechal et alj l|200(]|) . this cluster 
extension of the Hubbard-I approximation is obtained as 
the leading order approximation in a strong-coupling ex- 
pansion in the hopping amplitude and hence this method 
was named Cluster Perturbation Theory (CPT). 

Generally, cluster formalisms share the basic idea to 
approximate the effects of correlations in the infinite lat- 
tice problem with those on a finite size quantum cluster. 
We refer to this class of techniques as quantum cluster 
theories. In contrast to Finite System Simulations (FSS), 
these techniques are built for the thermodynamic limit. 
In this review we focus on the three most established 
quantum cluster approaches, the DCA, the CDMFT and 
the CPT formalisms. The CDMFT approach was origi- 
nally formulated for general, possibly non-orthogonal ba- 
sis sets. In this review we restrict the discussion to the 
usual, comple tely localized orthog onal basis set and refer 
the reader to lKotliar et al\ ()200l|) for the generalization 
to arbitrary basis sets. 

The organization of this article is as follows: To famil- 
iarize the reader with the concept of cluster approaches, 
we develop in section ILBI a cluster generalization of the 
Curie- Weiss mean-field theory for spin systems. Section 
nil sets up the theoretical framework of the CDMFT, 
DCA and CPT formalisms by presenting two derivations 
based on different starting philosophies. The derivation 
based on the locator expansion in Sec. III. Al ls analogous 
to the cluster generalization of the Curie- Weiss mean- 
field method and thus is physically very intuitive. The 
derivation based on the cluster approximation to dia- 
grams defining the grand potential in Sec. III.Bl is closely 
related to the reciprocal space derivation of the DMFT 
bv|Miillcr-Hartmann ( 1989b). The nature of the different 
quantum cluster approaches together with their advan- 
tages and weaknesses are assessed in Sec. III.CI Discus- 



sions of the effective cluster problem, generalizations to 
symmetry broken states and the calculation of response 
functions are presented in Sees. III.Dl III. El and lll.FI The 
remainder of this section is devoted to describe the ap- 
plication of the DCA formalism to disordered systems in 
Sec. lII.d and to a brief discussion of alternative methods 
proposed to introduce non-local corrections to the DMFT 
method in Sec. III. HI In Sec. IIIII we review the various 
perturbative and non-perturbative techniques available 
to solve the effective self-consistent cluster problem of 
quantum cluster approaches. We include a detailed as- 
sessment of their advantages and limitations. Although 
numerous applications of quantum cluster approaches to 
models of many-particle systems are found in the litera- 
ture, this field is still in its footsteps and currently very 
active. A large body of work has been concentrated on 
the Hubbard model. We review the progress made on 
this model in Section IIVI together with applications to 
several other strongly correlated models. Finally, Sec. El 
concludes the review by stressing the limitations of quan- 
tum cluster approaches and proposing possible directions 
for future research in this field. 



B. Corrections to Curie-Weiss theory 

As an intuitive example of the formalism developed in 
the next sections we consider a systematic cluster exten- 
sion of the Curie- Weiss mean-field theory for a lattice 
of classical interacting spins. This discussion is espe- 
cially helpful to illustrate many new aspects of cluster 
approaches as compared to finite size simulations. The 
quality of this approach, and its convergence and critical 
properties will be demonstrated with a simple example, 
the one-dimensional Ising model 

= CTiCTi+i - (Ti (1) 

i i 

where <Ti = ±1 are classical spins, h is an external mag- 
netic field and the exchange integral J > acts between 
nearest neighbors only, favoring ferromagnetism. The 
generalization of this approach to higher dimensions and 
quantum spin systems is straightforward. 

We start by dividing the infinite lattice into N/Nc clus- 
ters of size Nc (see Fig. O with origin x and the exchange 
integral Jy into intra- (Jc) and inter-cluster {S3) parts 

3{xt - Xj) = Jc(55,,5^. + SJ{xi - Xj) (2) 

where each of the terms is a matrix in the Nc cluster sites. 
The central approximation of cluster theories is to retain 
correlation effects within the cluster and neglect them 
between the clusters. A natural formalism to implement 
this approximation is the locator expansion. The spin- 
susceptibility Xij = Pii'^i'^j) ~ where (3 =1/T 
is the inverse temperature, can be written as a locator 
expansion in the inter-cluster part <5J of the exchange 
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interaction, around the cluster limit x° = xi^-^ = 0) 



X°^xi,xj + X° 2^ 53{xi - xi)x{xi - Xj) (3) 



where we used again a matrix notation in the Nc cluster 
sites. By using the translational invariance of quantities 
in the superlattice a;, this expression can be simplified in 
the reciprocal space q of i to 



x(«)-X°+X°'5J(g)x('?) 



(4) 



This locator expansion has two well-defined limits. For 
an infinite size cluster it recovers the exact result since 
the surface to volume ratio vanishes making S3 irrelevant, 
and thus x = X°- For a single site cluster, Nc = 1, 
it recovers the Curie- Weiss mean-field theory. This is 
intuitively clear since for Nc = 1 fluctuations between all 
sites are neglected. With the susceptibility of a single 
isolated site x° = ^/T and 5J{q = 0) = J{q = 0) = J, 
we obtain for the uniform susceptibility 



xiq - 0) 



1 



1 



1/X°-J(9 = 0) T-Tc 



(5) 



the mean-field result with critical temperature Tc — J. 

For cluster sizes larger than one, translational sym- 
metry within the cluster is violated since the clusters 
have open boundary conditions and 53 only couples sites 
on the surface of the clusters. As detailed in the next 
section, this shortcoming can be formally overcome and 
translational invariance restored by considering an anal- 
ogous expression to the locator expansion in the 
Fourier space Q of the cluster 



xiQ,q) = x°{Q) + x°{Q)5J{Q,q)x{QA) 
1 



l/x°(Q)-<5J(g,g) 



(6) 



with analogous relations for the intra- and inter-cluster 
parts of J 



(5J(Q,g) = J{Q + q)-J{Q) 



(7) 

(8) 



Here, g is a vector in the reciprocal space of x, and Q 
is a vector in the reciprocal space of the cluster sites. 
The Fourier transform of the exchange integral is given 
by J{Q + q) = Jcos{Q + q), the intra-cluster exchange 
is J{Q), while the inter-cluster exchange is SJ{Q,q). As 
we will see in the next section, the resulting formalism 
is analogous to the dynamical cluster approximation for 
itinerant fermion systems. 

In analogy to the Curie- Weiss theory, the lattice sys- 
tem can now be mapped onto an effective cluster model 
embedded in a mean-field since correlations between the 
clusters are neglected. The susceptibility restricted to 



cluster sites is obtained by averaging or coarse-graining 
over the superlattice wave- vectors q 



/x°{Q) - r(g) 



with the hybridization function 



r(Q) 



N 



Y.JJ\QA)x{QA) 



^ + 7tEqSJ{Q,q)x{Q,q) 



(9) 



(10) 



This follows from the fact that the isolated cluster suscep- 
tibility x°iQ) docs not depend on the integration variable 
q in Eq. (01 . 

This expression defines the effective cluster model 
He = - ^ J(Q)a(Q)a(-g) - /ia(Q = 0) 



-^<5J(Q,g)(7(Q)(a(-Q 
Q,q 



q)) , (11) 



where cr{Q) {(j{q)) denotes the cluster (lattice) Fourier 
transform of cTj and (. . .) the expectation value calcu- 
lated with respect to the cluster Hamiltonian Tic- As in 
the Curie- Weiss theory, the cluster model is used to self- 
consistently determine the order parameter {a{Q + q)) = 
{a{Q))d{q) in the ferromagnetic state. In the paramag- 
netic state, the susceptibility calculated in the cluster 
model takes the same form as the coarse-grained result 
Eq. obtained from the locator expansion. 

The uniform susceptibility xiQ = 0,g = 0) contains 
information about the nature of this cluster approach, its 
critical properties and its convergence with cluster size. 
The sum in Eq. © may be solved analytically 



J(Q = 0) = J{Nc/n) sin(7r/iVc) ■ 



(12) 



The isolated cluster susceptibility x°{Q) can also be cal- 
culated analytically by us ing the transfer matrix method 
to give (|GoIdenfeI3,[l99l 



X°(Q = 0)=/3exp(2if) 



1 - (tanh(ii:)) 
1 + (tanh(if ))^^ 



(13) 



where K = (3J{Q = 0) = /3J(7Vc/7r) sin(7r/iVc). 
With these expressions the uniform lattice susceptibility 
Eq. 10 becomes 



X(T) 



1 



(14) 



i/x°(Q = o)-<5J(0 = o,g = o) 
1 

~ 1/X°(0 = 0) - J (1 - (Nc/n) sm{TT/Nc)) ' 

The cluster estimate of the lattice susceptibility interpo- 
lates between the Curie- Weiss result and the exact lattice 
result as Nc increases. It may be used to reveal some of 
the properties of cluster approximations and to compare 
the cluster results to both the finite-size calculation and 
the exact result in the thermodynamic limit. 
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Cluster mean- field 
Finite size 



20 
N 



40 



Figure 1 The cluster and finite-size estimates of the uniform 
lattice susceptibility versus cluster size when J = 1 and T — 
0.7. 



First, both the cluster mean- field resuh Eq. H14|l and 
the finite-size result Eq. with K — (3 J may be re- 
garded as an approximation to the thermodynamic re- 
sult. However, as illustrated in Fig. ^ the cluster mean- 
field result converges more quickly as a function of cluster 
size Nc than the finite size result. This reflects the supe- 
rior starting point of the cluster approximation compared 
to the finite-size calculation. The cluster approximation 
is an expansion about the mean-field result, whereas the 
finite-size calculation is an expansion about the atomic 
limit. 

It is instructive to explore the convergence of the clus- 
ter result analytically. For large the character of the 
susceptibility Eq. (|14(l can be split into three regimes. At 
very high temperatures 



xiT) 



1 



for T > J 



(15) 



where 6 ~ 2J-t-^ (^) ' intermediate temperatures, 



X(T) « /3e^^^ (^1 - ^ (^^ j j for J»r»T,. 

(16) 

The true critical behavior of the system can be resolved 
by studying the properties of this intermediate temper- 
ature regime. At both high and intermediate tempera- 
tures, the susceptibility differs from the exact result by 
corrections of order 0{1/N^). In general, cluster meth- 
ods with periodic boundary conditions have corrections 
of order O^l/L^), where Lc = nI^^ is the linear size of 
the cluster. 

At low temperatures, very close to the transition to 
the ferromagnetic state, deviations from the exact result 
are far larger. Here, for large clusters 



X{T) 



N, 



T-Tc 



(17) 



with the critical temperature Tc > 0, whereas the exact 
susceptibility in this regime x{T) ~ /3exp (2/3 J) does not 
diverge until zero temperature. This discrepancy is ex- 
pected in cluster approximations, since they treat long 
length scales which drive the transition in a mean-field 
way and therefore neglect long wave-length modes which 
eventually suppress the transition. Hence, cluster ap- 
proximations generally predict finite transition tempera- 
tures independent of dimensionality due to their residual 
mean-field character. With increasing cluster size how- 
ever, the transitions are expected to be systematically 
suppressed by the inclusion of longer-ranged fluctuations. 

For cluster sizes larger than one, all three regions are 
evident in the plot of the cluster mean-field estimate of 
the inverse susceptibility, shown in the inset to Fig. |21 
For Nc = 8 and Nc = 16, the high and low tempera- 
ture parts are seen as straight lines on the plot in the 
inset, with the crossover region in between. In numeri- 
cal simulations with significant sources of numerical noise 
especially close to the transition, it is extremely difficult 
to resolve the true low-temperature mean-field behav- 
ior. Exponents extracted from fits to the susceptibility 
in these simulations will more likely reflect the properties 
of the intermediate temperature regime. 

Despite the large deviations of the cluster result from 
the exact result low temperatures, we may still ex- 
tract the correct physics through finite size extrapola- 
tion. In general, for a system where the correlations 

build like ^ 



T* - aL. 



This behavior 



we expect Tc 

where T* is the exact transition temperature, Lc is 
the li near cluster s ize, and a is a positive real con- 
stant ijSuzuki Il986tl . However, for the ID Ising system, 
^ \ exp (/3 J) , so more care must be taken. Fortu- 
nately, an analytic expression for the transition temper- 
ature may be extracted from Eq. (|14|l . For large clus- 

ters, T, « JiV, (^1 (^)'-^(^)^ 

is shown in the main frame of Fig. |21 with the circles de- 
picting the numerical values for Tc and the solid line their 
asymptotic behavior. 



II. QUANTUM CLUSTER THEORIES 

In this section we provide two derivations of quantum 
cluster approaches for systems with itinerant quantum 
degrees of freedom. The locator expansion in Sec. III. A! 
is analogous to the cluster extension of the Curie- Weiss 
mean-field theory developed in the preceding section. 
Sec. lII.Bl provides a microscopic derivation based on clus- 
ter approximations to the thermodynamic grand poten- 
tial. A detailed discussion of the nature of quantum 
cluster approaches and the effective cluster model is pre- 
sented in Sees. III. CI and III.DI Generalizations for sym- 
metry broken phases, the calculation of susceptibilities 
and the application to disordered systems is explained in 
Sees. III. El III.FI and III.GI and a brief discussion of alter- 
native cluster methods is presented in Sec. III. Hi 
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Figure 2 Transition temperature for the ID Ising model ver- 
sus the inverse cluster size when J = 1 obtained with the 
cluster mean-field approach. For large clusters, Tc ~ ^/Nc, 
shown as a solid line. Inset: The inverse susceptibility ver- 
sus temperature. The Curie- Weiss behavior at low T where 
l/x{T) is linear in T, illustrates that the transition is always 
mean-field like. The true critical behavior of the transition is 
seen in the crossover region near T = J. 



A. Cluster approximation to the locator expansion 

In this section, we derive a number of cluster for- 
malisms for itinerant many body systems using an anal- 
ogous approach to that discussed in Sec. II. Bl for classical 
spin systems. For simplicity we assume in this section 
that no symmetry breaking occurs; the treatment of sym- 
metry broken phases is discussed in Sec. III. El The basic 
idea is to write down a locator expansion, i.e. an expan- 
sion in space around a finite-size cluster. This approach 
is not only intuitive but also allows us to assess the nature 
of quantum cluster approximations. As with their classi- 
cal counterparts, quantum cluster theories approximate 
the lattice problem with many degrees of freedom by an 
effective cluster problem with few degrees of freedom em- 
bedded in an external bath or mean-field created from the 
remaining degrees of freedom. By neglecting correlations 
that extend beyond the cluster size, one can then for- 
mulate a theory in which the lattice system is replaced 
by an effective cluster embedded in a mean- field host. 
While the formalism derived here is analogous to the for- 
malism discussed in Sec. II.BI for spin systems, there are 
significant differences. Since we are dealing with itiner- 
ant fermions, the theory is built upon the single-particle 
Green function instead of the two-particle spin correla- 
tion function, and the mean-field is dynamical due to the 
itinerant nature of the particles. 

This derivation is illustrated on the example of the 
extended Hubbard model 



ij^a ij.aa-' 



ij ' ^la ' ''J (7' 



(18) 



Here i and j are lattice site indices, the operators 
(Cj^) create (destroy) an electron with spin a on site z. 



Tiia = c|g.Cj^ is their corresponding number density, and 
Uij denotes the Coulomb repulsion between electrons on 
sites i and j. The hopping amplitude between sites i and 
j is denoted by ty, its local contribution ta — to and 
its Fourier transform to reciprocal space is the dispersion 
Ek- In this section we limit the discussion to the regular 
Hubbard model with a purely local interaction Uij — 
U{1 — 5aa')5ij. The more general case of finite non-local 
interactions Uij for i ^ j is discussed in Sec. III.BI 

The central quantity upon which we build the loca- 
tor expansion is the single-particle thermodynamic Green 
function (r is the imaginary time, T-r the corresponding 
time ordering operator, (3 = \/T the inverse temperature 
and oJn are the fermionic Matsubara frequencies) 



Gij.cr{iuJn) = / dTe''^"''Gij^a{T) 



(2n + l)7r 



(19) 



(20) 



or respectively its analytical continuation Gij^(j{z) = 

{{cifj,c^j„)) z to complex frequencies z. 

To set up a suitable notation for cluster schemes, we 
divide the I?-dimensional lattice of N sites into a set of 
finite-size clusters each with Nc sites of linear size such 
that Nc = , and resolve the first Brillouin zone into 
a corresponding set of reduced zones which we call cells. 
This notation is illustrated in Fig. O for iVc = 4 site clus- 
ters. For larger cluster sizes and more complex cluster 
geometries we refer the reader to Jarrcll et al. (2001b). 
Care should be taken so that the point group symmetry 
of the clusters does not differ too gr eatly from that of 
the original lattice ijBetts and Stewa rt. 1997). We use 
the coordinate x to label the origin of the clusters and 
X to label the Nc sites within a cluster, so that the site 
indices of the original lattice x = X + x. The points x 
form a superlattice with a reciprocal space labeled by k. 
The reciprocal space corresponding to the sites X within 
a cluster shall be labeled K, with — na ■ 2Tr /Lc and 
integer Ua- Then the wave- vectors in the full Brillouin 
zone are given by k = K -I- k. 

With these conventions, the Fourier transforms of a 
given function /(X,x) for intra- and inter-cluster coor- 
dinates are defined as 



/(X,5c) = -^^e^k.x^(x,k) (21) 

k 

/(X,k) = ^e-*-/(x,x) (22) 

X 

/(X,k) ^ ^^e^(K+'')-^/(K,k) (23) 

^ K 

/(K,k) - 5]e-^(K+'^)-^/(X,k). (24) 



To separate out the cluster degrees of freedom, the 
hopping amplitude t and the self-energy E (defined from 
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Figure 3 Definition of tfie coordinates in real (left) and recip- 
rocal (right) space illustrated for Nc = 4: site clusters. The 
origin of a cluster is labeled by x, the sites within a clus- 
ter by X. The reciprocal space to X is labeled by K, the 
wave-vectors of the superlattice, i.e. within a cell by k. 



the Green function via the Dyson equation G 



Gn 



S vifith the non-interacting Green function Gq) is split 
into intra- and inter-cluster parts 



t(Xj - Xj) 



<5t(x, 



3) (25) 
S(xj-Xj,2;) Sc(z)(5i,,i^. + ^S(xi - Xj,2;) . (26) 



All the quantities are Nc x Nc matrices in the cluster sites, 
tc = t(x — 0) and Tlc{z) — S(x = 0, z) are the intra- 
cluster hopping and self-energy, while (5t(x) and 5S(x, z) 
are the corresponding inter-cluster quantities only finite 
for X ^ 0. 

With these definitions we write the Green function 
using a locator expansion, an expansion in 5t and 5Ti 
around the cluster limit. In matrix notation in the Nc 
cluster sites it reads 



G(x, 



g(^)^iiij + g(2) ^['5t(xj - Xi 

I 



+ -x,,z)]G(i,-x„z) (27) 

where the Nc x Nc matrix 

g(z) = [(z + ^)l-te-Se(z)]-i (28) 

is the Green function of the cluster decoupled from the re- 
mainder of the system (/z is the chemical potential) . Since 
translational invariance in the superlattice x is preserved, 
this expression may be simplified by Fourier transforming 
the inter-cluster coordinates to give 

G(k, z) = g(z) -f g(z)[(5t(k) + 5S(k, z)]G(k, z) . (29) 

The central approximation that unites all cluster for- 
malisms is to truncate the self-energy to the cluster, by 
neglecting (5S to arrive at 



G(k,z) = g(z)+g(z)<5t(k)G(k,z) 
= [g-i(z)-<5t(k)]-i. 



(30) 



As we discuss below, this approximation corresponds to 
truncating the potential energy to the cluster while keep- 
ing all the contributions to the kinetic energy. Therefore 



quantum cluster approaches are good approximations to 
systems with significant screening, where non-local corre- 
lations are expected to be short-ranged. Quantum clus- 
ter approaches are particularly powerful because the re- 
maining self-energy term Tic{z) implicitly contained in 
the propagator ^{z) in Eq. I|30|l is restricted to the clus- 
ter degrees of freedom. Hence it can be calculated non- 
perturbatively in an effective cluster model as a func- 
tional Sc(^) = -^[6(2)], where 



(31) 



is the k-averaged or coarse-grained G(k, z), i.e. the 
Green function restricted to the cluster. As discussed 
in the next section this approximation is consistent with 
neglecting inter-cluster momentum conservation, i.e. ne- 
glecting the phase factors e*'''^ on the vertices of the self- 
energy diagrams. 

Using the expression H3U|I for the lattice Green function 
and the fact that g(z) does not depend on k, the coarse- 
grained Green function G can be written as 



G(z) = [g-i(z)-r(z)]-i , 

with a hybridization function T defined by 



(32) 



r(z) = 



l + §E'^tWG(k, 



N 



§5]a(k)G(k,z)<5t(k) 



N 



(33) 



Its physical content is that of an effective amplitude for 
fermionic hopping processes from the cluster into the host 
and back again into the cluster. The denominator in 
Eq. H33(l is a correction that excludes the cluster from 
the effective medium. T{z) thus plays an analogous role 
to that of the internal magnetic field in mean-field ap- 
proximations of spin systems. However, due to the itin- 
erant character of the fermionic degrees of freedom, it is 
a dynamical quantity. 

Both the CPT and the CDMFT formalisms may be 
defined at this point. A self-consistent set of equations 
is formed from G as a functional of Sc using Eq. H3U|I 
together with Eq. H28|l. and with an appropriate choice 
of a cluster solver (see Sec. IIII|) . Sc as a functional of 
G. In the CDMFT approximation the hybridization T is 
determined self-consistently with Eq. (|33|l . i.e. from the 
translational invariance of the super-lattice. The result- 
ing self-consistency cycle is discussed in Sec. III. B. 41 The 
CPT formalism is obtained when T is neglected. The 
Green function g(z) then becomes the Green function 
of an isolated cluster and the CPT result for the lat- 
tice Green function is obtained immediately via Eq. H3U|) 
without self-consistency. Thus the renormalization of the 
cluster degrees of freedom due to the coupling to the host 
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described by T is neglected in the CPT but included in 
the CDMFT formahsm. 

The DCA formalism may be motivated by the de- 
mand to restore translational invariance within the clus- 
ter. Since the inter-cluster hopping i5t(k) is finite for 
sites on the surface of the cluster and zero for bulk sites, 
only surface sites hybridize with the host. Hence trans- 
lational invariance with respect to the cluster sites X is 
violated. The cause of this violation can be seen by rep- 
resenting the hopping integral t(k) as the intra-cluster 
Fourier transform of the dispersion using Eq. (|23|l . 



[t(k)]x.x. 



, «(K+k)-(X,-X,) 



^ K 



(34) 



The violation of translational symmetry is caused by the 
phase factors e*'' '-^'"-^^' associated with the superlattice 
wave-vectors k. Thus translational symmetry can be re- 
stored by neglecting these phase-factors, or equivalently, 
by multiplying [t(k)]xiXj with the k-dependent phase 

g-ik-(Xi-X,)^ 

[ti5CA(k)]x.x, = [t(k)]x.x,e-'''-(^'-^^) 

^^,.K,x.-x,),^^^. (35) 



K 



Since t^cA is fully cycHc in the cluster sites, the DCA 
intra- and inter-cluster hopping integrals can be written 
as cluster Fourier transforms 



tc.DCAjXiX,- 



(36) 



5t 



DCA 



(k) 



XiX,- 



K 

^ ±^e'K.(x.-x,)^^(K + kp7) 



with 



K 



TV Z^^K+k 
k 



5i(K + k) = fK+k-^K. 



(38) 
(39) 



Since the DCA intra- and inter-cluster hopping inte- 
grals retain translational invariance within the cluster, 
the DCA cluster self-energy Sc and hybridization func- 
tion r are translationally invariant. The lattice Green 
function, Eq. (|3U|I hence becomes diagonal in cluster 
Fourier space 

G(K + k,z) = 5(K,z) + 5(K,z)(5i(K + k)G(K + k,z) 

- ^ ^ (40) 

(7-i(K,z)-<5<(K + k) 

with the Green function decoupled from the host 

g(K,z) = [z-eK + Ai-Sc(K,z)]-i. (41) 



Along the lines presented above, the DCA cluster self- 
energy S]c(K, z) is calculated as a functional of the 
coarse-grained Green function 



N 



1 



ff-i(K,z)-r(K,z) 



which defines the DCA hybridization function 



r(K,2) 



f Ek'^^'(K + k)G(K + k,z) 
l + ^Ek'5^(K + k)G(K + k,z 



(42) 



(43) 



The self-consistent procedure to determine the DCA clus- 
ter self-energy Sc(K, z) is analogous to CDMFT and dis- 
cussed in detail in Sec. III. B. 41 

This locator expansion yields a very natural physi- 
cal interpretation of cluster approximations. We note 
that the potential energy may be written as Tr (SG) 
l|Fetter and WaleckaL Il97l|) . where the trace runs over 
cluster sites, superlattice wave- vectors, frequency and 
spin. As detailed above, the central approxima- 
tion of cluster expansions is the neglect of the term 
(5S(k, 2;)G(k, z) in Eq. (|25J). Thus the approximation 
(5S(z) = essentially neglects the inter-cluster correc- 
tions to the potential energy in all calculated lattice 
quantities. On the other hand, the kinetic energy is iden- 
tified as Tr(tG). Since its inter-cluster contribution is 
not neglected, the kinetic and potential energy contribu- 
tions are not treated on equal footing. Indeed this is the 
essential difference between cluster mean-field approxi- 
mations and finite size calculations. In the former the 
potential energy of the lattice is truncated to that of the 
cluster whereas the kinetic energy is not. This leads to a 
self-consistent theory, generally (but not always) with a 
single-particle coupling between the cluster and the host. 
In the latter both the kinetic and potential energies of 
the lattice are truncated to their cluster counterparts. 
Therefore we might expect cluster methods to converge 
more quickly as a function of cluster size, compared to 
finite size techniques, for metallic systems with extended 
states and significant screening. That this is indeed the 
case was illustrated in the previous section for classical 
spin systems (see Fig.^l. 

For completeness we note that non-local interaction 
terms, e.g. a nearest neighbor Coulomb repulsion IJij can 
be treated in a similar manner by splitting it into intra- 
(t/c) and inter-cluster ((5t/) parts. A similar locator ex- 
pansion to Eq. H3U|I or respectively Eq. H4U|) is then writ- 
ten down with respect to the inter-cluster part W for 
the corresponding susceptibility (in this case this would 
be the charge susceptibility). As a result, an additional 
coarse-grained interaction Uc acts within the cluster (see 
also Sec. III. B() . and the locator expansion in bU leads to 
an additional self-consistency on the two-particle level. 
For a cluster size of iVc = 1 this formalism corresponds 
to the extended DMFT l|Smith and Sll2000() . A first ap- 
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plication of this extended cluster al gorithm to th e 2D t-J 
model for A^^, > 1 was discussed in iMaied ll2003|) . 



B. Cluster approximation to the grand potential 

In this section we provide a microscopic derivation of 
the CPT, the CDMFT and the DCA formalisms based 
on different cluster approximations to the diagrammatic 
expression for the grand potential. The advantage of 
this approach is that it allows us to employ almost all of 
the diagrammatic technology which has been developed 
in the past several decades to a new set of cluster for- 
malisms. Furthermore we are able to asses the quality of 
cluster approximations regarding their thermodynamic 
properties. The most significant disadvantage is that the 
formalism developed here only applies to systems which 
are amenable to a diagrammatic expansion. 

The following ideas will be illustrated on the extended 
Hubbard model Eq. H18|l . We use the notation introduced 
in Sec. III.AI Fig. O i.e. the cluster centers are denoted 
by X and sites within the cluster by X. The wave- vectors 
k and K are t heir respe ctive c onjugates. 

iBavm and~Ka danofli l)l96l|^ (see also iBavmL Il962l) 
showed that thermodynamically consistent approxima- 
tions may be constructed by requiring that the single- 
particle self-energy S fulfills 



G 



(5$[G] 
SG 



(44) 



i.e. is obtained as a functional derivative of the Baym- 
Kadanoff ^-functional with respect to the Green function 
G and that the approximation is self-consistent (via the 
left hand identity) . The Baym-Kadanoff generating func- 
tional <i>[G,U] is diagrammatically defined as a skeletal 
graph sum over all distinct compact closed connected di- 
agrams constructed from the Green function G and the 
interaction U. Thus, the diagrammatic form of the ap- 
proximate generating functional together with an appro- 
priate set of Dyson and Bethe-Salpeter equations, com- 
pletely defines the diagrammatic forma lism. 

As described in standard textbooks l|Abrikosov et all 
Il963|) the relation between the grand potential functional 
f2 and the ^-functional is expressed in terms of the linked 
cluster expansion as 

n[G,V] = -fcsr{$[G,U] -Trln(-G) -Tr(SG)} , 

(45) 

where the trace indicates summation over cluster sites X, 
superlattice wave-vectors k, frequency and spin. With 
the condition H44|) , the grand potential is stationary with 
respect to G, i.e. Sfl/5G — 0. Such approximations 
are thermodynamically consistent, i.e. observables cal- 
culated from the Green function G agree with those cal- 
cu lated as deriv atives of the grand potential fi. As shown 
bv lBavmll)l962|) the requirement (|44ll together with mo- 
mentum and energy conservation at the vertices also as- 
sures that the approximation preserves Ward identities, 
i.e. satisfies conservation laws. 



Prominent examples of conserving approximations in- 
clude the Hartree-Foc k theory and the fl uctuation ex- 
change approximation ijBickers et al\ . ll98^ . As exempli- 
fied by these theories, the typical approach to construct 
a conserving approximation is to restrict the diagrams in 
$ to a certain sub-class, usually the lowest-order (in the 
interaction U) diagrams. The resulting weak-coupling 
approximations however usually fail for systems where 
the interaction U is of the same order or larger than the 
bandwidth. 

Quantum cluster approaches go a different route: In- 
stead of neglecting classes of diagrams in $, quantum 
cluster approaches reduce the infinite number of degrees 
of freedom over which $ is evaluated to those of a fi- 
nite size cluster. In contrast to perturbative approaches 
however, all classes of diagrams are kept. 



1. Cluster perturbation theory 

The simplest way to reduce the degrees of freedom in 
$[G] is to replace the full lattice Green function G(k, z) 
by the Green function g{z) = [{z + /i)l — tc — Sc(z)]~^ 
of an isolated cluster of size Nc- Consequently, the self- 
energy Sc = (5$[g]/(5g obtained from $ is the self-energy 
of an isolated finite size cluster. This however leads to a 
theory which lacks self-consistency. Moreover one has to 
make the ad-hoc assumption that the lattice self-energy 
is identical to the one obtained from the cluster, Sc- The 
left hand side of Eq. then yields the form for the CPT 
lattice Green function 



G(k,z) = [Goi(k,z)) 



(46) 



where all the quantities are Nc x Nc matrices in the clus- 
ter sites. Since the bare lattice Green function is given 
by Go(k, z) = [{z + fi)t — t(k)]~^ and the hopping can 
be split into intra- and inter-cluster parts (see Eq. H25(l ^. 
t(k) = tc -I- (5t(k), we obtain 



G(k,z) = [g-i(2)-5t(k)]- 



(47) 



with g^^(z) = (z + /i)l — tc — Sc(z). This form was de- 
rived in Sec. III. A1 from the locator expansion Eq. (|30|l by 
ignoring the hybridization T between cluster and host. 
According to Eq. (|46|l . the CPT can be viewed as the ap- 
proximation that is obtained by replacing the self-energy 
in the Dyson equation of the lattice Green function G by 
the self-energy of an isolated cluster S^. T his idea was 
first developed bv iGros and Valeiitil l)l994() and applied 
to the 3-band Hubbard mode l. A different a pproach to 
derive the CPT was taken bvlPairault et al\ jl998,) (see 
also iPairault et all \200(l ISenechal et aJl 120021) . Thev 
showed that Eq. (|47|l is obtained as the leading order 
term in a strong coupling expansion in the hopping St 
between sites on different clusters. This derivation of the 
CPT provided a fundamental theoretical basis to assess 
the nature of the approximation as well as to systemat- 
ically improve the quality of the approach by including 
higher order terms in the perturbative expansion. 
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The CPT becomes exact in the weak couphng hmit 
U/t = and the strong-couphng hmit t/ U = as weh as 
in th e infinite cluster size limit iVc — > oo (|Senechal et all 
l2002j) . The limit t/U = is reproduced exactly since the 
CPT is the pcrturbative result in the hopping. In this 
limit, all the sites in the lattice are decoupled, and the 
system is solved exactly by the single-site Green func- 
tion gc- In the opposite limit U/t = 0, the cluster self- 
energy Sc in Eq. (|46|l vanishes and G(k) = Go(k) is 
the exact solution. In the limit Nc — > oo, the clus- 
ter Green function gc becomes the exact Green func- 
tion of the full system. At finite t/U and cluster size 
Nc = 1, the CPT re covers the Hubbard-I approxima- 
tion llHubbardL Il96l where the self-energy is approx- 
imated by the self- energy U'^/4uj (at half-fil ling) of an 
isolated at om l|Gros and Valentil 119941: ISenechal et aii . 

According to the derivation of the CPT, the cluster 
Green function gc is to be calculated on a cluster with 
open boundary conditions. Since the hopping tc between 
sites inside the cluster is treated exactly whereas the 
inter-cluster hopping St between surface sites on different 
clusters is treated perturbatively, translational invariance 
for sites X in the cluster is violated while it is preserved 
for sites x in the superlattice. As a result, the cluster 
wave-vector K is not a good quantum number and we 
have as a generalization of the Fourier transform Eq. H24(l 
(we omit the frequency dependence for convenience) 

^ Q X,,Xj 

X S(k - k' - Q) . (48) 

where k and k' are wave- vectors in the full Brillouin zone 
and Q is a wave-vector in the cluster reciprocal space. 
Here we used the relation G(k) ~ G(k) which follows 
from Eq. (|22l) by replacing k by k = k -|- K. To restore 
translational invariance in the full lattice Green function, 
the CPT approximates G(k, k') by the Q = contribu- 
tion to obtain 

GcPT{k) = ^ J2 e-^''(^--^^)G(X,,X„k), (49) 

as the translational invariant propagator used to calcu- 
late spectra. With this approximation, the CPT provides 
a very economical method to calculate the lattice Green 
function of an infinite size [N — *■ oo) Hubbard-like model 
from the Green function (or equivalently self-energy) of 
an isolated cluster of finite size Nc ^ N. From GcPT(k) 
one can calculate single-particle quantities such as pho- 
toemission spectra, kinetic and potential energies, double 
occupancy, etc. 

To reduce the numerical cost, it was suggested to use 
periodic boundary conditions on the CPT cluster by 
adding the appropriate hopping terms to the intra-cluster 
hopping tc and subtracting them from the inter-cluster 



hopping 5t ("Pahn ken et all l2002l) . However, as dis- 
cussed bv .Senechal et all l|2002|) . periodic boundary con- 
ditions lead to less accurate spectra for the ID Hubbard 
model than open boundary conditions. This a-posteriori 
argument for open boundary conditions is substantiated 
by calculations within Potthoff's self-energy functional 
approach (see Sec. III.H|I which show that the grand po- 
tential of the system is o nly stationary in the limit of 
open boundary conditions ijPotthoff et a/.l 12003(1 . 



2. Cellular dynamical mean-field theory 

A superior approximation may be obtained if, instead 
of the isolated cluster Green function g, the full lattice 
Green function G restricted to cluster sites is used to 
evaluate the functional This approximation can be 
motivated microscopically by approximating the momen- 
tum conservation on internal vertices in the diagrams 
defining Momentum conservation at each vertex is 
described by the Laue function 

(50) 

where ki, k2 (k j , k9) are the mornenta e ntering (leav- 
ing) the vertex. iMiiller-HartmannI (|l989b) showed that 
the DMFT may be derived by completely ignoring mo- 
mentum conservation at each internal vertex by setting 
A = 1. Then one may freely sum over all of the in- 
ternal momentum labels, and the Green functions in 
the diagrams are replaced by the local Green function 
G.. = l/iVEkG(k). 

The CDMFT and DCA (see below) techniques may 
also be defined by their respective approximations to the 
Laue function. In the CDMFT the Laue function is ap- 
proximated by 

AcDMFT = ^e*^-(^i+^i+^=^+^=^ + -'^'i-'''i-^^-ki — ■) 
X 

(51) 

Thus the CDMFT omits the phase factors e^^'^ resulting 
from the position of the cluster in the original lattice, 
but keeps the phase factors e^^'-^. The latter are directly 
responsible for the violation of translational invariance. 
Consequently, all quantities in the CDMFT are functions 
of two cluster momenta Ki, K2 or two sites Xi, X2 
respectively. 

If the CDMFT Laue function Eq. IHl is applied to 
diagrams in <I>, each Green function leg is replaced by the 
CDMFT coarse-grained Green function (the frequency 
dependence is dropped for notational convenience) 

G(Xi,X2) =G(Xi,X2;i-0) = 

^.(K,+k,)-x.G(Ki,K2;ki,k2)e-(K^+'^^)-^^ = 

K1.K2 

ki,k2 
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^ ^ C 



J2 G(Xi,X2,ki,k2), 



(52) 



ki,k2 



or in matrix notation for the cluster sites Xi and X2 



N 



(53) 



since G is diagonal in ki , k2 due to the translational in- 
variance of the superlattice. Similarly each interaction 
line is replaced by its coarse-grained result 



N 



(54) 



The summations over the cluster sites X within each 
diagram remain to be performed. As a consequence of 
coarse-graining the propagators in <i>, the CDMFT self- 
energy 



5G{z) 



(55) 



is restricted to cluster sites and consequently independent 
of k. Note that by definition, G and tj are truncated out- 
side the cluster, i.e. if the interaction U is non-local, U 
includes only interactions within, but not between clus- 
ters. 

The CDMFT estimate of the lattice grand potential is 
obtained by substituting the CDMFT approximate gen- 
erating functional <i>[G,U] into Eq. (|45|l . From the con- 
dition that the grand potential is stationary with respect 
to the lattice Green function, S^l/SG — 0, one obtains 
a relation between the lattice self-energy and the cluster 
self-energy 

E(Ki,K2;ki,k2) = 

^ g-»(Ki+kO-Xij.^(Xi,X2)e'(^^+''^)-^=. (56) 

Xi ,X2 

With Eq. the left hand side of Eq. then becomes 
the coarse-graining relation 

= ^Et^o '(k,z) - S,(z)]-i (57) 



with the bare Green function Go(z) = [(z+/i)l — t(k)] ^ . 



3. Dynamical cluster approximation 

In the DCA the phase factors e'*' -'^ are omitted too, 
so that the DCA approximation to the Laue function 
becomes 



^DCA — A^c'^Ki+K2 + ---,Ki+K^ + - 



(58) 



and Green function legs in $ are replaced by the DCA 
coarse grained Green function 



(59) 



since Green functions can be freely summed over the 
N/Nc wave-vectors k of the superlattice. Similarly, the 
interactions are replaced by the DCA coarse grained in- 
teraction 



^(K) = §E^(K + k). 



N 



(60) 



As with the CDMFT, the effect of coarse-graining the 
interaction is to reduce the effect of non-local interactions 
to within the cluster. This collapse of the diagrams in the 
<f> functional onto those of an effective cluster problem is 
illustrated in Fig. ^ for a second order contribution. 



K ■ q 





Ki+Ka,Ks + Kj 



U{Q) 



K' 




K'-l-Q 



Figure 4 A second-order term in the generating functional 
$ of the Hubbard model. The dashed line represents the 
interaction U, and the solid line on the left hand side (right 
hand side) the lattice (coarse-grained) single-particle Green 
function G (G). When the DCA Laue function is used to 
describe momentum conservation at the internal vertices, the 
wave- vectors collapse onto those of the cluster and each lattice 
Green function is replaced by its coarse-grained average. 

The resulting compact graphs are functionals of the 
coarse grained Green function G(K) and interaction 
C/(K), and thus depend on the cluster momenta K only. 
For example, when Nc — 1, only the local part of the 
interaction survives the coarse graining. As with the 
CDMFT, within the DCA it is important that both the 
interaction and the Green function are coarse-grained 
(Hcttler et ai, 2000^). As a consequence of the collapse 
of the <l>-diagrams, the DCA self-energy 



Sc(K,z) 



SG(K,z) 



(61) 



only depends on the cluster momenta K. 

To obtain the DCA estimate of the lattice grand po- 
tential, we substitute the DCA approximate generating 
functional ^G{K),U{K)] into Eq. The grand po- 

tential is stationary with respect to G when 



6G{k) 



Ee(K)-E(k) =0, (62) 
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which means that S(k) = Ec(K) is the proper ap- 
proximation for the lattice self-energy corresponding to 
$[G'(K),J7]. The self-consistency condition on the left 
hand side of Eq. (|44|l then becomes the coarse-graining 
relation 

^(K, z) = ^ ^[G„-i(K + k) - S,(K, . (63) 

k 

with the bare Green function Go (K + k, z) = [z — e^+k + 

Both the CDMFT and the DCA have well defined 
limits. In the infinite size cluster limit oo, the 

CDMFT and DCA approximations to the Laue function 
recover the exact Laue function, $ is evaluated with the 
full lattice Green function and interactions and thus the 
exact result is recovered. When A'c = 1, both the Laue 
functions reduce to A = 1, $ is evaluated with the local 
Green function and local contributions of the interac- 
tions, and the DMFT result is recovered. 

4. Self-consistency scheme 

The two equations H44|l form a non-linear set of equa- 
tions which have to be solved self-consistently to deter- 
mine the cluster self-energy Sc[G] with the use of a suit- 
able cluster solver. For cluster solvers that sum up all 
diagrams of Sc, i-e. in contrast to a skeletal expansion of 
Sc, an additional step is necessary in the self-consistent 
cycle. In order to not overcount self-energy diagrams, Sc 
is to be calculated as a functional of the corresponding 
bare propagator to G, the cluster excluded Green func- 
tion 

a(z) = [G-i(z)-f (64) 

This equation^ unambiguously defines the self-consistent 
iteration procedure illustrated in Fig. [S] 

1. The iteration is started by guessing an initial clus- 
ter self-energy Hdz), usually zero or the result from 
second order perturbation theory, to 

2. Calculate the coarse grained quantities 

= ^E[Go"'(k,z)-S,(z)]-i , U = ^^U(k) 

ic k 

3. The effective cluster problem is then set up with 
the cluster excluded Green function G{z) and U. 



^ A unifying matrix notation is used. In the CDMFT, the 
quantities are matrices in the N^. cluster sites and in partic- 
ular [Go^(k,2)]x,x,- = (z-Hm)5x,x, - [t{k)]x,x,. For the 
DCA, the matrices are diagonal in the cluster momenta K and 



4. The self-energy Sc(-z) or respectively the cluster 
Green function Gc{z) is calculated in the effective 
cluster model (see Sec. III.D|I by using any of the 
quantum cluster solvers discussed in Sec. IIIII 

5. For techniques that produce the cluster Green func- 
tion Gc rather than the self-energy, the new es- 
timate of the cluster self-energy is calculated as 

i:ciz)=g-Hz)-G,-\z). 

The iteration closes by re-calculating the coarse-grained 
Green function G(z) in step (2) with the new estimate 
of the cluster self-energy. This procedure is repeated un- 
til the cluster Green function Gc{z) equals the coarse- 
grained Green function G{z) to within the desired accu- 
racy. 



Cluster Solver 




g Gc 
t I 

\ / 

G-f Ek[Go'(k)-Sj-i 

Figure 5 Sketch of the CDMFT and DCA embedded cluster 
algorithms. The iteration starts with computing the coarse 
grained Green function G using an initial guess for the cluster 
self-energy Sc. The cluster excluded Green function Q is then 
used to define the effective cluster problem which yields a new 
estimate of Sc. 



C. Discussion 

In contrast to the DMFT, a unique setup for the em- 
bedded cluster theory does not exist. Depending on the 
trea tment of e.g. bou ndary conditions (see Sec. III.AI 
and lBiroli and Kotliail 12002) differences in the coupling 
of the cluster to its environment arise (see comparison 
in Sec. III. C. 41 below). In fact, there exist infinitely 
many realizations of embed ded cluster theories for 
any given mode l Hamiltonian ("Biroli and Kotliar, '200^1 
Okamoto efall. 12003: Potthofi. 20iI3b.: Potthoff et al. . 
20031) . two of which we focus on in this review. 

The fundamental approximation common to all ap- 
proaches is that they try to go beyond conven- 
tional mean-field approximations and introduce non-local 
physics by replacing the unsolvable lattice Hamiltonian 
by some manageable finite portion - possibly with ef- 
fective model parameters - and reintroduce the thermo- 
dynamic limit by a mean-field type treatment of the re- 
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maining system. Thus, the influence of truly long-ranged 
correlations, i.e. those beyond the cluster size, is still not 
incorporated, but short-ranged correlations and in par- 
ticular the local dynamics with respect to the cluster can 
ideally be treated exactly. That this can already lead 
to substantial renormalizations has been demonstrated 
in Sec. II. Bl and is also well-known from the DMFT. 
The combination of short- to medium-ranged correlations 
with mean-field treatment of long-ranged physics enables 
the investigation of a system's tendency to certain types 
of order not accessible by conventional mean-field theory. 
One can thus more clearly identify the possible existence 
of long-ranged correlations which are normally hard to 
establish in conventional finite-system calculations. On 
the other hand, it is this mean-field likeness that prac- 
tically disables a proof for the appearance of a phase 
transition in the real model, even though the behavior 
of transition temperatures etc. with cluster size can give 
valuable hints about the system (see Sec. lI.B|l . 



1. Conservation and thermodynamic consistency 

Both the DCA and CDMFT approximations are self- 
consistent and are ^-derivable since they satisfy Eq. (|44|l . 
Thus, they are thermodynamically consistent in the 
Baym-Kadanoff sense. Observables calculated from the 
lattice Green function G agree with those calculated as 
derivatives of the lattice grand potential fl. Since the 
CPT is not self-consistent, it is not thermodynamically 
consistent. However, none of these approaches is con- 
serving in the Baym-Kadanoff sense since they all vio- 
late momentum conservation. Thus, each of these ap- 
proa ches is likely to viol ate some set of the Ward identi- 
ties ()Hettler et all l200(]() . 



2. Causality 

One problem in any formation of an embedded cluster 
theory is the manifestation of causality, i.e. a physical 
Green function cannot have poles anywhere except on the 
real axis. In particular, for the fundamental quantity of 
the theory, the single particle Green function, this means 
that the proper self-energy in momentum space has to 
obey 9mE(k, u; -I- jO"*") < 0. Ea rly attempts to formu- 
late cluster extensions to DMFT l)Schiller and Ineersend . 
^95) ran into exactly this problem, which e.g. mani- 
fests itself in negative single-particle spectral functions. 
Explicit proof s for c ausality can be give n for the DCA 
llHettler et all \200(t\ and the CDMFT l)Kothar et all 
1200 It see also iBiroli et all \200^ . More precisely, any 
embedded cluster theory consistent with the locator ex- 
pansion (or a suitably defined cavity construction) obeys, 
due to Eqs. H28() and H32() . causality. A closer look at Eqs. 
(|28|l and H32|) also reveals how problems can arise. It is, 
for example, tempting to replace the well defined cluster 
quantity ^^(z) in H28|) by some approximation to S(z), 



i.e. the full self-energy. As ha s been discussed in some de- 
tail bv lOkamoto et al\ l)2003() (see also Sec. III.H.2'|) . such 
a procedure will in general introduce ringing phenomena 
and lead to acausal behavior. How strongly such a viola- 
tion of causality will eventually infiuence the interesting 
low-energy results is quite likely a question of the model 
parameters under investigation. In any case, it must be 
taken as serious reason to at least doubt the quantitative 
accuracy of such results. 



3. Reducible and irreducible quantities 

Fundamental quantities like the one-particle self- 
energy S(2), or its many-particle counterparts entering 
e.g. susceptibilities, carry the whole fragile information 
about the many-body physics of a given model. In the 
language of diagrammatic perturbation theory they are 
built of so-called irreducible diagrams. Hence they are 
also frequently called irreducible quantities; in contrast, 
the single-particle Green function or a susceptibility is a 
reducible quantity. It is an important aspect, that the 
cluster theories discussed use approximants for these ir- 
reducible quantities only, i.e. the quantities obtained as 
derivatives of the Baym-Kadanoff $ functional. In fact, 
in the formulation in Sec. III.BI any attempt to replace 
reducible quantities like the one-particle Green function 
directly by approximants is at least dangerous. 

To see this, consider the grand potential functional 
Eq. (|45|) . It is expressed as a sum over all closed 
connected distinct graphs constructed from the Green 
function G and interaction U. The subset of compact 
graphs comprise the Baym-Kadanoff generating func- 
tional <I>[G,U] which is expressed as a skeletal graph 
sum over all distinct compact closed connected graphs. 
Compact diagrams have no internal parts representing 
corrections to the Green function G. 





Figure 6 Non-compact (left) and compact (right) non-local 
corrections to the grand potential. Here the upper (lower) 
circle is meant to represent a set of graphs which are closed 
except for the external lines shown, and restricted to site n 
(the origin). 

In quantum cluster theories the graphs for are 
approximated by their cluster counterparts. As an 
example, consider the limit of infinite dimensions, 
used by Metzner and V oUhardt to derive the DMFT 
llMetzner and VoUhardti 11989.) . In this limit most closed 



connected graphs are local since the Green function falls 
off quickly with distance r, G{r) ~ D^^/^. In fact, only 
a small set of graphs, corresponding to non-compact cor- 
rections, remain finite. To see this, consider the simplest 
non-local corrections to non-compact and compact parts 
of the grand potential of a Hubbard-like model, illus- 
trated in Fig. Here the upper (lower) circle is a set of 
graphs composed of intra-site propagators restricted to 
site n (the origin). Consider all such non-local corrections 
on the shell of sites which are n mutually orthogonal unit 
translations from the origin. In the limit of high dimen- 
sions, there are 2"Z3!/((Z3 — n)l n\) ^ 0(D"') such sites . 
Since as D^oo, G{r) - £)-'^/2|Met zner and VollhardtL 
^89), the legs on the compact correction contribute a 
factor 0{D^^'^) whereas those on the non-compact cor- 
rection contribute 0(Z3^"). Therefore the compact non- 
local correction falls as D^" and vanishes as 13 — > cx); 
whereas, the non-compact correction remains of order 
one, regardless of how far site n is from the origin. As 
a result, the generating functional, which is composed of 
only compact graphs, is a functional of the local Green 
function and interaction in this limit 

$[G,U] = $[G„;, 0(l/i?). (65) 

A similar analysis was done for the DCA cluster problem 
byErvanoour et al. (2002). They find that the correc- 
tions from sites outside the cluster associated with com- 
pact diagrams are quite small (i.e. high order in the linear 
cluster size l/^c) justifying the approximation 

1>[G,U]«<i>[G(K),c7(K)], (66) 

while those associated with non-compact diagrams are 
large and cannot be neglected. The same analysis may 
be done for the CDMFT, simply by replacing the DCA 
graphs by those for the CDMFT. 

Thus, the essential approximation for the DMFT, the 
DCA and the CDMFT is to approximate the lattice gen- 
erating functional by its cluster counterpart in the 
estimate of the lattice grand potential, Eq. H45|l . Con- 
committantly the derivatives of <&, i.e. the lattice self- 
energy and vertex functions are approximated by their 
respective cluster counterparts. This once more under- 
lines why in embedded cluster theories it is necessary to 
always work with the cluster irreducible quantities; they 
are the only quantities compatible with the cluster and 
any attempt to include more information by e.g. Fourier 
transformation will lead to consistency problems, which 
eventually express themselves as causality violation. Irre- 
ducible quantities are also those important to discuss the 
convergence behavior of a cluster method with increasing 
cluster size as discussed in the next section. 



4. Comparison 

Detailed comparisons of t he DCA and CDMFT al - 
gorithms were presented in iMaier and Jarrelll l)2002() . 
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Maier et al\ ll2002bl) . iBiroli and Kotliail l)2002|) and 
Biroh et al\ l^2003^ . Both approximations share the un- 
derlying idea and general algorithm, and differ only in the 
form used for the hopping matrix t (see Eq. (|35f) ). The 
purpose of this section is to convey the consequences of 
this, at first sight, small difference on the effective clus- 
ter problem, convergence properties and the calculation 
of the lattice self-energy. 



a. Nature of effective cluster problem. For notational con- 
venience, we use a ID model with nearest neighbor hop- 
ping t only, set the on-site energy eo = and denote the 
cluster size by Lc- The generalization to higher dimen- 
sion or longer-ranged hopping terms is straightforward. 

The CDMFT uses the original form for the hopping 
matrix t(fc) which is obtained e.g. as an inter-cluster 
Fourier transform (see Eq. H22II ') of t(i).Only entries be- 
tween neighboring sites inside the cluster [t{k)]xi,Xi±i = 
—t and between neighboring sites on the surface of the 
cluster ['t{k)]xi,Xi±{L^-i) = —te^^'^^" are finite. The for- 
mer entries form the intra-cluster hopping matrix tc = 
Nc/NJ2k'^W while the latter entries form the inter- 
cluster hopping matrix 5t{k) = t(k) — tc- Both ampli- 
tudes are given by the original hopping t. For the effec- 
tive cluster problem this translates to the fact that only 
sites on the surface of the cluster couple to the effective 
medium, while sites inside the cluster only couple to their 
neighboring sites in the cluster. Hence the cluster prob- 
lem has open boundary conditions and translational in- 
variance is violated within the cluster. The lattice Green 
function G{k) (see Eq. (I30|l ) is a matrix in the cluster 
sites and cannot be diagonalized by going over to cluster 
if-space. Therefore the coarse-graining step Eq. is 
done in real space. 

The DCA restores translational invariance by setting 
[tDCACk)]x^x, = [t(fc)]x.x,e-'^-(^'-^^) (see Eq. 
As a consequence, its matrix elements become identi- 
cal and are given by — ie^*'^ between sites Xj = AT^ =F 1 
and Xj — Xi ± (Lc — 1). Hence the DCA hopping ma- 
trix tjjcAik) is fully cyclic with respect to the cluster 
sites and the lattice Green function G(fc) is diagonal- 
ized by going over to cluster if -space. The DCA intra- 
cluster hopping matrix tc.DCA = Nc/N J2k'^DCA{k) is 
also cyclic with finite matrix-elements 

ltc,DCA\x,.x, = -t — sin— (67) 
n Lc 

between sites Xj = Xi±l and Xj = Xi zL {Lc — 1) and 
the DCA cluster problem therefore has periodic bound- 
ary conditions. At finite cluster size Lc, the intra-cluster 
hopping is reduced by the factor l/6(7r/Lc)^+C((7r/Lc)^) 
compared to its lattice counterpart t. In the infinite clus- 
ter size limit it becomes t. This reduction in the intra- 
cluster coupling is compensated by the inter-cluster hop- 
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ping which is of long-ranged nature, 
r., sin[(i =F l)7i-/-^c 



sin(7r/£c) 
tt/Lc 



(68) 

between sites Xj — Xi^l and Xj = Xi± {Lc — !)■ It is 
important to note that 6tocA couples all the sites in the 
cluster to sites in the effective medium. It vanishes for 
5; = and decreases as 1/x with increasing x. We also 
notice that 



St 



DC A 



(69) 



for large linear cluster sizes Lc and emphasize that this 
result holds generally in any dimension D. 

The restoration of translational invariance in the DCA 
is achieved by mapping the lattice onto a cluster with 
periodic boundary conditions with reduced hopping in- 
tegrals tc^DCA and coupling every site in the cluster to 
a neighboring site in the effective medium through long- 
ranged hopping integrals dtBCAix). The sum of all fi- 
nite intra- and inter-cluster couplings for a given site is 
again given by the original value t. Similar conclusions 
about the nature of the effective CDMFT and DCA clus- 
ter problems were reached in a study of the large U limit 
of the Fahcov-KimbaU m odel (FKM see Eg . ifTS^ I. i.e. 
the classical Ising model llBiroli et aLl . l2003j) . 

We stress that clusters with linear size — 2 are 
special. Here both terms in the intra-cluster hop- 
ping Eq. 1671) give a contribution to the same matrix- 
element. Hence the nearest-neighbor hopping is given by 
—2t (Lc/tt) sin (tt/Lc), i-e. with the prefactor —2t instead 
of —t for larger clusters. This reflects the fact that every 
site sees its nearest neighbor twice due to the periodic 
boundary conditions. Non-local fluctuations are thus en- 
hanced in clusters with linear size Lc = 2 as seen e.g. in 
an over-proportional suppression of transition tempera- 
tures (see Secs. ll\rBl and lIV.D.2t . 



b. Convergence with cluster size. The differences in 
boundary conditions translate directly to different 
asymptotic behaviors for large cluster sizes Nc- In lead- 
ing order the hybridization with the mean-field, T van- 
ishes like St"^ as the cluster size increases (see Eqs. 
and lO). In the CDMFT the magnitude of 6t is of or- 
der one for the sites on the surface of the cluster and zero 
otherwise. The average hybridization per cluster site in 
the CDMFT thus scales like 

^CDMFT — T^Tr [Tcdmft] ~ — ) , (70) 

where the trace runs over cluster sites and frequency. 



D-l\ 



sites 



This behavior is evident since only the 2D{L 
on the surface contribute to the sum and Nc ~ L^ . In 
the DCA, 5t = 0(1/ Lc) (see Eq. ^). The average 
hybridization of the DCA cluster to the effective medium 
hence scales faster to zero as 



^DCA 



1 



Tr [Tdca] 



(71) 



since each of the Nc terms contributes a term of the order 
Oil/Ll). 

It can further be shown that T acts as the small pa- 
rameter in these theories: The approximation performed 
by the DCA and the CDMFT is to replace the_lattice 
Green function G by its coarse-grained quantity G in di- 
agrams for the generating functional $ (see Sec. III.B|I . 
Both Green functions differ by the inter-cluster hop- 
ping (5t(k), the self-energy and the hybridization fmic- 
tion r. Since the diagrams in $ are summed over k, 
the terms having the same order as (5t(k) vanish since 
Nc/NJ2k^'^0^) = 0- If assume that the self-energy 
has corrections of the same or higher order in 1/Lc as 
r, the convergence of $ is entirely determined by T. 
With the scaling relations H7()|l and 171|l we find for the 
CDMFT that $cdmft « $ + 0{l/Lc) while the DCA 
converges like ^dca « ^ + 0{l/Ll). Since S = S^/SG, 
it converges with Lc as the corresponding <i> confirming 
the assumption. The generating functional $ and hence 
the grand potential thus converges faster in the DCA 
than in the CDMFT with increasing cluster size. 



c. Numerical comparison. The scaling behavior Eqs. H70() 
and (jZU of the CDMFT and DCA average hybridiza- 
ti on strengths was verified numeric ally in the ID FKM 
in lMaier and JarreJ (2002^ (see also' Maier et a^J . l2002b(l . 
Here we review the effects of the different scaling behav- 
iors of the average hybridization on the phase transition 
in this model. The Hamiltonian of the FKM is discussed 
in Sec. IIV.BI Eq. (|182|l . It can be considered as a sim- 
plified Hubbard model with only one spin-species being 
allowed to hop. However it still shows a complex phase 
diagram including a Mott gap for large U at half fill- 
ing, an Ising-like charge ordering with the corresponding 
transition temperature Tc being zero in ID, and phase 
separation in all dimensions. Since the ID FKM is in the 
ID Ising universality class we expect similar scaling be- 
havior as observed in the results for the ID Ising model 
in Sec. ll.BI In particular, we expect finite transition tem- 
peratures within both cluster approximations due to their 
residual mean-field character. The CDMFT and DCA ef- 
fective clust er models were solve d with a QMC approach 
described in lHettler eTaI\ l|200(]|) . 

The DCA transition temperature Tc was determined 
from the divergence of the lattice charge susceptibility 
xiQ) calculated from the particle-hole correlation func- 
tion as detailed in Sec. HTfI In the CDMFT formal- 
ism the calculation of lattice susceptibilities is difficult if 
not impossible due to the lack of translational invariance. 
Here Tc is determined from the calculation of the charge 
order parameter m as detailed in Sec. III. El For the DCA 
both techniques are illustrated in the inset to Fig. [3 

As for the ID Ising model (see Fig.|2in Sec. II. Bp . the 
DCA result for Tc scales to zero almost linearly in l/Nc 
for large Nc- Moreover, Tc obtained from DCA is smaller 
and thus closer to the exact result than the Tc obtained 
from CDMFT. The CDMFT does not show any scaling 
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Figure 7 Transition temperature versus inverse cluster size 
calculated with DCA/QMC (circles) and CDMFT/QMC 
(squares) when U — W — 41 = 1. Inset: DCA order pa- 
rameter m{T) and inverse charge susceptibility x{Q)~^ versus 
temperature. Taken from lMaier et all i2002bD . 



behavior and in fact seems to tend to a finite value for Tc 
as Nc — > oo. As explained above, this striking difference 
can be attributed to the different boundary conditions. 
The open boundary conditions of the CDMFT cluster 
result in a large surface contribution so that Tcdmft > 
^DCA- This engenders pronounced mean-field behavior 
that stabilizes the finite temperature transition for the 
cluster sizes treated in Fig. [7| For larger cluster sizes the 
bulk contribution to the CDMFT grand potential should 
dominate so that Tc is expected to fall to zero. 

Complementary results are found in simulations of 
finite size systems. In general, systems with open 
boundary conditions are expected to have a surface 
contribution in the gran d potential of order 0{1/Lc) 
l|Fisher and Barb'eil Il972ft . This term is absent in sys- 
tems with periodic boundary conditions. As a result, 
simulations of finite size systems with periodic boundary 
conditions converge much more quickly th an those with 
open boundary conditions (jLandaulll97(? ). 

The DCA converges faster than the CDMFT for criti- 
cal properties as well as extended cluster quantities due 
to the different boundary conditions and the coupling 
to the mean-field. As detailed above, each site in the 
DCA cluster experiences the same coupling to the effec- 
tive medium, while in the CDMFT only the sites on the 
boundary of the cluster couple to the mean- field host. 
Provided that the system is far from a transition, the 
sites in the center of the CDMFT cluster couple to the 
mean-field only through propagators which fall exponen- 
tially with distance. Local results such as the single- 
particle density of states thus converge more quickly 
in the CDMFT wh en measured on a central site (see 
iMaier et a/.l.l2002b(l . 



d. Calculation of the lattice self-energy. Another signifi- 
cant difference between the two cluster techniques ap- 



pears in the calculation on the lattice self-energy. The 
DCA approximates the lattice self-energy by a constant 
within a DCA cell in momentum space, I](K + k, z) = 
Sc(K,z). Therefore the self-energy is a step function 
in k-space. In order to obtain smooth non-local quan- 
tities such as the Fermi surface or the band structure, 
an interpolated S(k, z) may be used. Bilinear interpo- 
lation in 2D is guaranteed to preserve the sign of the 
function however leads to kinks in E(k, z). Yielding the 
smoothest possible interpolation of 5](K, z), the use of an 
Akima spline which does not overshoot is consistent with 
the DCA assumption that the self-energy is a smoothly 
varying function in k-space. However it is important to 
note that this interpolated self-energy should not be used 
in the self-consistent loop as this can lead to violations 
of causality as discussed above. 

In the CDMFT, the lattice self-energy is given by the 
Fourier transform of the cluster self-energy (see Eq. (|56|l 'l 

5](k,k',z) = E e-'''-'''S,(X„X„z)e*'-^^ 

^ Q X,,Xj 

X S{k - k' - Q) . (72) 

Since the CDMFT cluster violates translational invari- 
ance, the lattice self-energy depends on two momenta k 
and k' which can differ by a wave-vector Q of the clus- 
ter reciprocal space. To restore translational invariance, 
the CDMFT approximates the lattice self-energ y by the 
Q = contribution to give ijKotliar et all . l200lj) 



(73) 



Xi,X,- 



In real space, the lattice self-energy 

S(x,,Xj,z) = I]c(X,,Xj)^x.-x,,x,-x, (74) 



Xi.Xj 

is thus obtained by averaging over those cluster self- 
energy elements I](Xi,Xj) where the distance X^ — 
X,- equals the d i stance — Xj. As explained in 
iBiroli and Kotliail l)2002(l . the factor 1/Nc leads to an 
underestimation of non-local self-energy contributions at 
small cluster sizes, since the number of contributions 
for fixed x^ — Xj > in the sum Eq. H74|) is always 
smaller than Nc. As a po ssible solution to this problem, 
iBiroli and Kotliail l)2002|) suggested to replace the form 
(|74|l for the lattice self-energy by a weighted sum which 
preserves causality. One could e.g. weight the terms in 
the sum by their number instead of Nc to achieve better 
results. 

It is important to note that as in the DCA, the lat- 
tice self-energy Eg. (1751 or (ITHl d oes not enter the self- 
consistent loop. IBiroli et all l)2003j) however realized that 
a translational invariant formulation of the CDMFT algo- 
rithm can be obtained by replacing the cluster self-energy 
Sc(-z) by the translationally invariant lattice self-energy 
S(k), Eq. l(7^ in the coarse-graining step Eq. ifS^ . De- 
spite the dependence on k, this form of S(k ) can be 
shown to preserve causality ijBiroli et all l2003(l . 
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D. Effective cluster model 

Quantum cluster approaches reduce the complexity of 
the lattice problem with infinite degrees of freedom to 
a (numerically) solvable system with Nc degrees of free- 
dom. As detailed in Sec. III.BI this is achieved through 
the approximation of <i>[G, U], the exact Baym-Kadanoff 
functional of the exact Green function G and interaction 
U, by a spatially localized quantity <i>[G,U] which is a 
functional by the corresponding (coarse-grained) quanti- 
ties restricted to the cluster sites, G — Nc/NJ2i_G{'k) 

andU = iVe/A^EkU(k). 

$[G,U] may be calculated non-perturbatively as the 
solution of a quantum cluster model 

He = i/c.o + H,j . (75) 

Tic consists of a non-interacting term Hcfi describing the 
bare cluster degrees of freedom and their coupling to a 
host. The interacting term Hcj is related to the corre- 
sponding term in the original lattice model through the 
coarse-grained interaction U. This ensures that the func- 
tional dependencies of the cluster functional $c and its 
lattice counterpart <I> are identical. 

The non-interacting term iJc,o is fixed by the require- 
ment that the Green function Gc of the cluster model 
equals the coarse-grained Green function G of the origi- 
nal model 

Ge = G= [6?-i-S,]-i, (76) 

and hence is specified by the cluster-excluded Green func- 
tion g (see Eg. 

For Green function or respectively action based cluster 
solvers, like e.g. perturbation theory or QMC, Hc.o can 
hence be encoded in the cluster-excluded Green function 
Q. The corresponding cluster action Sc for the fermionic 
cluster degrees of freedom represented by the Grassman 
variables 7, 7* reads in imaginary time and cluster real 
space 

+ /^^ E ^7*.(r)7;.'(r)7,.'(T)7..(r) , (77) 

where we used the short hand notation i, j for the cluster 
sites Xj,_Xj. Note that for the CDMFT the quantities 
Qij and Uij are given by Eqs. and respectively, 
while for the DCA they are given by the cluster Fourier 
transforms of g(K) = [G-i(K,z) + Ec(K,z)]-i and re- 
spectively of f/(K) (see Eq. 

For Hamiltonian-based techniques, like e.g. the non- 
crossing approximation, exact diagonalization or numer- 
ical renormalization group, the need for an explicit for- 
mulation of Hc.o is inevitable. To setup the bare part 



Hcfl, it is convenient to use Eq. for the CDMFT or 
respectively (1^ for the DCA together with Eq. iTT^ to 
represent the cluster excluded Green function Q. In the 
CDMFT, we have with Eq. ^ 

g{z) = [{z + ^i)t-tc-T{z)]-\ (78) 

and the matrix-elements of the intra-cluster hopping tc 
are given by the hopping amplitudes of the original lat- 
tice, tij as detailed in Sec. III. C. 41 The non-interacting 
problem is thus split into two parts, cluster degrees of 
freedom with hopping integrals tij and their coupling to 
a dynamic host described by the hybridization function 
r(z). The CDMFT cluster mo del can hence be written 
as (see also iBolech et oil l2003^ 

nc = Hc,o + Hc,i = ^(ty - Ai<5»j)ctcj^ 

ij,a 

+ X..,~\a\^ a.,- + h.c] 

+ ^[y,,.(k)cLa^.^^ + /..c.] 

■ij,k,CT 

+ E ^4sl-S.'C,,. (79) 

ij,aa' 

The first part describes the isolated cluster degrees of 
freedom with fermionic creation (destruction) operators 
cL (Cia)- The second term simulates the host degrees 
of freedom as a non-interacting conduction band with 
the help of auxiliary operators a'- {a.r ) and energy 
levels X-J^^. The coupling between the cluster states and 

the bath with an amplitude VijCk) is described by the 
third term and the interacting term is given by the last 
term. The sums over k run over the N /Nc wave- vectors 
of the superlattice. Integrating out the auxiliary degrees 
of freedom yields an action of the form (|77|l with 

g,j{z) = [(z + /i)l-t,-r,(2)]r.i (80) 
= ^l^,r(k)[^l-A(k)]r„V™,(k). (81) 

^m,k 

Self-consistency then requires that the auxiliary param- 
eters Ayk and Vij (k) are chosen in a way such that the 
cluster hybridization function Tc{z) is identical to its lat- 
tice counterpart T{z) defined in Eq. H78|) . It is important 
to note that, while the isolated cluster parameter tc can 
be deduced directly from the non-interacting part of the 
lattice system, the energy levels Ak and coupling con- 
stants Vij{k) of the auxiliary particles are not known 
a-priori, but determined by the self-consistency condi- 
tion Tc = r. Since F^ is only finite on the surface of 
the cluster (see discussion in Sec. III.C^ . the coupfing 
Vij (k) between the cluster and the host is only finite for 
sites i on the surface of the cluster which effectively re- 
duces the number of baths. This was numerically verified 
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inCDMFT exact diagonalization studies bv lBolech et al\ 

il2nn,'j) . 

For the DCA we have with Eqs. HH) and (gSJ 

g{K,z)^[z-e-K + fi-T{K,z)]-\ (82) 

and hence the DCA effective ciuster modei is best repre- 
sented in cluster K-space as 

T^c = Hcfi + Hcj 

K,(T k,(T 

+ E [^K(k)cL«K+k. + ^-C-] 
K,k,(T 

+ ^K+Oa(^K'-Oa'^'K'<T'(''Ka ■ (83) 

Q 

The operators c^^ (cko-) create (destroy) an electron 
with momentum K and spin a on the cluster. t/(Q) is 
the Coulomb repulsion between electrons on the cluster 
defined in Eq. (|60|) and the sum over k in the coupling 
term again is restricted to the N/Nc wave- vectors of the 
superlattice. Analogous to the CDMFT case it is easy to 
see that the DCA cluster model yields an action of the 
form H77|l (in cluster Fourier space) with a. Q oi the form 
()82|1 and the cluster hybridization function 

The auxiliary parameters of the DCA cluster model are 
then determined by the condition rc(K) = r(K). 

For Nc = 1 both the CDMFT and the DCA cluster 
models reduce to the single-impurity Anderson model. 
If self-consistency is also established on the two-particle 
level (see discussion at the end of Sec. III.A|) via a sus- 
ceptibility, an additional coupling to a bosonic field rep- 
resenting the corresponding fluctuations has to be con - 
sidered in the cluster model (for details see lMaieill2003(l . 
For A^c = 1 the cluster model then reduces to the ef- 
fective impuri t y mo del used in the EDMFT approach 
llSmith and SilboOd) . 

As detailed in Sec. III.AI the CPT formalism sets the 
hybridization function T to zero, i.e considers an isolated 
cluster without the coupling to a host. Thus the CPT 
cluster model is identical to the original lattice model 
restricted to cluster sites, i.e. given by the first and last 
terms of the CDMFT cluster model Eq. (f75|l . 

E. Phases with broken symmetry 

For simplicity, in the preceding sections the self- 
consistent equations of quantum cluster theories have 
been derived assuming the absence of symmetry- 
breaking. In Sec. III.FI we review how instabilities to 



ordered phases can be identified by the computation of 
response functions. However, to explore the nature as 
well as possible coexistence of symmetry-broken states, 
generalizations of the cluster algorithms that explicitly 
account for the possibility of symmetry-breaking on the 
single-particle level are necessary. 

The applicability and modifications required to treat 
symmetry broken phases depend on the cluster approach. 
The CPT formalism is not amenable to the description 
of ordered phases because its self-energy is that of a 
finite isolated cluster in which spontaneous symmetry 
breaking cannot occur. However, Dahn ken et al. devel- 
oped a variational extension of the CPT ijDahnken et all 
ba sed on the self-energy functional approach by 
Potthofj l|2003bl) which yields a self-consistent scheme to 
study ordered phases (for details see Sec. III. Hp . The 
CDMFT formalism can describe ordered phases which 
are identifiable by a broken translational symmetry (such 
as antiferromagnetism) by construction, since the trans- 
lational sy mmetry of the CDMFT cluster is already 
broken (see iBiroli' g^l l2003t iMaier and Jarrell l2002t 
iMaier et aLLl2002h(l . Hence translational invariant solu- 
tions are often found to be unstable against the ordered 
one ijBiroh et all l2003f) . The DCA formalism is trans- 
lationally invariant by construction, and therefore gen- 
eralizations of the algorithm are necessary to treat or- 
dered phases. To keep this section concise, we exemplify 
the necessary generalizations of the DCA formalism to 
a selection of relevant types of symmetry-broken phases 
along with the mapping onto the corresponding cluster 
models. The adoption of the presented concepts to the 
CDMFT approach is straightforward. 

Once the equations are generalized to account for sym- 
metry breaking, the requisite algorithmic changes are rel- 
atively simple. A finite conjugate external field is used to 
initialize the calculation and break the symmetry. The 
field is then switched off after a few iterations and the 
system relaxes to its equilibrium state in the absence of 
external fields. On the other hand, if the field remains 
small and finite, the dependence of the order parameter 
on the field can be determined and used as an alternate 
way to calculate the susceptibility (by extrapolation to 
zero field). This approach is especially useful for clus- 
ter solvers such as the non-crossing approximation or the 
fluctuation-exchange approximation where the computa- 
tion of two-particle correlation functions is numerically 
too expensive. 



1. Uniform magnetic field - Ferromagnetism 

We flrst consider the formalism necessary to treat fer- 
romagnetism. A flnite homogeneous external magnetic 
fleld h is introduced which acts on the spin cr of the 
fermions according to the Zeeman term 

-hJ2<^clc,a- (85) 

i,cr 
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The effect of h on the motion of the spatial degree of 
freedom of the electrons, i.e. the diamagnetic term, is 
neglected for simplicity^. 

In the presence of finite h or a uniform magnetization, 
the single-particle Green functions for up- and down- 
electrons are not equivalent. To account for this SU{2) 
symmetry-breaking, the spin-index of the Green function 
Go-, self-energy Yi„ and effective medium Q„ (and hence 
Va) in the derivation of the DCA-equations has to be re- 
tained. For a finite uniform magnetic field h the DCA 
lattice Green function reads 

G.(k, z) = — — ^- ^ , (86) 

and the coarse grained and corresponding cluster- 
excluded Green function 



G^{K,z) = ^^G<,(K + k,2), 



N 



g;1(k,z)-he,,,(k,z) 



(87) 
(88) 



become spin-dependent. 

The action of the effective cluster model is identical 
to the action in the paramagnetic state, Eq. H77|l . but 
the spin indices have to be explicitly retained. It then 
describes electrons in an external magnetic field h cou- 
pled to a spin-dependent host and self-consistency is es- 
tablished by equating the Green function of the cluster 
model with the coarse-grained Green function (|87|) . 

In analogy, for Hamiltonian based cluster solvers, an 
additional term 



K,(T 



(89) 



is added to the cluster Hamiltonian, Eq. I|83|) , in the pres- 
ence of a finite external magnetic field h. The coarse- 
grained Green function G, Eq. (|86|l . is then used to cal- 
culate the magnetization to :— i/N'^^^a{nia-) according 



to 



TO 



^^(7a(K,T = o-) 



(90) 



and after analytical continuation 



Kcr 



dLuf{u)^mGa{K,uj + i6). (91) 



2. Superconductivity 

In this and the next section we generalize the DCA 
formalism to treat phases with superconducting and an- 
tiferromagnetic order, respectively. For better read- 
ability, we refrain from discussing the description of 



^ In 2D systems the magnetic field can be applied parallel to the 
plane to avoid orbital effects. 



phases with coexisting superconducting and antiferro- 
magnetic order. The extension to an integrated for- 
malism is straightforward and ha s been discussed in 
iLichtenstein and Ka tsnclson' ('2000) . 

We consider superconducting phases where the elec- 
trons are paired in spin singlet or triplet states with 
Sz = indicated by finite values of the order param- 
eter Ak := (cktc_kx) for some k. In addition to the 
normal Green function G(k, r) it is therefore necessary 
to introduce the anomalous Green function F(k, r) = 
— {TTC]i-\{T)c-i!_i). The spatial and temporal symmetry 
of the pairing state can then be inferred from the sym- 
metries of F. Since F describes the pairing of fermions, 
it necessarily is antisymmetric under the exchange of two 
particles. The spatial symmetry of the pairing state 
is determined by the k-dependence of the anomalous 
Green function F(k, r). If we assume conventional even- 
frequency pairing i^(k, — r) = i^(k, r), in the case of 
spin-singlet pairing, F has to be symmetric in k, i.e. 
F(— k, r) — i^(k, r) as is the case for even parity or- 
der parameters such as e.g. s-wave and d-wave. In the 
case of spin-triplet pairing F is antisymmetric in k. i.e. 
F(— k, r) = — F(k, r) as e.g. in a p-wave state. 

The allowed symmetry of the pairing state is re- 
stricted by the cluster geometry. It depends upon the k- 
dependence of the dispersion Ck and the K-dependence of 
the cluster self-energy S(K, r). When Nc — 1, is local 
and the k-dependence of F(k, r) is entirely through Ck- 
Hence only pairing states with the symmetry of the lat- 
tice such as s-wave and extended s - wave can be described 
by th is formalism ijjarrell Il992l: iJarrell and Pruschkel 
119931) ■ Larger cluster sizes are necessary to describe 
order-parameters with a symmetry less than the lattice 
symmetry, e.g. simulations with Nc — 4 are necessary 
to describe phases with a da;2_y2-wave order parameter 
which transforms according to cosfc^, — cosfcj,. 

By utilizing the concept of Nanibu-spinors 



(CkT'' 



(92) 



the self-consistent equations can be written in a more 
compact form, since the corresponding Green function 
matrix in Nambu space 



G(k,z) :=((vl'k;*l». 



G(k,z) F(k,z) 
F*(k,z*) -G*(-k,-z*) 



(93) 

contains information about both the normal and anoma- 
lous Green functions. In the presence of an external pair- 
ing field ?7(k) = ^^'(k) +ir]"{'k) which couples to c_kxCkT> 
the non-interacting part of the Hubbard Hamiltonian can 

be written as Hq = Ek ^kfekcrs - ?7'(k)cri + ?7"(k)cr2]*k 
so the the DCA lattice Green function in the supercon- 
ducting state becomes 

G(k, z) = [z(Jo~ {e\^~ y)(7z 

-77'(k)ai - 7,"(k)a2 - S,(K, z)]-^ , (94) 

where k = K -I- k and ai are the Pauli-spin matrices. 
The diagonal parts of the Nambu-matrix Sc(K, z) de- 
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scribe quasiparticle renormalizations and the ofF-diagonal 
parts contain information about the K- and frequency de- 
pendence of the pairing state. Again, the coarse-grained 
Green function 

N, 



(95) 



N 

k 

G{K,z) _F{K,z) 
F*{K,z*) -G*{-K,~z*) 

is used to calculate the cluster-excluded Green function 

g(K,z) = [G-i(K,z) + Se(K,z)]-i (96) 

which together with the coarse-grained interaction C/(K) 
defines the action of the corresponding effective cluster 
model 



Sc 



r')i>Ar') 



^ 2 



V'l(T)a3^,(T)l \{^\(T)o^i,AA{'^l) 



As in the normal state case (see Eq. ((TTJ), the cluster 
action is represented in cluster real space and all the 
quantities are cluster Fourier transforms of the corre- 
sponding quantities in cluster K-space. The spinors of 

Grassmann- variables V'l = '^i — (V'l)^ gen- 

erate coherent states corresponding to the cluster Fourier 
transform of the Nambu-spinors Eq. 192|) . In analogy, 
the corresponding cluster Hamiltonian in the supercon- 
ducting state is obtained from the normal state cluster 
Hamiltonian Eq. (|83|l by representing it with the Nambu- 
spinors H92() and adding a U(l) symmetry breaking term 



K 



X]*l,[77'(K)ai-r?"(K)a2]*K, 
iVc/iV ^j^?7(K -(- k) is the coarse-grained 



where fy(K) 
pair-field. 

After self-consistency is established by requiring that 
the Green function of the effective cluster model calcu- 
lated with the action Eq. H97I) equals the coarse-grained 
Green function Eq. H95(l , the order parameter Ak can be 
calculated. Within the DCA the resolution in k-space is 
restricted to the cluster K-points and the order parame- 
ter is coarse-grained 



3. Antiferromagnetic order 

In this section we derive the DCA cluster formalism for 
antiferromagnetism on a bipartite lattice. This formalism 
is appropriate when Nc > 1 . A formalism suitable for the 
case w hen A/'c = 1 is discussed in detail in lGeorges et al\ 
l)l996j) in the context of DMFT. 

The antiferromagnetically ordered state is character- 
ized by a spatial variation of the magnetization according 
to 



m(x) 



(101) 



where x denotes the sites in the lattice and Q is the 
antiferromagnetic wave- vector (Q = (tt, tt) in 2D). Hence 
bipartite lattices can be divided into two inequivalent 
sublattices according to 



1 : X e A sub-lattice 
-1 : X e B sub-lattice 



(102) 



The magnetic ordering thus reduces the translational 
symmetry of the original lattice. The volume of the mag- 
netic unit cell is doubled compared to the structural unit 
cell. Accordingly, the volume of the first Brillouin zone 
in the antiferromagnetic state is reduced to half the orig- 
inal volume and Q = (tt, tt) becomes a reciprocal lattice 
vector. Hence in neutron scattering experiments Bragg 
peaks are found at the wave-vector Q corresponding to 
a period of two lattice spacings. 

As a consequence of the translational symmetry break- 
ing, the correlation function 



G„(k,k + Q;r) = -(T,Ck,(T)4 



k+Qcr/ 



(103) 



becomes finite. Along the lines of the formalism for the 
superconducting state it is thus convenient to introduce 
spinors 



t 

kcr 



I' kcr' ''k+Q(T 



); * 



k(T 



(104) 



for the antiferromagnetic state. In the presence of a 
staggered external magnetic field /i(x) = /iexp(zQ • x) 
the non-interacting part of the Hamiltonian for bipar- 
tite lattices with ek+Q = ^^k then becomes = 

X^kCT ^k(T(^kf3 — f t'''i)^ko- where the prime on the sum 
indicates summation over the reduced Brillouin zone 
only. The corresponding Green function 



N ^ 



(cK+kTC-(K+k)i) = ^(K, T = 0+) (99) G.(k, z) = *L»^ 



(105) 



and given by the equal-time coarse-grained anomalous 
Green function F as (after analytical continuation) 

i?(K,r = 0+) = i^i^(K,iu;„) 



with k = K + k, is coarse-grained over the DCA cells 



/3 



C^-(K) = §5]G,(K + k) 



du) tanh 



/3w 



N 



(106) 



'^mF{K,uj + iS). (100) 



G,(K,K) G,(K,K + Q) \ 
G,(K + Q,K) G,(K + Q,K + Q) J ' 
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where we dropped the frequency argument for conve- 
nience. The cluster excluded Green function 

gA^,z) = [G;i(K,z) + Se,.(K,z)]-i. (107) 

then has two elements, Ga(K., K) and Ga(K., K + Q). As 
a result, its Fourier transform to real space 

- -l5:e'^-(^'-^^)e.(K) (108) 

^ K 

breaks translational symmetry. The action of the cor- 
responding cluster model in the antiferromagnetic state 
is then formally identical to the action in the paramag- 
netic state, Eq. (|77|l . As in the superconducting state the 
corresponding cluster model is obtained from the para- 
magnetic cluster model Eq. H83() by changing the repre- 
sentation to the Nambu spinors Eq. (|104|l and adding a 
symmetry breaking term 

Ka 

to account for the external staggered magnetic field h. 

After convergence, the sub-lattice magnetization m = 
i/N^^^ e'^^'^^(j{nia-) may be calculated from the off- 
diagonal component of the Green function matrix 
Eq. H106() according to 

TO = -^EaG,(K,K + Q;r = 0-) (110) 

Kcr 

4-00 

= --^^y / dijf{u;)a':smGaiK.K + Cl-Lj + i5). 

A/cTT j-^ J 

F. Calculation of susceptibilities 

A convenient way to identify continuous phase tran- 
sitions is to search for divergences of susceptibili- 
ties. One particular advantage of quantum cluster 
theories is that they allow us to consistently calcu- 
late these sus ceptibilities fr o m the corresponding cluster 
susceDtibilitv l|Hettler et aj] . 1200(1) . Unfortunately, the 
calculation of two-particle correlation functions in the 
CDMFT formalism is strongly hampered by the viola- 
tion of translational invariance on the cluster. So, in this 
section, we will restrict our attention to the calculation of 
two-p articl e quantities in th DC A following lHettler et al\ 
ll2000D and lJarrell fiF^ ll2001hD . 

As a specific example, we describe here the calculation 
of the two-particle Green function 

rP rP 

Xaa'{q,k,k') = / / / / dTidT2dT3dT4 

Jo Jo Jo Jo 

^ gi(((i)„+!^)Ti-W„T2+W„/T3-(w„/+l')T4) 

^ (7Vc^q^(Tl)cka(T2)cj,,^,(T3)ck'+q^'(T4)) , 



where we a d opt the conventional notation 
l|Abrikosov et all \l96^ k = (k,iw„), k' = (k,iw„/), 
q = (q, w„) and Tj- is the time ordering operator. 

X(ycr'{(l,k^k') and the irreducible two-particle vertex 
function T^a' (q, k, k') (not to be confused with the single- 
particle hybridization function) are related to each other 
through the Bethe-Salpeter equation 

X.a' {q, fc, k') - xL' (<?> k, k') + xL" (9, k, k") 

X Tcr"a"'{q,k",k"')Xa"'a'{q,k"',k') (111) 

where Xaa'il^^^^") non-interacting susceptibility 

constructed from a pair of fully-dressed single-particle 
Green functions. As usual, a summation is to be made 
for repeated indices. 

We now make the DCA substitution Fo-o-' (q, K + 
k,K' + k') -> Tc^^, (q,K,K') in Eq. (fTTT|l where Fc is 
the irreducible two-particle vertex calculated on the clus- 
ter (frequency labels have been suppressed). Note that 
only the bare and dressed two-particle Green functions 
X depend upon the superlattice wave-vectors k. Since x 
and x° in the product on the RHS of Eq. (|111|1 share no 
common momentum labels, we may freely sum over the 
wave-vectors k, yielding 

Xaa' {q, K, K') = xL' (9, K, K') + xL" K") 
X T,,„,,„{q,K",K"')x.'".'{q,K"',K'). (112) 

By coarse-graining the Bethe-Salpeter equation, we have 
greatly reduced its complexity; each of the matrices 
above is sufficiently small that they may be easily ma- 
nipulated using standard techniques. 

In contrast with the single-particle case where the 
coarse-grained quantities are identical to those of the 
cluster, the cluster quantity Xco-o-' (?) K') is not equal 
to x<y<y'{qjK,K'). This is because the self-consistency is 
established only at the single-particle level. Unlike the 
single-particle case where both Yjc{K) and G{K) are di- 
rectly calculated, neither Fco-o-' (?, ^4^, -?^') nor the coarse- 
grained susceptibility Xo-o-' (?, K, K') are calculated dur- 
ing the self-consistency. Instead, the coarse-grained non- 
interacting susceptibility x^o-' {q^ calculated in 
a separate program after the DCA converges using the 
relation 

X ^ E G^(K -f k, ia;„)G„(K + k + q, iw„ Un) (113) 

fe 

The corresponding cluster susceptibility is calculated by 
the cluster solver, e.g. the QMC process, as discussed in 
Sec. nil. (Ill and the vertex function is extracted by in- 
verting the cluster two-particle Bethe-Salpeter equation 

Xe..' (9, K, K') = XcL' iq. K, K') + xZ,', [q, K, K") 

xF,,„,,„ {q, K", K"')xc.,„,, [q, K"\ if 1114) 

If we combine Eqs. I|114() and H112(l . then the coarse- 
grained susceptibility may be obtained after elimination 
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of Tcca' K') between the two equations. It reads 

X°", (115) 



X Xc Xc 



where, for example, x is the matrix formed from 
Xo-.cr' [q, K, K') for fixed q. The charge (ch) and spin (sp) 
susceptibiUties Xch,spiq, T) are deduced from x according 
to 



Xch,spiq,T) 



(kBT? 

m 



E 

KK'aa' 



Ka'Xaa'{q,K,K') 



where Ao-u' ~ 1 for the charge channel and \cci = aa' 
for the spin channel. The calculation of particle-particle 
(i.e. pairing) susceptibilities follows from a straightfor- 
ward generalization of this formalism. The reader is re- 
ferred to prior articles on cluster QMC for more details 
on these topics ((.Tarrell el a,L . .2()0lK] . 



G. Disordered systems 

In this section, we describe cluster approximations for 
disordered systems. The mapping between the lattice 
and the cluster in the CDMFT and the molecular CPA 
(MCPA) is identical, and when the CDMFT is apphed 
to disordered systems, it becomes identical to the MCPA 
l|Duca,ste]1et iTsukadaL ^^M)- The MCPA and re- 

lated techn iques for diso rdered systems are extensively 
reviewed in IConisI lll992|) and will not be discussed here. 

In the remainder of this section we describe a DCA for 
disordered systems that recovers the CPA for Nc ~ 1 and 
becomes exact when Nc oo. This is an extension of 
the fo rmalism introduced bv lJarrell and Krishnamurthvl 
, modified to include the effects of correlated dis- 
order l).Tarrell and John son'. "2004) . We consider an An- 
derson model with diagonal disorder, described by the 
Hamiltonian 

^^-^ (4cj- + c]<,c,,)+^(F,~M)n,, (117) 



where cj^ creates a quasiparticle on site i with spin a and 
nia = c\^Cicr. The disorder occurs in the site occupancies 
labeled by and in the local orbital energies Vi which 
describe the electrostatic potential at site i. We assume 
that Vi are quenched random variables distributed ac- 
cording to some specified probability distribution which 
includes the effects of intersite correlations. 

The effect of the disorder potential ViUi^a- may be 
described using standard diagrammatic perturbation the- 
ory (although we eventually sum to all orders). We per- 
turb around the ordered state, described by the first term 
in Eq. ((TTTji . 

Translational invariance and momentum conservation 
are restored by averaging over all allowed values of the 
site occupancies {^i} and the corresponding disorder po- 
tentials Vi. We study this effect, e.g., on the sixth graph 
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Figure 8 A few low-order diagrams in the irreducible self- 
energy of a quenched diagonally disordered system. Each cir- 
cle represents the scattering of a state k from sites denoted 
by X and a dotted line. The correlations between the elec- 
trostatic potentials on different sites are denoted by the wavy 
line. 



shown in Fig. |H1 which makes a contribution to the self- 
energy 

E (^^.'T^')G(ki)G(k2)G(k3) 

j#J,ki,k2,k3 

^ gixi-(ki-k+k3-k2)gixj-(k2-ki+k'-k3)^-|^-[^g-j 

After averaging over the disorder configurations, (V^^V^) 
becomes a function of Xi — Xj . We identify this average 
as Dfj. With translational invariance restored, we may 
complete the Fourier transform, and obtain 

J2 i?'^(q)G(ki)G(k2)G(k3) 



ki,k2,k3,q 



(A^Vk2+k,ki+k3 - 1) • (119) 



It is easy to extend this argument to all orders in per- 
turbation theory. All graphs are composed of sums of 
products of Laue functions (e.g. A = iV(5q+k2+k,ki+k3), 
and Green functions G(k) and D"™(q), where Z)"™(q) 
is the Fourier transform of DfJ^ = {V^'Vp). 

A hierarchy of approximations may then be con- 
structed by approximating the Laue functions within the 
graphs. These include the CPA where conservation of 
the internal momentum labels is completely neglected. 
Here, all of the Laue functions involving the internal mo- 
mentum labels are set to one, A = iV(5ki+k2- - 1- In 
this case we may freely sum over all internal momentum 
labels, and all terms describing non-local correlations, 
such as those on the bottom of Fig. |H1 vanish, whereas 
the CPA graphs shown on top remain. Different clus- 
ter approximations, including the MCPA and the DCA 
may be constructed by systematically restoring momen- 
tum conservation through the appropriate choice of Laue 
function, as discussed in Sec. III.Bl 

For the DCA, we divide the lattice into clusters as 
described in Fig.|31 write wave- vectors as k = K -I- k and 
employ the DCA Laue function Eq. (|58|) . 



A 



DCA 



NcS- 



Ki+K2,K3+K4 



(120) 



such that momentum conservation is preserved only 
within the cluster reciprocal space K. With this choice 
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for the Laue function, we may freely sum over superlat- 
tice wave- vectors k within each DCA coarse-graining cell 
(see Fig. . This leads to the replacement of the lat- 
tice propagators G(k) and I?"™(k) with coarse-grained 
propagators G(K) and Z)""'(K), respectively, 



N, 



^(K'^) = ^EG(K-t-k,z), 



(K) = §E^""(K + k), 



N 



(121) 



(122) 



The first of these sums is straightforward; however, the 
second requires some investigation due to the powers of 
the potential. 

Here, we calculate Z)"'"(K) for a binary A,B alloy 
where the concentration of A atoms is c and that of B 
atoms is 1 — c. This calculation can be generalized to mul- 
ticomponent alloys and may easily be extended for more 
complex alloys. For the binary alloy, we employ an idem- 
potent formalism where the idempotent f i = 1 indicating 
that the site is occupied by an A atom, or = indicat- 
ing a B atom site. If we associate Vi — Va OJ: Vb for an 
A or B atom, respectively, then Vi — (.iVa + — £,i)VB ■ 
Then, since = 6 and (1 - e»)" = 1-6, 

{V^Vn = {{^^VA + {l-^^)VBn^,VA + (l-^,)VBr) 

= mvx + (1 - ^^)vsmvx' + (1 - CjWS")) 



(123) 



where repeated indices in the last line are summed over 
and 



4" 



(6(1-6)) 

((i-6)(i-6)) 



(124) 
(125) 
(126) 



are the joint probabilities for occupation of sites i and j 
by atoms of the designated types. 

This formalism may be generalized to a multicompo- 
nent alloy a = A, B,C, . . .. Since a site may only be 
occupied by an atom of one type only, in general 



fl7 



(127) 



{V"VJ") is a linear function of the probabilities g^^ for 
all n and m. Thus, the effect of coarse-graining _D"™(k) 
is equivalent to coarse-graining g^^Ck) for all n and m. 

As an example, consider a binary alloy with only near- 
neighbor configurational correlations. Here we may write 



9[ 



g^'"^ + a {25p^ - 1) ,5,;+,,, 



(128) 



where e indexes the near-neighbors to site i, and g^'^^ = 
(?, g^'^^ = (1 - cY and g°'^^ = c(l - c) are the joint 
probabilities for the occupation of different sites for a 



system without configurational correlations. On a hy- 
percubic lattice of dimension D 



/-^(k) = g°'^^(5k,o + a (2<507 - 1) ^ cos{h) . (129) 

1=1 

The corresponding coarse-grained result is 

D 

/^(K) = .g*'^''^<5K,o +«(2<5^7 - l)i?^cos(i^O (130) 

1=1 

where R — (Lc/tt) sin(7r/ic) is a coarse-graining factor 
{Lc = nI^^ is the linear cluster size). If we transform 
back to the cluster coordinates, then 
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■"^-^ ^g^'P^ + Ra{26p^-l)5, 



(131) 



are the configurational probabilities for the cluster, where 
e labels the sites adjacent to i. In the CPA limit, Lc = 1 
and thus R vanishes indicating the lack of any configu- 
rational correlations. In the limit as Nc = — > oo, 
-R = 1, so correlations are systematically restored. 

The cluster problem generated by the substitution 
A — ^ ^DCA may be solved numerically. Each of the 
diagrams in Fig. |S1 representing independent scat- 
ters, or less, remain finite; however, scattering diagrams 
for greater than Nc scatters vanish. The complexity of 
the problem is further reduced since the nontrivial sums 
involve only the cluster momenta K ( numbering Nc 
instead of N). Furthermore, since these diagrams are 
the same as those from a finite-sized periodic cluster of 
Nc sites, we can easily sum this series to all orders by 
numerically solving the corresponding cluster problem. 
The resulting algorithm is identical to that presented in 
Sec. III. except that: 1. In the coarse-graining step we 
must calculate both the coarse-grained correlation func- 
tion 



N, 



/-(K) = i^^/^(K + k). 



and the coarse-grained cluster Green function 



1 



N 



Sc(K,z) 



(132) 



(133) 



2. We then solve the cluster problem by performing a 
weighted average of the cluster Green function (in matrix 
notation in the cluster sites) 



(134) 



over all disorder configurations. The weighting of each 
configuration is determined by the Fourier transform of 
^/37(K) to obtain the cluster configurational probabilities 
g^^ . After convergence is reached, the irreducible cluster 
quantities may be used to calculate the properties of the 
lattice. 
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H. Alternative cluster methods 

In this review we decided to focus on the three, in 
our view most estabhshed cluster methods. This section 
reviews several other ideas proposed to introduce non- 
local correlations, some of which are complementary to 
the approaches discussed. 



1. Self-energy functional theory 

The self-ener gy-functi qnal approach deve loped by 
iPotthofj l|2nn3b(l (see also lPotthoff et allhOO^ is a very 
general unifying concept for existing cluster approaches 
and in addition provides the power to construct novel 
cluster algorithms. Similar to the formalism presented in 
Sec. III.BI this approach views the grand potential as 
the central quantity. Here, the self-energy S is considered 
the basic dynamic variable and a self-energy functional 
17[S] = Tr;n(-(Go + is constructed from 

the Legendre-transform F[S] = $[G[E]] - Tr(SG[S]) of 
the <&[G] functional. This approach proceeds with set- 
ting up a general variational scheme to use dynamical 
information from a (numerically) solvable reference sys- 
tem H' = i?o(t') + ^^i(U) to approximate the physics 
of the original system H = Ho{t) + _ffi(U). While the 
single-particle parts Hq are generally different, the in- 
teraction part Hi{\J) is kept fixed, to ensure that the 
functionals F[E] of the reference system and the origi- 
nal system share the same functional dependence. It can 
then be shown that the grand potential of the original 
system Qt [S] can be evaluated exactly according to 

^t[S(t')] = ^^' + Trln(-(Goi-£(t'))-i) 

- Trln(-G') (135) 

from the grand potential fi', the trial self-energy S(t') 
and Green function G' of the reference system. Varia- 
tion is then performed with respect to the single-particle 
parameters t' of the reference system and the stationary 
point is determined by dV,t['S{t')]/dt' = 0. 

The choice of a suitable and solvable reference system 
H' given the original Hamiltonian H is in principle arbi- 
trary. Within the SFT, cluster approximations are con- 
structed from a reference system H' that represents de- 
coupled clusters of size Nc which are optionally coupled 
to additional ni, uncorrelated bath sites. A schematic 
illustration for the Hubbard model is shown in Fig. O 
A set of decoupled sites as a reference system H' (not 
shown) yields the Hubbard-I approximation. Non-local 
correlations can be included by considering a system of 
decoupled clusters of size A'^c > 1- This approach is iden- 
tical to the CPT when the intra-cluster parameters t' are 
fixed to the original values t : The self-energy of the ref- 
erence system, S(t') = Gg^^ — G'~"^ is calculated once 
and used to approximate the Green function of the orig- 
inal lattice model, G'^ = G^^ - S(t'). This yields the 
expression (071 of the CPT, G'^ = G'"^ - St, since the 
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SFT 
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Figure 9 Taken from iPotthoff et al\ J2003D . Various refer- 
ence systems H' for the 2D Hubbard model on a square lat- 
tice. The solid and dashed lines represent the intra- and inter- 
cluster hopping respectively, the filled circles the on-site in- 
teraction U. Additional rij, bath sites are depicted by open 
circles, while the shaded circles represent Ub = oo bath sites. 



bare Green-functions of the original and reference sys- 
tems only differ by the inter-cluster hopping St. 

Alternatively within the SFT the intra-cluster hopping 
t' can be treate d variationally. Calc ulations for the ID 
Hubbard model l)Potthoff et a/.U2003(l however show that 
the SFT functional Eq. H135|l is stationary for t' slightly 
larger but very close to the original value t, with the dif- 
ference becoming smaller as Nc increases. The SFT can 
also be used to find the correct boundary conditions of 
the CPT cluster. Calculations with fixed intra-cluster 
hopping t' = t and additional variational hopping tr 
between the surface sites to simulate periodic boundary 
conditions show that the self-energy functional is only 
stationary at tr = (open boundary con ditions), not at 
tr = t (periodic boundary conditions) l|Potthoff et all 
120031) . 

One shortcoming of the CPT is its inability to describe 
phases with broken symmetry due to the finite size of the 
clusters. The SFT can be used to construct a variational 
CPT, an extension of t he CP T to account for long-range 
order (iDahnken et all l2003|) . This is done by adding a 
fictitious symmetry breaking single-particle term to the 
reference system and subtracting it in the inter-cluster 
part so that it has no effect on the original lattice Hamil- 
tonian. The magnitude of this term is then optimized in 
the SFT variational treatment. Calculations for the 2D 
Hubbard model show very good agreement with varia- 
tional Monte Carlo and auxiliary field Quantum Monte 
Carlo results for the ground-state energy and the stag- 
gered magnetization ijDahnken et all\200^ . 

Embedded cluster approximations are constructed if 
additional nt uncorrelated bath sites are introduced. The 
DMFT is obtained by choosing as a reference system a set 
of decoupled sites coupled to nf, = oo uncorrelated bath 
sites. The CDMFT is identical to a system of decoupled 
clusters coupled to rif, = oo bath sites with the intra- 
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cluster parameters being fixed at the original values t and 
the bath param eters being treated as variational. It can 
then be shown ijPntthnff e.t nlV l2nn.-^ that the SFT self- 
energy functional is stationary at the bath parameters 
that fulfill the CDMFT self-consistency. A corresponding 
analysis of the DCA formalism within this approach is 
still lacking. 

Intermediate approximations can be constructed by 
considering a finite number < n;, < oo of bath sites. 
This approach is superior to the CPT and inferior to the 
CDMFT. For a reference system with A^c = 1 and ni, = 1 
this approach can be used to s t udy th e Mott transition in 
the Hubbard model l|Potthofj . EinSah analytically and it 
was shown to yield the same qualitative pi cture as the full 
DMFT {Nc = 1, nf, = oo) analysis fsee iGeorges et all . 



2. Fictive impurity models 

Similar to the SFT approach discus s ed abo ve, this very 
general approach by lOkamoto et al\ l)200,'^ is centered 
on the self-energy E as the basic dynamic variable. It 
is based on the idea that the cluster model is merely 
an algorithm to calculate coefficients in an orthogonal 
function expansion of the momentum dependence of the 
electronic self-energy 

(k,z)= $.(k)S,(z). (136) 

i— 0...n 

The coefficients Ei(z) can then be obtained from the so- 
lution of a n-|-l-site fictive impurity model involving n+1 
mean-fields which are fixed by the requirement that the 
impurity model Green functions Gi equal the correspond- 
ing integrals over the lattice Green function 

G,(z) -^$,(k)[G;i(k,z)-Eapprox(k,z)]-^ (137) 

k 

To include local and nearest neighbor correlations, the 
momentum dependence of the self-energy may be ex- 
panded up to second order using the orthogonal functions 
<I'o(k) = 1 and $i(k) = e''^". Since in general the orthog- 
onal functions <&i(k) change sign over the Brillouin zone 
except for the local term i = 0, causality is not guar- 
anteed when the expansion is truncated at low order. 
However it is shown that simple filtering of the higher- 
order terms m ay be used to circumvent these problems 
l)Okamoto et a/... .20031. 

The DGA may be viewed as a specific example of this 
approach, where the indices i correspond to the centers 
K of the DCA cells and the functions <E'K(k) are set to 1 
if k is contained in the cell represented by K and oth- 
erwise (see Fig. 13). Causality problems are thus avoided 
by using orthogonal functions which are non-negative ev- 
erywhere. The resulting approximants however have dis- 
continuities in momentum space. 



3. Non-local effects via spectral density approximation 

iLaad and van den Bossch j l)200Cl() proposed to include 
non-local 1 / Z?-corrections by combining DMFT with the 
Spectral Density Approximation (SDA) (^lioth . 1969) and 
applied this approach to the FKM (see Eq. H182|l 'l for 
which it becomes particularly simple. In the SDA, the 
moments of the spectral function are determined (via re- 
peated evaluation of commutators with the Hamiltonian) 
by complicated but static correlators. For the FKM, 
to order 1/D the SDA self-energy of the d electrons, 
So^(k, z) can be expressed in terms of the static sus- 
ceptibility of the /-electrons. This self-energy is purely 
real, but momentum dependent and is used to approxi- 
mate the bath self-energy in the hybridization function 

^ z - eu + fi- 2^o,d(k, z) 

for the effective impurity problem in = oo. The impu- 
rity self energy. Simp, d (2) (see Brandt and Miclsch, 198^ 
combines with the bath self-energy So,d(k, z) to a dynam- 
ical, non-local self-energy of the form 

T^diK z) = j:i,np,d{^) + ^oAK z)-^ i;o^d(k, z) . (139) 

k 

This finally determines the Green function of the mo- 
bile d- electrons of the usual form, Gd(k, z) = [z — + 
/i — I](i(k, z)]~^, which is used together with I]d(k, z) 
to estimate the susceptib i lity o f the /-electrons (see 
iLaad and van den Bosscha l2000l for details). I]o,d(k, z) 
can then be recalculated to close the self-consistency 
loop. 

With this method, Laad and van den Bossche stud- 
ied the DOS and the spectral function A(k, us) of the 
FKM on a 2D square lattice. Their resu lts agree well 
with known results and other studies like iHettler et al\ 
1)20001) . On the other hand, the FKM is a particularly 
gentle test bed for their method, due to its effective im- 
purity nature (mobile electrons in a static background). 
For a model with true dynamics like the Hubbard model, 
it is unclear whether the method is even feasible (due 
to the additional spin- flip and pair-hopping correlators), 
and the limitation of a purely real bath self-energy is 
likely to be too restrictive. 

4. Non-local corrections via projection technique 

Using the projection technique Tien has developed 
a cluster extension of the DMFT by taking into 
account both local and non-local contributions to 
the dynamics within a relevant subspace of Liouville 
or operator space (Min h-Ticn 1998, 1999a b, 200il 
iTanh-Hai and Minh-Tie'3[ I2OOII) . The information of 
a given subspace is stored in static susceptibility and 
frequency functions while the effects of the remain- 
ing subspace are collected in a dynamic memory func- 
tion. The idea of Tien's approach is to approximate 
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this memory function by a local quantity. Hence, non- 
local correlations are taken into account through static 
quantities (susceptibilities), while dynamical correlations 
are approximated by a local memory function. If the 
relevant subspace is spanned by operators acting on 
the same or nearest-neighbor sites only, these quanti- 
ties can be calculated in an effective impurity model 
(^inh-Tien, 1998, 1 999a|) or , in an imp r oved version , 
in a two-site cluster model llMinh-Tienl. Il999bl l200lt 
iTanh-Hai and. 

This approach has not been rigorously proven to 
be caus al, however its applicat ion to the Falicov- 
Kimball ("Minh-Tienl, Il998l Il999b|) and Hubbard model 
(Minh-Ticn, 2001) shows that the spectral function is 
positive definite and the sum-rules of the first few mo- 
ments of the spectral densities are preserved. The ne- 
glect of non-local dynamical correlations however leads 
to spurious behavior at low temperatures: A Kondo res- 
onance emerges in the in 2D half-filled Hubbard model 
at low temperatures similar to the behavior observed in 
DMFT simulations but inconsistent with other cluster 
calculations which show a pseudogap down to the lowest 
temperatures (see Sec. IIV.D.2p . This shortcoming is due 
to the fact that the nearest-neighbor static correlations 
are proportional to the band-dispersion ek. At the Fermi 
energy, ek = 0, hence non-local correlations vanish and a 
Kondo-resonance is generated due to the local nature of 
the dynamical correlations. 



includes only on-site and nearest-neighbor correlations, 
it already captures important signatures of correlations 
consistent with DCA/QMC results for larger clusters 
such as the existence of a Mott-Hubbard gap in the 
2D Hubbard model fo r all v alues oi U > at low T 
l|Stanescu and Phillipgl. l200l[l and the emergence of a 
pseudogap in the density of states due to antiferromag- 
netic correlations (,Stanescu and PhillipsK .2003ai see also 
Sec iryTPll . 



III. QUANTUM CLUSTER SOLVERS 

Cluster techniques map the lattice system onto a self- 
consistently embedded quantum cluster model. This 
chapter discusses the most promising numerical ap- 
proaches used to solve this cluster problem. After stress- 
ing some general difficulties faced by potential cluster 
solvers in Sec. IHI.AI we present several perturbative 
techniques, including second order perturbation theory 
in Sec. IIII.B.ll the fluctuation-exchange approximation 
in Sec. IIII.B.2I and the non-crossing approximation in 
Sec. IHI.B.3I Section IIH.CI reviews the application of 
non-perturbative, (numerically) exact techniques includ- 
ing Quantum Monte Carlo in Sec. IIII.C.ll Exact Diago- 
nalization (ED) in Sec. IIII.C.2I and the numerical renor- 
malization group in Sec. IIII.C.3I 



5. Two-site correlations with composite operators 

A dynamical non-perturbative two-site approximation 
for the Hubbard model based on the composite opera- 
tor m ethod was developed by iMatsumoto and Mancinil 
l)l997j) a nd later adopted and improve d by Stanescu and 
Phillips l|Stanescu and PhillipsL 12001 . 2003a). By using 
Hubbard operators as a local basis which exactly di- 
agonalize the interaction part of the Hamiltonian, this 
approach recovers both the weak coupling U <^ t and 
strong coupling U ^ t limits of the Hubbard model. 
The memory function (5m(k,w) which collects the effects 
of dynamical correlations is expanded in a two-site ap- 
proximation as ^m(k, oj) = 5mo(w) -I- ak(5TOi(w) where 
ak = (cos fc^ + cosfcj,)/2 and Sttiq (Smi) are the lo- 
cal (nearest-neighbor) contributions. Hence this method 
contains on-site and nearest-neighbor dynamical correla- 
tions. Unknown quantities are expressed in terms of re- 
solvents for the eigenstates of a two-site impurity system. 
The renormalization of these resolvents due to their cou- 
pling to the surrounding of the two-site system is treated 
within the non-crossing approximation. 

The application of this technique to the Hubbard 
model s hows qualitative agreemen t of th e single-particle 
spectra jMatsumoto and Mancini', '1997^ in 2D with fi- 
nite size QMC results and high accuracy of specific heat 
results ((stanescu and Phillipstl2001i) as compared to the 
Bethe ansatz solution in ID. Although this technique 



A. General remarks 

The fundamental difference between a finite size clus- 
ter and the effective cluster problem of quantum cluster 
theories is the existence of additional quantum mechani- 
cal bath degrees of freedoms in the latter. The simplest 
realization of such a syste m is of course th e well-known 
Anderson impurity model l|Andersoii Il96l[l . In general, 
its ground state is a nontrivial many-body state, which is 
not perturbatively connected to states of simpler starting 
points like a non-interacting or a free impurity. Further- 
more, the excitations also are of many-body character 
and typically involve dynamically generated low-energy 
scales which depend non-analytically on system param- 
eters. Consequently, any perturbation theory is faced 
with severe limitations concerning its region of applica- 
bility and the most successful techniques used to solve 
this fundame ntal pr oblem of solid state theory are non- 
perturbative l|Hewso n. 1993). 

Nevertheless, a variety of tools to approximately or nu- 
merically sol ve this model h ave been developed over the 
last 25 years (|HewsonLll993l) . Since the physics o f the An- 
derson impurity model is very well understood ijHewsonl 
1993), this knowledge can be employed to judge the qual- 
ity of results and region of applicability of these various 
analytical or computational techniques at hand, a priori 
as well as a posteriori. This statement does also apply to 
a large extent to the DMFT, where an e ffective Ander- 
son impurity model plays the central role l|Georges et al\ . 
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ll99fiHPruschke pJ, a/llTool . 

The situation becomes much more involved for quan- 
tum cluster problems. First, from a purely technical 
point of view, the complexity of the system can limit the 
applicability of a method on principal grounds or even 
rule it out as a potential cluster solver altogether. This 
aspect will be discussed for every technique in detail in 
the following sections. 

Similarly important is a novel complexity of the physics 
entering already on the level of the bare cluster and 
largely complicating the task to judge the reliability and 
quality of results obtained with a certain method. With 
increasing cluster size not only high-energy states (charg- 
ing energies) - which are the only one present in the im- 
purity case - will couple to the quantum bath but also 
medium- and to some extent low-energy states (e.g. non- 
local spin fluctuations). Moreover, the latter already can 
and in general will contain interesting many-body effects. 
The coupling to the bath can then lead to a subtle inter- 
play or even competition between such intrinsic many- 
body effects and those introduced by quantum fluctu- 
ations in the external bath. A well-established exam- 
ple is the occurrence of a qua ntum phase-transition in 
the two-impurity Kondo model l|.Tones and Varm 
'Jone s et a/.Ul988|) driven by the competition between the 
RKKY magnetic interaction favoring a non-local singlet 
and the Kondo effect with its local singlet formation. 

In contrast to the single impurity case where the qual- 
itative results to be expected from a calculation can in 
principle be read off the input parameters, such a vi- 
tal possibility of a plausibility check for a given method 
does not exist in case of quantum cluster models. Ex- 
perience tells that even at first sight physically plau- 
sible results need to be checked carefully. Again, the 
two-impurity Kondo model can serve as a pedagogi- 
cal example. Th e original, quite intuitive result by 
iJones etai] l)l988|) has later been reexamined by others 
and fou nd to be valid only under very special circ um- 
stances fSaka i and ShimizuL Il99l ISakai erall . 11990*). A 
fully satisfactory understanding of the physics of this sim- 
plest quantum cluster model is in fact lacking until today. 
Things are further complicated for the case of cluster 
mean-field theories due to the "backflow" of the complex 
local physics via the effective bath. An a priori under- 
standing of the behavior of the system is, at least given 
the current level of knowledge, virtually impossible, but 
also an a posteriori plausibility check is rather based on 
subjective physical intuition than on solid understanding 
of the basic physics. 

Thus, in order to obtain a reliable and consistent un- 
derstanding of the physics of correlated electron systems 
in the framework of cluster mean-field theories, it is vital 
to employ a variety of complementary and possibly non- 
perturbative tools to solve the effective quantum cluster 
model. We therefore discuss in detail different computa- 
tional and analytical tools regarding their applicability 
and reliability. 



B. Perturbative techniques 

The numerical effort to compute the cluster self-energy 
increases rapidly with the cluster size, in principle ex- 
ponentially for exact methods. This calls for simpler 
methods for which the complexity of the problem can 
be reduced. This can be achieved by perturbation the- 
ory in its many variations. One example is the standard 
weak-coupling perturbation theory in the electronic in- 
teraction U. Second Order Perturbation Theory (SOFT) 
and the more elaborate Fluctuation Exchange Approxi- 
mation (FLEX) are discussed in the following sections. A 
complementary strong-coupling approach handles the in- 
teraction exactly, but treats the coupling to the host F in 
a perturbative expansion. An example for this approach 
is the Non-Crossing Approximation (NCA), discussed in 
Sec. liii.B.31 



1. Second order perturbation theory 

Second order perturbation theory has proven very use- 
ful to solve the effective impurity problem of the DMFT. 
This iterated perturbation theory (IPT) becomes ex- 
act in the weak coupling limit J7/t <C 1 and coinci- 
dently at half-filling even in th e strong coupling limit 
U/t > 1 l)Georges all Il996|) . The IPT was shown 
to give results in good qualitative agreement with non- 
perturbative DMFT results, e.g. it captures the Mott- 
transition at half-filling. 

The application of second order perturbation theory in 
the context of quantum cluster theories follows immedi- 
ately from the formulation in Sec. lII.ijl bv considering the 
first and second order diagrams in the cluster functional 
^[G] and neglecting higher-order diagrams. Here we il- 
lustrate the formalism for the Hubbard model with local 
interaction U. For the DCA, an expression for the cluster 
self-energy Sc = 6^[G]/SG may be obtained by introduc- 
ing the bare particle-hole susceptibility, i.e. the " bubble" . 
With the notation K = (K, icom), Q — (Q, ii^n) we have 

UT 

Xvh{Q) = -^Y.^{K + Q)G{K), (140) 

^ K 

where T is the temperature and G is the coarse-grained 
Green function defined in Eq. H59|) . Note that ujm is 
a fermionic Matsubara frequency, while Vn is of bosonic 
type. The cluster self-energy to second order in U is then 
given by 

UT 

Se(i^) - 1^ E Xph{Q)G(K - Q) . (141) 
Q 

We ignored the constant (first order) Hartree term which 
just shifts the chemical potential. Physically, the self- 
energy accounts for the interaction of electrons with 
virtually excited electron-hole fluctuations described by 
Xph- Eq. H141|l is a convolution of the susceptibility and 
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the Green function in both momentum and Matsubara 
frequencies, and thus can be readily computed via Fast 
Fourier transform (FFT). FFTs scale only logarithmi- 
cally with size, so a much larger range of cluster sizes is 
available for perturbative studies of this kind. For real 
frequencies, the use of FFT is less appealing, since FFT 
also requires an equidistant grid of frequencies. Com- 
pared to the usual expression for a finite size system, in 
the DCA we use the coarse-grained Green function G for 
evaluating the susceptibility and the self-energy, so the 
momentum sums run over the Nc cluster momenta K, Q. 
For the CDMFT, one can formulate similar expressions, 
but, as discussed in Sec. ^ the self-energy and Green 
function are matrices in cluster real space. The matrix 
inversions needed during the iteration become costly as 
the cluster size is increased. 

It must be emphasized that simple minded application 
of FFT to Green functions (in particular with fermionic 
Matsubara frequencies) is likely to incur large errors, due 
the discontinuity at zero imaginary time that leads to the 
slow 1/iuJn decay of the Green functions and self-energies. 
The way to deal with this is to manually subtract the 
slowly decaying parts before application of the FFT. This 
is possible as the high frequency behavior can be analyti- 
cally determined. The subtraction leads to functions that 
decay sufficiently fast to have a numerically well defined 
transform, even though one has to cut off the Matsubara 
sums in Eqs. (|140I141|I at a finite va lue. Details on thi s 
subtraction technique can be found in lDeisz et a/.M2003(l . 



2. Fluctuation exchange approximation 



FLEX ^Bickers et all 11989: Bick ers and Whitel 1199(1^) 

extends the second order perturbation theory through 
higher order diagrams. It rests on the assumption that 
the interacting electron system can be considered as a 
problem of electrons exchanging self-consistently deter- 
mined fluctuations of various kinds, i.e. density, spin 
and pair fluctuations. In the context of quantum clus- 
ter tlieori^Sj_tlie^]LEXw^ applied within the DCA 
bv lArvanpour et al\ ll2003al) . Here we briefly review the 
corresponding formalism. The generating functional of 
the FLEX is a sum o f the three fluctuation contributions 
l|Bickers et a?J . ll989|) 



where 
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and the trace Tr denotes (T/Nc) J2q J2n- addition to 
the particle- hole bubble Eq. (|14()|l . the particle-particle 
susceptibility Xpp also appears and is defined by 



Xpp( 



(146) 
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With the generating functional of Eq. H142I) the cluster 
self-energy can be written as 

^c{K) = —Y,[v('^'^\Q)G{K-Q) 



-V''PP){Q)G{Q-K)], (147) 



where the so called FLEX potentials are given by 
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The FLEX potentials constitute an infinite sum of fluctu- 
ation diagrams similar to the series of density fluctuations 
also known as the Random Phase Approximation, which 
is in fact a subset of the FLEX. However, the FLEX 
also includes significant spin and pair fluctuations. Note 
that second order perturbation theory is reproduced by 
expanding Eq. (|148|l to first order in Xph- 

Naturally, the drawback of using these weak-coupling 
methods lies in the fact that some strong coupling phe- 
nomena such as the opening of a Mott gap at half-filling 
simply cannot be described by e.g. the FLEX approach. 
Also, quite often the self-consistent evaluation simply 
does not converge, particularly at low temperatures. This 
is due to divergences showing up in perturbatively eval- 
uated response functions. As a result, the variety of sys- 
tems that can be studied is limited. 

To obtain real time (or frequency) data, an analytic 
continuation of the Green function or self-energy is neces- 
sary. For the FLEX, this is achieved by m eans of Fade ap- 
proximation ijVidberg and Sereii3 . ll977(l . Real frequency 
results and details on the im plementation of the Fade a p- 
proximation can be found in lArvanpour et al\ l)2003bf) . 



(142) 3. Non-crossing approximation 



The Non-Crossing App r oxima ti on (NCA) (iGrewe , 
198^; 'Keiter and Czvcholl, Il983t iKeiter and Kimbal , 



1971; Kuramoto. J.983) was originally developed for the 
single-impurity Anderson model a nd is based on a di- 
agrammatic perturbation theory ijKeiter and Kimball 
1970; Kuramoto, 1983) around the atomic limit of this 
model. A comprehensive and detailed description of 
the NC A and it s limitations can be found in the re- 
view by iBickersI (|l987). The NCA has been exten- 
sively applied as an approximate solution of the effec- 
tive impurity model in the DMFT context to study the 
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NCA has been extend ed to solve t he effe ctive cluster 
problem of the DC A bv lMaier et'oB l|2000bj) . The aim of 
this section is to convey the general concepts of the NCA 
method and discuss the extensions necessary to solve the 
effective cluster model of quantum cluster theories. 

The perturbation theory used to construct the NCA 
is based on a resolvent technique and is in general 
applicable to problems where a discrete level system 
not representable by standard fermionic or bosonic op- 
erators is coupled to continuous degrees of freedom 
l|Becker and Kelleilll987: .Kuramota. l98a) . The discrete 
level system in the context of the DCA is given by the 
eigenstates of the local part of the cluster Hamiltonian, 
Eq. 



Hc,loc = Xl'^'^K -A*)ck^' 



Ka 



(150) 
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The perturbative expansion is performed with respect to 
the coupling 

Hhyb = [^K(k)cL«K+k. + ^■'^■l (151) 

to the continuous degrees of freedom, i.e. the auxil- 
iary non- interacting fermions a^.^. (al^^). Hence the NCA 
should be an especially useful approximation where the 
local energy scales (in this case U) exceeds the magnitude 
of the coupling V to the host. 

To apply the resolvent technique, the fermionic clus- 
ter operators c^^, (c^^.) are expanded in terms of the 
Hubbard-operators X„in — \m){n\ 



(152) 



where are the eigenstates of the local Hamiltonian, 

Eq. linnil . 



Hc,loc ~ ^ ^ -^m-^n 



(153) 



with eigenenergies Em and — (r7i|cKcr|«). In this 

representation, the hybridization term Eq. (|151|l becomes 

Hhyb= E E[^K(k)Ci?;'*X„™aK+k. + /^-c.]. (154) 

Since the Hubbard operators X^n do not obey stan- 
dard fermionic commutation relations, the conventional 



Feynman diagram technique cannot be used for a pertur- 
bational expansion and the concept of r esolvents must 
be introduced instead ijKuramotoL The matrix- 

elements of these resolvents in the space of the local 
eigenstates have the form 



(155) 



where the resolvent self-energy I]^(z) collects renormal- 
ization effects of the individual molecular states 
due to the hybridization Eq. H154(l to the host. We note 
that E-'^(z) and hence R{z) reduce to diagonal matrices 
if the hybridization term Hhyb commutes with the lo- 
cal part Hc.ioc of the cluster Hamiltonian. While this 
is the case in the paramagnetic state we stress that the 
non-diagonal el ements become essential in the sym metry- 
broken phases (|Maieil l200lt iMaier et a.l\ . l2000a|) . 




Figure 10 Non-crossing diagrams for the resolvent self-energy 
E^. The dashed line stands for the dressed resolvent R de- 
fined in Eq. 115511 . the full line for the host electrons. 

Keeping only the lowest order, self-consistent diagrams 
for the resolvent self-energy illustrated in Fig.llOlde- 
fines the NCA. Since the hybridization with the host de- 
grees of freedom enters the expressions only in the form 
of a density of states 
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there is no need to explicitly calculate the auxiliary host 
parameters l^(k) and Ak of the cluster Hamiltonian, 
Eq. H83|) . To ensure that the solution of the cluster prob- 
lem is also the solution of the coarse-grained lattice prob- 
lem, Fc(K, z) is replaced with the hybridization function 
r(K,z) (see Eq. (|82ll ) determined from coarse-graining 
the lattice Green function. After observing energy, mo- 
mentum and spin conservation, the analytical expression 
for the resolvent self-energy within the NCA reads 

K,(T 

def{e)^mr{K,e)C^!,Ru'{z + e)C^l* 

JO 

+ J def{-e)Qmr{K,e)C^l*Ru'{z-e)C^", 

— oo 

The coupled equations (| 15 5(1 and ((157|l are solved self- 
consistently to compute the matrix-elements of the resol- 
vent R(z). By using the expansion Eq. H152|) . the cluster 
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Green function Gc(K, z) = {{cj^cr^ '^K.a)) z is written in 
terms of Hubbard operators 

(158) 

Within the NCA, the correlation function on the right 
hand side reduces to a convolution between resolvents 

+00 

{{Xmm']Xn'n))z ^ dee"^^ [p„m (e)i?m'„' (e + 

— 00 

-Pm'n' (e)Km " ^)] : (159) 

where pnm{ijj) = —l/ir'^mRnmi^ + *<5) is the spectral 
density of the resolvents, Zc = X^n /_t^ '^se~''^p„„(£) 
the cluster partition function and /? = 1/T the inverse 
temperature. 

Although the application of the NCA was illustrated 
for the DCA cluster model, the NCA can in principle 
also be applied to solve the CDMFT cluster problem, 
Eq. 173). The fact that in the CDMFT the hybridization 
function F is a matrix in the Nc cluster sites however 
complicates this task considerably. 

One of the great advantages of this resolvent pertur- 
bation theory and the NCA as approximation is that it 
can be solved for real frequencies, i.e. a cumbersome and 
problematic analytic continuation as in the case of QMC 
simulations is not necessary. Moreover, one can show 
analytically that it can produce a low-energy scale de- 

t mding non-analyt ically on the bare model parameters 
/liiller-HartmannL Il984j) and that certain general fea- 
tures like universality in this energy scale are present, 
too fFischcr, 1997). Surprisingly, the NCA is capable 
of recovering a variety of complex features in thermo- 
dynamic and dynamic quan tities of the sin gle impurity 
Anderson model at U — oo (jBickersl Il987j) . However it 
was realized rather quickly from th e numerical solution 
of the NCA equations (^Grewe. 1983: iKuramotol 11983^1 as 
well as an approximate a nalytical treatment at T = 
l|Muller-Hartmannl Il984|) that the NCA shows serious 
deficiencies in its low temperature/low energy behav- 
ior. More precisely, it generates unphysical power-law 
struct ures in the low-energy be havior of physical quan- 
tities l)Miiller-HartmannL Il984|) . Nevertheless there ex- 
ists a finite window from high down to temperatures 
of the order of an effective Fermi liquid scale where 
the NCA pro duces reliable and physically sensib le re- 
sults f Bickers, '198/; "Co x and Gre!^ll987t iFischeil Il997t 
iPruschke et ai . IMSb) . 

A rather different question is whether these results are 
only qualitatively or even quantitatively correct. Here, 
the answer critically depends on the system under con- 
sideration. Going beyond a single Anderson impurity 
at U = oo introduces the necess ity to include dia- 
grams beyond the lowest ord er NCA llBickers et a^J . ll98^ 
IPruschke and Crewel IT989>1 to at least ensure a correct 
reproduction of energy scales or exactly known limits 



(like for instance uncoupled sites in a cluster calcula- 
tion). An inclusion of these so-called vertex corrections 
to the bare NCA defined by Eq. I|157ll has up to now been 
possible only for the single impurity model with finite U 
l|Pruschke and Crewe.. 1989.) , because the complexity and 
number of the a dditional diagrams r ather fast exceeds all 
practical limits ijHeindl et a/.l . [2000(l . 

Neglecting vertex corrections in practice means that 
the magnitude of a possible Fermi liquid scale will be 
grossly underestimated. However, in the case of strong 
non-local fluctuations, as one would expect in 2D prob- 
lems, such a scale will in general be extremely small or 
even zero anyway. Here, the influence of vertex correc- 
tions can be expected to be less severe and the NCA even 
quantitatively reliable. 

Even with these restrictions, the NCA still becomes 
formidable as the cluster size Nc increases. The number 
of eigenstatcs and hence number of coupled equations 
(|155l I157|l that have to be solved self-consistently grows 
like Mj^p, where riimp is the number of Fock states of 
the isolated impurity (e.g. Uimp = 4 in the Hubbard 
model, Uimp = 3 in the t-J model). Thus, although the 
study of larger clusters is in principle possible, the NCA 
technique is very limited in cluster size and so far has 
only been applied to A^c = 4 size clusters. 



C. Non-perturbative techniques 

Non-perturbative techniques solve the effective cluster 
problem (numerically) exact. This advantage comes at 
the expense that only small cluster sizes can be treated. 
The size restriction on the cluster QMC technique de- 
scribed in the following section is basically the same as in 
the case of ESS QMC. For the ED and NRG approaches 
discussed in the next sections, however, a Hamiltonian 
for both the cluster and host degrees of freedom needs to 
be explicitly simulated, further limiting the usefulness of 
these methods. 



1. Quantum Monte Carlo 

a. Introduction. Quantum Monte Carlo (QMC) is a 
powerful and general technique for quantum cluster prob- 
lems. QMC has several advantages including the ability 
to treat relatively large clusters, the simplicity of the re- 
quired code, and the fact that only the cluster excluded 
Green function Q and the coarse-grained interaction U 
are required as inputs. QMC is also numerically exact 
with small and controllable sources of systematic and 
statistical error. Its disadvantages include the fact that 
there is a minus sign problem which is unpredictable, 
difhculties in calculating real- frequency results, and the 
numerical expense of the approach. 

The QMC algorithm for clusters is based on the Hirsch 
and Eye algorithm which was developed to simulate the 
Anderson impurity problem ijHirsch and Fvel Il986|) . It 
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was later generalized to s olve the DMF T an d DC A efFec- 
tive imp urity problems bv lJarrelll l)l992() and lJarrell et al\ 
l|2n01tJ) . respectively. Although the algorithm is formu- 
lated using a path integral in imaginary time, the Maxi- 
mum Entropy Method may be used to analytically con- 
tinue the QMC data to obta in real-frequency spectra 
l)Jarrell and GubernatilllQQefl . 

The QMC algorithm uses the action of the effective 
cluster model, Eq. lf77|) as a starting point and there- 
fore can be equally applied within the CDMFT and the 
DCA. It thus requires as inputs the initial bare Green 
function C/, and the form of the coarse-grained interac- 
tion"^ U. This is an advantage, since due to the required 
self-consistency of embedded cluster techniques, we gen- 
erally do not know the auxiliary parameters of the cluster 
Hamiltonian. 

Several steps are required to evaluate the path integrals 
of this action using QMC. We first introduce Hubbard- 
Stratonovich (HS) fields required to decouple the inter- 
action, transforming a problem of interacting electrons 
and bosons to one of non-interacting particles coupled 
to time-dependent HS fields. The fermionic and bosonic 
fields are then integrated out, and the integrals over the 
decoupling fields are performed with a Monte Carlo algo- 
rithm. Measurements of any diagrammatic quantity are 
accomplished by decomposing the associated operators 
using Wick's theorem (both connected and disconnected 
contractions must be included) and then averaging the 
result over the Monte Carlo generated field configura- 
tions. 



b. QMC for the simple Hubbard Model. As a specific ex- 
ample, we consider the simple Hubbard model with a 
local interaction of strength U. Since this interaction 
is local, it is unaffected by coarse-graining, so U — U. 
To approximate the time-integrals in the action, we in- 
troduce a discrete time grid of length L and time step 
At = (3/L. The interacting part of the action is then 
decoupled by mapping it to an auxiliary Ising field via 
a discre te Hirsch-Hub bard-Stratonovich (HHS) transfor- 
mation l)HirschL Il983|) , 

g-Ar;7i:,K-r-i/2)(".i-i/2) 

^ ig-ArC//4-Q ^ gas,KT-"u)^ (160) 

^ i Si=±l 

where cosh(Q;) — ^'^'^^ 1"^ and the index i denotes the 
cluster sites. With this change, the cluster action takes 
the form 



^ It is also possible to treat t he interaction in a self-consistent 
manner, as was described in iHettler ef "oil 120001) 



" XlXUI"'^M,'Tf^SM7M-i,<T (162) 

i=l l = \ a 

where I denotes the time slice t/ and Q{i,l;i',l') = 
Gii' [ti — Ti' ) the cluster excluded Green function defined 
in Eq. H7fci|) . Now we integrate out the remaining clus- 
ter Grassmann variables. The partition function then 
becomes 

= Tr{s^^,}l[det{C,,^y^ (163) 

a 

where factors which are fixed during the QMC process 
have been ignored. {Gca-^su) inverse cluster 

Green function matrix with elements 

{Gca;sa)lli,i> = Gi^liM ~ S^-jSi, j^iaas^^i . (164) 

A Monte Carlo algorithm is used to perform the re- 
maining integral over the HHS fields. The Markov 
process in this algorithm proceeds by suggesting local 
changes of the HHS fields at one point in space-time. 
These changes are accepted or rejected according to the 
change in their Boltzmann weight, the argument of the 
trace in Eq. H163|l . If the change is accepted, the Green 
function must be updated accordingly. Several approxi- 
mations, i.e. changes to these equations that are beyond 
linear order in At, are necessary to obtain an efficient 
algorithm. First, we re-exponentiate the first term on 
the right-hand side of Eq. 1164|) . obtaining in a simple 
matrix notation in space-time 

G,-' ^g-^ - Tie""' -1) , (165) 

where T is SijSi-i^i and Va-{i,l) = asija. Note that the 
term in the parenthesis is beyond zeroth order in At. 
Therefore, to first order in At, we may write 

Gc-' = + {g-' - 1) (e^- - 1) , (166) 

since g^^ — 1 = — T + 0{At). Therefore, the inverse 
Green functions for two different field configurations, 
{sii} and {sjj}, are related by 

G.r'e-^^ = G^-'e-""- - e"^" + e"^^ . (167) 

Or, after multiplying by e^^ , and collecting terms 

G,r' - Gc-' = {G,-' ~ l)e-^"(e^^ - e^") . (168) 

Multiplying from the left by Gc and from the right by 
Gc , we find 

G,; = G,„ + (G,, - l)(e^<^-^' - l)Ge; (169) 

or 

Gc.G^r' = 1 + (1 - G,,)(e^^~^" - 1) . (170) 
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The QMC algorithm proposes changes in the Hubbard- 
Stratonovich field configuration {si^} {s'n}, and ac- 
cepts these changes with the transition probability Ps~^s' ■ 
Thus, to implement the algorithm, we need Ps-ts' and a 
relation between the cluster Green functions Gc and Gc 
for the two different auxiliary field configurations. To 
simplify the notation, we introduce a combined space- 
time index i = and consider only local changes 
in the fields Sm — > Kn = "Sm- As can be inferred 
from Eq. pt)3|l . the probability of a configuration {si} 
is Ps (X dct{Gc^g.y) det{GcJ^^ y); ^ on the other hand 
detailed balance requires Ps'Ps'^s — PsPs^s' for all s' . 
We may satisfy this requirement either by defining the 
transition probability Ps'^s = ^/(l + R)i where 



n _ det(G'4)det(G/J 
Ps' " det(Gct)det(Gex) 



(171) 



is the relative weight of two configurations, or by setting 
Ps'^s — minimum(_R, 1). The first choice is called the 
"heat bath" algorithm, and the second the "Metropolis" 
algorithm which is generally, but not always, superior. If 
the difference between two configuration is due to a flip of 
a single Hubbard Stratonovich field at the mth location in 
the cluster space-time (^Hirsch and Fva.l98fi) . we obtain 
from Eq. ifTTHjl 



l)]-i.(172) 



For either the Metropolis or the heat bath algorithm, if 
the change is accepted, we must update the Green func- 
tion accordingly. The relationship between G^ and Gj 
is given by Eq. Hl()9|l 



G, 



1 + (1 - G,.,„™)(e— (^"-^;J - 1)^-^'^'*''^ 

The QMC procedure is initialized by setting Gca,ij ~ 
Qij, and choosing the corresponding field configuration 
with all Si = 0. Then we use Eq. H173|l to create a 
Green function corresponding to a meaningful field con- 
figuration (i.e. Si = ±1, for each i = or the {si} 
from a previous run or iteration) and proceed by step- 
ping through the space-time of the cluster, proposing lo- 
cal changes Si — > —Si. We accept the change if Pg'^s is 
greater than a random number between zero and one and 
update the Green function according to Eq. H173|l . After 
roughly one hundred warm-up sweeps through the space- 
time lattice of the cluster, the system generally comes 
into equilibrium and we begin to make measurements. 



* If Ps is not positive definite, then \Ps\ is used as the sampling 
weight, and its sign is appended to the measurement. I.e. for a 
measurement m and sign S, (in) p = {mS)^p^ / (S)^p^. 



c. Measurements. Several points must be considered 
when making measurements in the QMC procedure. 
First, for a given configuration of the HHS fields, the 
problem is non-interacting. Thus, the estimators of any 
correlation function may be constructed by taking all al- 
lowed Wick's contractions (both connected and discon- 
nected). Therefore, any quantity which may be repre- 
sented in terms of the Green functions (and perhaps the 
HHS fields themselves), may be measured. Second, great 
care must be taken when constructing the estimators of 
the measurements. Different estimators may yield differ- 
ent results due the systematic and statistical error in the 
QMC procedure. It is important to choose the optimal 
form of the estimator of each measurement, and then use 
all the prior knowledge (exact limits, symmetries, etc.) 
that we have to reduce the error. 

For example, one difficulty encountered with the QMC 
algorithm is that a reliable transform from imaginary- 
time quantities to Matsubara frequencies is required. A 
careful treatment of the frequency summation or the 
imaginary-time integration is crucial in order to ensure 
the accuracy and the stability of the algorithm and 
to maintain the correct high-frequency behavior of the 
Green function. We need to evaluate the following inte- 
gral 



Gc(K,zc^„)= / dTe*"""Ge(K,r) 
Jo 



(174) 



But from the QMC, we know the function Gc(K, r) only 
at a discrete subset of the interval [0, (3] . As may be read- 
ily seen by discretizing the above equation, the estimation 
of Gc(K,ia;„) becomes inaccurate at high-frequencies. 
This is formalized by Nyquist's theorem which tells us 
that above the frequency Uc — unpredictable results 
are produced by conventional quadrature techniques. For 
example, a rectangular approximation to the integral in 
Eq. H174|) yields a Gc(K, iw„) that is periodic in a;„. This 
presents a difficulty since causality requires that 



lim Gc(K, iujn) = 



(175) 



This problem may be avoided by using the high fre- 
quency information from other sources. For example, the 
second-order perturbation theory result, Eq. (|141|l has 
the correct asymptotic behavior. Alternatively, we can 
use the Green function from the previous iteration which 
has the correct high frequency behavior. We then com- 
pute the Matsubara-frequency Green function fro m the 
imag inary-time QMC Green function as (,Jarrell et all . 
119931) 

Gc(K,iun) = G,,pt(K,ic^„)+ (176) 
dTe'""^(G,(K,r)-Ge,pt(K,T)) . 



/ 



where Gc,pt is any Green function with the correct high- 
frequency behavior. The integral is computed by first 
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interpolating the difference G'c(K, r) — GcptCKiT) us- 
ing an Akima spline, and then integrating the spline 
(a technique often called oversampling) . The smooth 
Akima spline suppresses spurious high frequency behav- 
ior so that G'c(K,iw„) = Gc,pt (K, 107,1 ) when ijj„ ^ ujc- 
The resulting self-energy, extracted from Sc(K, zcj„) = 
l/QCKjiujn) — l/Gc(K,iw„) is still not accurate at high 
frequencies since it reflects either the perturbation the- 
ory or the previous iteration. However, the exact high 
frequency form is known. Thus, we can fit the high fre- 
quency self energy, at frequencies at or below the Nyquist 
cutoff 7r/(AT), to the form 

E,(K, iuj,,) « a(K)/(ic^„) + b{K)/iul . (177) 

In some cases a and b arc known exactly. We may 
then append this form onto the self-energy at frequen- 
cies higher than 7r/(Ar). 

2. Exact diagonalization 

The Lanczos ED method llHavdock et all \197^ pro- 
vides a numerical way to diagonalize the Hamiltonian of 
a finite size system. Since the cluster model in the CPT 
formalism is identical to the model of a finite size system 
with open boundary conditions, the regular Lanczos ED 
method can be applied in this case without modification 
ijSenechal et al.. . 2QQQ) ■ For the cluster model of embed- 
ded cluster theories however, it needs to be generalized 
to account for the self-consistent coupling to the host. 
ED has been used to solve the impurity problem o f the 
DMFT llCaffarel and KrauthL IToM ISi et all Il994l) . It 
is natural that it can be extended to be used as a clus- 
ter solver i n quantum cluster t heories. This was demon- 
strated bv lBolech et all l)2003(l in an application of the 
CDMFT to the ID Hubb ard model. The gener al method 
of exact diagonalization ijHavdock et all 1975|) need not 
be reviewed here. We only outline the specific implemen- 
tation necessary for a cluster theory. 

ED is a wave function based method, i.e. it is applied 
to diagonalize the effective cluster Hamiltonian, e.g. the 
Hamiltonian in Eq. (|79|l . To this end, the hybridization 
function T{z) (see Eq. (|78f) ') is fitted to the form 

r^,(^) = E - A(k)]rJiC,(k) , (178) 

!m,k 

to obtain estimates for the auxiliary parameters of the 
host, \ij(k) and VijCk). In order to apply ED to the 
cluster Hamiltonian, the auxiliary host degrees of free- 
dom, i.e. the sum over k, must be discretized to a finite 
set of Nh orbitals. Applying ED to the resulting Hamil- 
tonian, one can then compute a cluster self-energy S'^ 
llCaffarel and KrauthL [199^ . and from the cluster self- 
energy a new estimate of the coarse grained Green func- 
tion G to close the self-consistency loop. 

The Hilbert space of the resulting Hamiltonian in- 
creases exponentially with the cluster size Nc and the 



number of wave- vectors Nh representing the host. For 
ED to be feasible, the total number of "sites" N = 
Nc + Nh must be of the order of '--^ 20. Furthermore, 
for the distinction of the lattice into cluster and host to 
make sense, Nh > Nc. This size limitation is possibly 
prohibitive for anything but ID systems. Moreover, by 
applying ED to solve the cluster model one abandons 
the thermodynamic limit, i.e. one of the advantages of 
cluster mean-field theories. Still, ED is an exact method 
and can deal with more complicated interactions than 
the simple on-site Hubbard repulsion, as exemplified in 
iBolech et aO tOO^ . where nearest neighbor interactions 
were considered. 



3. Wilson's numerical renormalization group 

Over the past ten ye ars, Wilson's numeri cal renor- 
malization gr oup (NRG) l)Krishnamurthv et aL 1980a„bl: 
IWilsonLlT97l has become a major computational tool to 
study quantum impurity problems. Its advantages are 
that it is (i) non-perturbative, (ii) can handle exponen- 
tially small energy scales with unprecedented accuracy 
and (iii) does not suffer from any principle limitations 
regarding the parameter space of the model. In addition, 
through the inspection of the flow of the spectrum of the 
Hamiltonian under the renormalization group transfor- 
mation, the method provides direct very detailed infor- 
mation about the structure of the low-energy spectrum. 

Orig inally, its application w as limite d to stat ic quan - 
tities (iKrishnamurthv et all Il980albl: IWilsonl. Il975|) . 
but the calculation of dynamical properties is also 
possible (for a comprehensi ve overview of the early 
work see e.g. iHewso 1 I1993D an d thus its application 
to th e DMFT i mpurity problem ijSakai and Kuramotcl 
119941: IShimizu a nd Sakai, 1995). The direct calcu- 
lation of the one- particle self-energy from the NRG 
l|Bulla et all Il998|) triggered a variety of applications 
within the DMFT framework jBuiiaL Il999t iBulla et all 
200lt iPriisc hke et al'.. "2000: iPruschke and Zitzleii l2003t 
Zitzler et aL . .2002l .200,'1) at both T = and finite tem- 
peratures. 








Figure 11 Schematic setup of the NRG. The energy axis is 
partitioned into exponentially decreasing intervals. 

Given this tremendous success, the obvious question is 
if one can extend the NRG to become useful for quan- 
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turn cluster problems. Let us therefore briefly repeat the 
basic steps in the NRG scheme. The procedure starts by 
picking a number A > 1 and an interval on the energy 
axis that contains the support of the density of states of 
the bath degrees of freedom. As depicted in Fig. II II this 
interval is then partitioned into exponentially decreas- 
ing intervals ]A~("+^), A~"] for e > with n = 0, 1, . . .; 
likewise for e < 0. Within each interval, a Fourier decom- 
position of the bath operators is performed; the funda- 
mental approximation of the NRG consists of neglecting 
all Fourier components except for the homogeneous one. 
After tri-diagonalization, the resulting Hamiltonian reads 

i^QI ~ iflmp+ (179) 

OO 

(7 n—1 

where i?imp denotes the Hamiltonian of the "impurity" 

and the operators Unl represent the bath degrees of 
freedom. Note that the Hamiltonian l|179|l represents 
a semi- i nfinit e chain with the impurity at its left end. 
IWilsonI l|l97,^ showed that the quantities t„ cx A~"/^. If 
tn = y/A^tn and e„ = ^/A"en, one can cast the Hamil- 
tonian into a recurrency form 

Hn+1 = VAHn +^eN+ia^r^+i„aj^_^^^ (180) 

(7 

a 

Given the eigenvalues and eigenvectors of B.N-, -f^Af+i can 
be constructed straightforwardly. Since at each of these 
steps the dimension of the Hilbert space increases by a 
factor of four, the practical use would be limited even 
more severely as conventional exact diagonalization, be- 
cause the construction of Hn+i requires the knowledge 
of the whole set of eigenvectors and eigenvalues of Hn. 
However, the multiplication of Hn with -n/A > 1 expands 
the bandwidth of the spectrum and, because = 0(1), 
the low energy properties of -ff^v+i are determined by a 
restricted set of low-lying states of Hn only. This obser- 
vation is put into a practical computational scheme by 
the following algorithm: (i) At step TV > {N = cor- 
responds to -ffimp), diagonalize Hn and calculate all in- 
teresting local properties for that particular chain length, 
(ii) Use a suitable number A^^nrq of the lowest eigenstates 
of Hn to construct the next Hamiltonian Hn+i accord- 
ing to H18()(l . (iii) Continue with (i) until the desired 
accuracy for the ground state is reached. Note that step 
(ii) ensures that, no matter how long the chain becomes, 
the dimension of the Hamilton matrix to diagonalize can 
be fixed to a manageable number. 

As was discussed extensively in the literatur e ( Hewson, 
Il993t iKrishnamurthv et al l ll98nalbl:IWilsonill97a . the 
knowledge of the local properties at chain length TV pro- 
vides a means to calculate physical interesting quantities 
(thermodynamics and dynamics) at energy respectively 



temperature scales cx V . This exponential decrease 
of successive energy scales explains why the NRG is suit- 
able for studies of problems with extremely small dynam- 
ical energy scales. This however comes at the pr i ce of a 
loss of accuracy in high energy features l)HewsonLll993|) . 

Up to now the discussion was restricted to spin degen- 
eracy only. However, it is obvious that additional degrees 
of freedom do not influence the general lines of argument. 
In fact, for a problem with L > 1 internal degrees of free- 
dom in addition to the spin (like cluster sites or several 
orbitals per site), the result (|179|l acquires the form 

Hqi w i?imp + (181) 

L oo 
1=1 a n=l 

where iJjmp is again of arbitrary structure. Likewise, 
the recurrency relation (|18U|) can be set up and the al- 
gorithm extended. However, a simple example demon- 
strates that for L > 1 the technique can easily become 
useless. For a typical application to the single impu- 
rity Anderson model (SIAM) one chooses a A = 2 and 
A'nrg = 1000. This is sufhcient to obtain very accu- 
rate results for all relevant physical quantities. Let us 
now consider the next step, i.e. an orbitally degenerate 
problem with L = 2. Without coupling between the or- 
bitals, this corresponds to two independent SIAM we try 
to solve in a single calculation. Obviously, to obtain the 
same accuracy as in the true single impurity case, one 
needs at least TVnrg — 1000^ or A = 2^ (for a more de- 
tailed discussion of the issue of the accuracy of the NRG 
see Paula et al. (1923} and references therein). 

While A^NRG = 1000^ is beyond all numerical pos- 
sibilities, a strongly increasing A introduces a huge 
loss in accuracy, b o th a t high and low energies 
(jol iveira and OliveiraL Il994t) . As long as L = 2 the 
problem can be partially compensate d by respecting the 
additional symmetries in t he system llSaka i et a|J . |l989l) 
or other numerical tricks dol iveira and Oliveiral |1994). 
However, as a generic tool to solve the effective quantum 
cluster problem arising in embedded cluster techniques 
much larger values of L must be accessible. Thus, while 
for a study of qualitative effects of non-local correlations 
with Nc — 2 the NRG is still applicable (see results in 
section llV.D.l|l . the preceding discussion makes it clear 
that already for Nc= A the method is currently not fea- 
sible. 



IV. APPLICATIONS TO STRONGLY CORRELATED 
MODELS 

In this section we review the application of various 
quantum cluster approaches to a selection of standard 
models of strongly correlated electron systems. We put 
special emphasis on the capabilities and advantages of 
these techniques over both finite system simulations and 
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DMFT. In Sec. IIV.A1 we show that quantum cluster ap- 
proaches are complementary to FSS, i.e. that taken to- 
gether the information obtained from both techniques 
can yield conclusive results. The effects of non-local cor- 
relations on single- and two-particle spectra as well as 
phase diagrams are emphasized throughout the discus- 
sions of the Falicov-Kimball model in Sec. IIV.BI 



the ID 

Hubbard model in Sec. lIV.Cl and the 2D Hubbard model 
in Sec. IIV.DI Due to space restrictions we have to omit 
recent applications of the CPT and DCA to electron- 
phonon systems and re fer t he reader to t he articles by 
iHohenadler et al\ l)2003|) and lHagud l)2003|) . respectively. 



A. Complementarity of finite size and quantum cluster 
simulations 

FSS and quantum cluster approaches yield exact solu- 
tions in the infinite cluster size limit. At finite cluster size 
Nc quantum cluster approaches differ from FSS by the 
coupling to a self-consistent dynamic host. At cluster size 
Nc = 1 this difference is most pronounced: While FSS 
reduce to the atomic limit, quantum cluster approaches 
reduce to the DMFT, i.e. a highly non-trivial approxi- 
mation to the infinite size lattice. Thus it is instructive 
to compare quantum cluster results with those obtained 
from FSS systematically as a function of cluster size. 
This has been done in the half-filled 2D Hubbard model 
using the DCA/QMC algorithm, i.e. DCA combined 
wit h QMC as a cluster so l ver, b y lHuscroft et al\ ll2001 



an d iMoukouri and Jarrelll l|200ll) and with DC A/FLEX 
bv lArvanpour et al\ l|2003a l. 

In the 2D Hubbard model at half-filling, the antiferro- 
magnetic correlation length C increases with decreasing 
temperature and diverges at T = 0. In FSS, i.e. in 
a finite size lattice with periodic boundary conditions, 
the system freezes when the correlation length exceeds 
the system size and a gap to excitations opens. As the 
system size is increased this tendency is reduced. Corre- 
lation induced gaps are thus generally overestimated in 
FSS for smaller clusters since the system is artificially 
closer to criticality. In contrast, in quantum cluster ap- 
proaches the system is in the thermodynamic limit with 
correlations restricted to the cluster size. Hence the sys- 
tem never freezes. As the cluster size is increased, longer 
ranged correlations are progressively included. The ef- 
fects of correlations therefore increase with cluster size. 

This behavior is illustrated in Fig. ^1 where we 
reproduce the results obtained from DCA/QMC and 
FSS/QMC simulations for the imaginary time Green 
function G(k, t) at the Fermi wave- vector k = (tt, 0) 
for different cluster sizes Nc- This quantity has a more 
rapid decay from its maximum at t = /3/2 when correla- 
tion effects are stronger so the gap is more pronounced. 
With increasing cluster size, the DCA and FSS results 
for G'(k, r) converge from opposite directions. Consistent 
with the expectation that correlation effects are overesti- 
mated in FSS but underestimated in the DCA, the decay 




Figure 12 The imaginary time Green function G(k, r) at 
k = (tt, 0) in the half-filled 2D Hubbard model, cal- 
culated with finite size QMC (filled symbols) and DCA 
(emp ty symbols) when U = A At and 13 = 4/t. Taken 
from [ Moukouri and Jarrell ll200lf) . See also lAryanpour et ali 
tOOM) . 
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Figure 13 The imaginary time Green function G(k,r) at k = 
(tt, 0),t — (3/2 in the half-filled 2D Hubbard model, calculated 
with FSS/FLEX (circles) and DCA/FLEX (squares) when 
U — 1.57t for various cluster sizes Nc = Lx L). Both methods 
scale as and conv e rge to a single value as L — > oo. Taken 
from lAryanpour et al\ l2003al) . 



is stronger for smaller system sizes in FSS, while the DCA 
results show the opposite behavior. Since the two tech- 
niques become identical in the infinite cluster size limit, 
the exact G(k, t) curve is bracketed by the FSS and DCA 
curves. For methods like the FLEX (see Sec. lIII.B.2|) for 
which calculations with much larger clusters are feasi- 
ble, scaling of the results to the infinite system is possi- 
ble (see Fig. 1121). The extrapolations of FSS/FLEX and 
DCA/FLEX to the infinite system (l/L^ 0) coincide 
within numerical uncertainties, thus allowing the deter- 
mination of the infinite lattice result with unprecedented 
accuracy. 

This complementarity is also seen in results for the 
spectral gap in the 2D half- filled Hubbard model (see 
iHuscroft et all l200ll and lArvanpour et al 1 l2n03aH . In 
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the DCA the gaps converge from small to large as the 
cluster size increases, while the converse occurs in FSS. 
Although we have only shown results of the DCA, we ex- 
pect the CDMFT to show similar size dependence, since 
DCA and CDMFT share the same nature of the approx- 
imation (see discussion in Sec. III. ('(I . Results obtained 
with the CPT algorithm however can be viewed to some 
extent as a periodic continuation of FSS. Thus it is an 
open question if the CPT method shows similar comple- 
mentarity. 



B. 2D Falicov-Kimball model 

The usefulness of the discussed cluster theories was 
first demons trated by the application o f the DCA to the 
2D FKM bv lHettler et al ] l|200(lll998|) . While the FKM 
is a particularly gentle test bed for novel approaches, it 
allows the study of the effects of non-local fluctuations. 

The FKM can be considered as a simplified Hubbard 
model in which one spin species has zero hopping ampli- 
tude. The Hamiltonian reads 



H = -tY, 4d,+Uj2{r 



(182) 



with nf — dld^ and n{ — fjf- . For a 2D square lattice 
the bandwidth of the noninteracting system is W = 8t. 
At half filling and D > 2 the system has a second or- 
der phase transition from a homogeneous high tempera- 
ture phase to a charge density wa ve with ordering vec- 
tor Q = (tt, tt) for any nonzero U ijBrandt and Schmidd 
[1986, 1987). The universahty class is that of the 2D 
Ising mo del, the strong coupling limit U/t 1 of the 
FKM. In lHettler et aZ.I l)2000|) the FKM was evaluated 
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Figure 14 DOS of the 2D FKM for a 4 x 4-cluster at vari- 
ous temperatures calculated with DCA. The DOS develops 
a pseudogap as the temperature approaches Tc « 0.189t 
{U = 4t) due to the non-local CDW fluctuations present in 
the DCA (iVc > 1). In the DMFT {Nc = 1), t here is no T- 
dependence of the DOS above Tc- Figure from lHettler et all 

Eood). 



within the DCA by a combination of QMC methods and 
exact enumeration for small clusters. Since the DMFT 
is a single-site theory {Nc — 1), it yields an unphysical 
temperature-ind ependent DOS of the mobile d-electrons 
ijBrandt and Mielsch. 1989) due to the constraint of half 
filling (one electron per site of either the d or f vari- 
ety). In a cluster theory with TVc > 1 this artifact is 
absent, since a redistribution of Boltzmann weight with 
temperature is possible among the various configurations 
of d or / electrons on the cluster sites, while maintain- 
ing the condition of half filling on average. This tem- 
perature dependence is demonstrated in Fig. 1141 where 
a pseudogap develops in the local DOS with decreasing 
temperature. This pseudogap can be interpreted as a 
precursor of the eventual transition to a CDW phase, 
which features a full gap at the Fermi level. In addition 
to the gap, there is fine structure in the DOS, related to 
an exchange energy Jeff - This is better observed in the 
momentum resolved spectral function ( see Hcttlcr et 
I2OOOI: iLaad and van den Bosschd l200of) . 

As stated above, the FKM has an instability to a phase 
with CDW order. As discussed in Sec. II. Bl embedded 
cluster theories exhibit phase transitions at some tem- 
perature that, due to their residual mean-field charac- 
ter, lies above the exact Tc of the infinite system. As 
the cluster size increases, one expects the effect of the 
mean-field to decrease, leading to a decreasing Tc with 
increasing cluster size. The Nc dependence of the tran- 
sition temperature Tc is shown in Fig. 1151 together with 
a comparison with the Tc obt ained from FSS methods 
l|de Vries et all Il993albl. Il994j) and Tc of the 2D Ising 
model with exchange coupling J = l/(2J7).The extrap- 
olated cluster results agree with the FSS estimates and, 
for large values of U , also with the results obtained from 
the 2D Ising model. For smaller U however, charge fiuc- 
tuations begin to play a larger role, suppressing the Tc 
compared to that of the Ising model which lacks charge 
fluctuations. 

The effect of different boundary conditions in the DCA 
cluster is illustrated in the inset^. In small clusters the 
effect is strong, but already in a 6 x 6- cluster the bulk of 
the cluster dominates and the boundaries play a minor 
role. 

When Nc — 1 only charge fluctuations with an en- 
ergy scale U are present. In a cluster theory the non- 
local "spin" fluctuations with an effective energy scale 
Jeff (oc 1/2U for large U) must also be observable in 
thermodynami c quantities li ke the entropy and the spe- 
cific heat. In .Hettler et al\ l)2000() the entropy and spe- 
cific heat was computed from the tota l energy of the 
cluster via a maximum-entropy method llHuscroft et alV . 
For better comparison, the calculations were per- 



^ According to the derivation of the DCA formalism (see Sees. III. XI 
111. Bl and Fig.|2J the DCA cluster has periodic boundary condi- 
tions. By shifting the set of cluster K-points however, different 
boundary conditions can be simulated. 
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Figure 15 Tc in the 2D FKM as a function of the inverse Unear 
cluster size for U/t = 8,12,16 calculated with DCA. The 
Ising limit, and FSS estimates of Tc are shown for comparison. 
The inset shows that the influence of the cluster boundary 
conditions on Tc disappears rapidly with increasing cluster 
size. Figure from lHettler et al. (.2000.') . 

formed in the uniform phase, even at temperatures below 
the CDW ordering Tc by not allowing for the symmetry 
breaking. The results for a 2 x 2- cluster are shown in 
Fig. ^1 where the ratio of specific heat C over the tem- 
perature T is plotted for the DMFT and the DCA on a 
2x2- cluster. The appearance of the a second peak at 
lower temperature is a clear indication for the additional 
non-local fluctuations present on the cluster. The effect 
on the entropy is also strong, as shown in the inset. 

C. ID Hubbard model 

In this section we discuss the application of quantum 
cluster approaches to the ID Hubbard model (in usual 
notation) 
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Figure 16 Specific heat versus temperature for single site and 
2x2- clusters calculated with DCA and exact enumeration 
when JJ = 8t. For Nc — 1, there is a single peak with in- 
tegrated weight ln(2) associated with the suppression of lo- 
cal charge fluctuations. For Nc = 4, there is an additional 
peak at lower temperatures associated with critical fluctua- 
tions near the charge ordering transition temperature. Tc for 
A'^c = 4 is indicated by an arrow. In the inset the entropy 
S(T) = frT dT'^^ is shown in units of ln(2). Figure from 
iHettler ei all J200d) . 



the results with those obtained from the DMFT and 
DMRG approaches. As an example of this study we show 
in Fig. El a comparison of the exact result with that ob- 
tained from DMFT (referred to as LISA for "Local Im- 
purity Self-consistent Approximation") and CDMFT for 
the single-particle spectral gap A{U) as a function of the 
on-site Coulomb repulsion U in the half-filled case. The 
total number of sites in the effective cluster model in- 
cluding cluster and bath sites was fixed to six in the ED 
approach (see Sec. IIII.C.2|l . 



n = -tY,i4+laC^. + h.c.) + UY^ n,Tn,| , (183) 

i.fT i 

which provides a non-trivial test ground for these tech- 
niques. In ID quantum fluctuations are stronger than in 
higher dimensions. Hence quantum cluster approaches 
which cut off correlations beyond the length-scale set 
by the cluster size are expected to be less efficient than 
in higher dimensions. Therefore if quantum cluster ap- 
proaches accurately describe the physics in ID, they are 
highly likely to capture the physics in 2D and 3D. In 
addition, since the exact ground-state of the ID Hub- 
bard model is known from the Bethe ansatz solution 
iLieb arid Wu, 1968), a quantitative comparison of cer- 
tain static quantities is possible. Fairly reliable results 
for dynamical quantities can be obtained from the den- 
si ty matrix ren o rmali zation method. 

liBolech et al\ l|2003(l applied the CDMFT/ED method 
to the ID Hubbard model and systematically compared 




U/l 



Figure 17 Spectral gap in the half-filled ID Hubbard model as 
a function of U/t calculated with DMFT (LISA) and CDMFT 
for a 2-site cluster, Nc = 2 c ompared to the exact result. 
Taken from iBolech al\ i2003ri . 
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The behavior of the cluster results (CDMFT for = 
2) is fundamentally different from the single-impurity re- 
sult (LISA, Nc = 1): When Nc = 1, the spectral gap is 
reduced to values much smaller than the exact value and 
even vanishes at finite U. The existence of this Mott- 
transition in ID for A^c = 1 is a well known artifact 
of the DMFT method, reflecting the physics of infinite 
dimensions wh e re the Mott-transition is indeed present 
l|Georges et all Il996fl . In contrast, the Nc = 2 result 
follows the exact result quite accurately yielding an in- 
sulating solution for all finite values of [/ > in agree- 
ment with the well known physics of the ID model. This 
excellent reproduction of the exact gap for A'c = 2 is 
quite encouraging since larger cluster sizes should pro- 
duce even better results. And indeed, a comparison of 
the nearest neighbor Green function with DMRG results 
sho ws systematic impro vements with increasing cluster 
size l)Bolech et all \200^ . 

A related fundamental feature in ID correlated sys- 
tems is the breakdown of the Fermi-liquid picture be- 
cause of spin-charge separation as described in the con- 
cept of Luttinger-liquids (see e.g. Voit. 1 99^ . In a Fermi- 
liquid, the spectral weight A(k, to) is centered around a 
single quasiparticle peak at w = ek, while in a Luttinger- 
liquid, A(k, uj) is distributed between two singularities 
associated respectivel y with spin- and charg e-excitations 
(spinous and holons). ISenechal et a,l\ l|2nnffl have calcu- 
lated A(k, uj) in the ID Hubbard model using the CPT 
formalism. Fig. 1181 shows a comparison of this quantity 
at the Fermi wave- vector fc = 7r/2 as calculated by ordi- 
nary ED and its infinite lattice extension within the CPT 
method for various cluster sizes A^ when U = 8t (bottom) 
and U = m (top). 
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Figure 18 The spectral function A{k,ijj) at fc = tt/2 in the 
half-filled ID Hubbard model when U = 16t (top) and U = 8t 
(bottom), calculated with ordinary ED a nd with CPT fo r 
cluster sizes A = 4, 8, 12. Taken from .Senechal et al\ i2000tl . 

While no sign of spin-charge separation is seen in 
the pure ED results, the CPT method reveals the two 



branches of the spectral weight indicative of spin-charge 
separation. The two-peak structure resolves more clearly 
as the cluster size increases. Since propagation between 
clusters requires the spinon and holon to recombine, spin- 
charge separation can only exist on length- and time- 
scales limited by the cluster size. Consequently it takes 
fairly large clusters to clearly resolve this feature. 
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Figure 19 The spectral function A(k,Lj) in the half-filled ID 
Hubbard model when U = At, calculat ed with CPT for Nc = 
12. Taken from (iSenechal et aLll2000D . 



The analysis of the fc-dependent spectral function in 
Fig. El reveals the dispersion of the two branches. The 
spinon- and holon-branches show different dispersions 
and can be clearly identified together with the gap at 
the chemical potential at k — Tr/2. These T — simu- 
lat ions were extended to finite low temperatures T > 
bv lAichhorn et all l)2003t ) who implemented a novel low- 
temperature Lanzcos algorithm which connects the ex- 
act ground-state Lanczos method with the established fi- 
nite t emperature Lanzcos method (jjaklic and Prelovsekl 
l200Cl|) . They showed that the spin-charge separation ob- 
served in the T — results persists at finite, low temper- 
atures. 

As exemplified by these studies, quantum cluster ap- 
proaches have been proven very useful to explore the 
complex behavior of the ID Hubbard model. Even stud- 
ies with small cluster sizes are consistent with well-known 
results such as the opening of a Mott gap at half-filling 
for all ?7 > and the existence of spin-charge separation 
reflected in the two branches in the electronic dispersion. 
This success in the ID case shows great promise for ap- 
plications of cluster methods in 2D or 3D where cluster 
methods are expected to be even more efficient since cor- 
relations are less pronounced. 
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D. 2D Hubbard model 

The interest in the 2D Hubbard model (in usual nota- 
tion) 

W = E ^^AaC.a + UJ2 "»T«a (184) 
ij^(T i 

has been revived recently in particular since it is believed 
to capture the physics of the superconducting planes in 
high-t e mperature super conductors (HTSC) (Andersoiij, 
Il997bl: IZhang and Ric5 Il988l [1990). This section re- 
views various applications of quantum cluster approaches 
to the 2D Hubbard model at half-filling and at finite dop- 
ing, including results for a possible Mott-Hubbard tran- 
sition, antiferromagnetism and its precursors, pseudogap 
phenomena and superconductivity. 



1. Metal-insulator transition 

The study of a possible metal-insulator transition 
(MIT) in the 2D Hubbard model at half-filling {e^ = 
Ui = ~U/2) is under active research. This problem 
was studied in the unfrustrated model, i.e. with only 
nearest neigh bor hopping ta = enSa — tS(^ij^ within 
DCA/QMC bv lMoukouri and Jarreifl)200l|) and with the 
two-site composite operator method (see Sec. IH.H|) by 
IStanescu and Phillips! l)200lj) . The frustrated case (addi- 
tional nex t-nearest neighbor hop ping t'(5((y))) was inves- 
tigated bv lParcollet et al\ l)2003|) . 

Numerical calculations have shown that the ground 
state of the unfrustrated model is an antiferromagnetic 
insulator with the Neel temperature Tn = constrained 
by the Mermin- Wagner theorem. A spectral gap hence 
exists at r = 0. However the central question about the 
origin of the gap and its relation to antiferromagnetic or- 
dering is less understood. Is the gap a direct consequence 
of the antiferromagnetic ordering at T = or does it arise 
from strong correlations at higher temperatures? 

To appreciate the significance of this issue it is impor- 
tant to understand the fundamental difference between 
antiferromagnetic insulators, i.e. insulators due to mag- 
netic ordering and Mott insulators, i.e. insulators due to 
electronic correlations. Antiferromagnetic insulators re- 
sult from the doubling of the unit cell in the ordered state 
and are therefore adiabatically connected to band insu- 
lators which have an even number of electrons per unit 
cell. In contrast, paramagnetic Mott insulators have an 
odd number of electrons per unit cell and are therefore 
fundamentally different from band insulators. 

At strong coupling {U ^ W) the situation is well un- 
derstood: A charge gap of order U develops in the spec- 
trum below temperatures T k U due entirely to strong 
electronic correlations. The spins are coupled by the ex- 
change interaction J = /U and govern the low-energy 
physics. As a result spin and charge are separated. Sys- 
tems in this regime are hence Mott insulators and the 



antiferromagnetic ordering at T = is merely the result 
of the Mott-transition at higher temperatures. 

Different scenarios however exist for the weak-coupling 
regime {U <^ W): In the weak coupling point of view, a 
spin density wave forms at T = due to the nesting of the 
Fermi surface and leads to the doubling of the unit cell. 
Hence the gap in the spectrum is a direct consequence 
of the antiferromagnetism at T = 0. This perturbative 
point of view is referred to as Slater mechanism. It is in 
contrast to the second opinion due to'AndersoiJ l)l997alb|) 
who has argued that the 2D half-filled Hubbard model is 
always in the strong-coupling regime, so that a Mott gap 
is present for all ?7 > as in ID (see Sec. IIV.CJ)) . As the 
temperature decreases, local moments develop because of 
the opening of the Mott gap which then order at T = 0. 
Thus the antiferromagnet at T = is a consequence of 
the Mott transition. 

The MIT in the half-filled Hubbard model has been 
extensively studied within the DMFT (for a review see 
iGeorges et a^J . IT99^ . In the DMFT one can easily dis- 
entangle the effects leading to antiferromagnetic and 
Mott gaps. The DMFT equations for the paramag- 
netic state of the bipartite Bethe lattice are identical 
with the equations of the fully frustrated infinite dimen- 
sional model (Gcorges_et a/., 1996). This justifies the 
study of the paramagnetic solution within the antifer- 
romagnetic phase of the unfrustrated model which shows 
a first order Mott MIT end i ng at a finite temperature 
critical point ijGeorges et a/.l . ll996j) . Although this justi- 
fication does not hold for Nc > 1, one can still study 
the paramagnetic solution by enforcing the spin sym- 
metry and hence avoiding the opening of a full spectral 
gap due entirely to magnetic orderin g. Following this 
approach. iMoukouri and Jarrelll l)200l|) studied the MIT 
in the unfrustrated 2D half-filled Hubbard model using 
DCA/QMC systematically as a function of cluster size 
Nc- The MIT can be identified by analyzing the behav- 
ior of the double occupancy (D) = (n|n|). This quantity 
is shown in Fig. |^ for different values of the Coulomb 
repulsion U and cluster size Nc- 

When Nc — I the double occupancy displays evidence 
for a MIT when U is of the order of the bandwidth 
W = 8t = 2: (D) is monotonically increasing with 
temperature when [/ = 3, but displays a minimum for 
U = 0.5 and 1 indicating the emergence of quasiparti- 
cle states at the chemical potential at low temperatures. 
When A'c > 1 the situation is radically different: In the 
strong coupling regime {U = 3), local fiuctuations dom- 
inate and (D) is essentially independent of Nc except at 
very low temperatures. In contrast, in the weak coupling 
regime, the minimum found for Nc = 1 flattens progres- 
sively as Nc increases from 8 to 16. When iVc > 36 a 
downturn in (D) appears at low temperatures. By open- 
ing a gap and hence localizing the moments, the system 
can gain free energy by taking advantage of the short- 
ranged magnetic correlations. 

The results for the DOS shown in Fig. |^ support 
the evidence from the double occupancy. The quasi- 
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Figure 20 Double occupancy as a function of temperature 
for different values of U in the half-filled 2D Hubbard model 
calculated with DCA/QMC for various cluster sizes Nc- Here, 
t = 0.25, so the lin es for U = 0.5, 1,2 co r respo nd to U/t — 
2,4, 12. Taken from iMoukouri and JarTelH J2001D . 
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Figure 21 DOS when U/t = 4 for different temperatures in 
the half-filled 2D Hubbard model calculated with DCA/QMC 
for va rious cluster sizes Nc. Taken from iMoukouri and Jarrelli 
ll200lD. 
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Figure 22 Tiling of the first Brillouin zone for Nc = 2. 



The appropriate tiling of the Brillouin zone is shown in 
Fig. 122 and the resulting coarse-grained spectral func- 
tions 74(K,w) = — i3mG(K, ly) for the two cluster K- 
points are shown in Fig. 1231 Even for small values of U 
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Figure 23 DCA coarse-grained spectral functions for different 
values of U obtained from Nc — 2 DCA/NRG calculations at 
T — 0. Apparently, no finite critical value of U for a conven- 
tional Mott-Hubbard transition seems to exist. In contrast, 
for all 17 > a gap exists at the Fermi energy. 



particle Kondo-like resonance at the chemical potential 
for Nc = 1 is destroyed by non-local correlations when 
iVc > 1. As A^c increases, a gap opens at the chemical 
potential and the Hubbard sidebands become more pro- 
nounced consistent with the suppression of (D). Given 
the fact that DCA always underestimates correlation in- 
duced spectral gaps (see Sec. IIV.A|) . these simulations 
indicate the absence of a weak-coupling regime in the 
unfrustrated 2D Hubbard model at half-filling consistent 
with Anderson's point of view. 

Another interesting aspect is whether this result 
changes at zero temperature, T = and the MIT pre- 
dicted by DMFT returns. As evidence that even at 
T = one must expect a gap in the spectra for any 
U > 0, we show here first DCA results at T = ob- 
tained with Wilson's NRG for a cluster size of Nc — 2. 



a well defined gap exists in the spectrum at the Fermi 
energy. Note that all calculations were done in the para- 
magnetic phase, i.e. the concept of a Slater insulator does 
not apply here. The gap quickly increases with increasing 
U and at the same time the system gains more spectral 
weight in the incoherent parts of the spectrum, i.e. begins 
to resemble what one expects from Mott localized states. 

The D CA results wer e confirmed in a similar cluster 
study bv lStanescu and PhiUips (,2001) using the two-site 
composite operator method discussed in Sec. lH.HI Figure 
1241 shows the results of this study for the temperature de- 
pendence of the value of the DOS at the chemical poten- 
tial, p(0) in the ID and 2D models for small U = 2t. As 
a consequence of the shape of the non-interacting DOS, 
the 2D result for p{0) is enhanced over the ID result. 
However, in both ID and 2D, p(0) falls to zero as the 
temperature decreases indicating the absence of a metal- 
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Figure 24 DOS as a function of temperature at the chemical 
potential in the ID and 2D half-filled Hubbard model when 
U = 2t calculated wi th a two-site composite opera tor approx- 
imation. Taken from lStanescu and Philli'p^ ll200lf) . 



lie state at half-filling even for small U. 

To elucidate the role of antiferromagnetic correlations 
in the opening of the Mott gap, the frustrated Hub- 
bard model may be studied. In the t-t' Hubbard model, 
a next-nearest neighbor hopping t' between sub-lattices 
strongly frustrates a ntiferromagnetic corre lations. This 
model was studied bv'Parc ollet et al\ l)2003(l for t' = t on 
a 2 X 2 cluster using CDMFT/QMC. We reproduce their 
results for the U dependence of the double occupancy 
docc = l/4X]i=i("-iT"'ji) different temperatures in 
Fig. [23 Similar to the behavior found in DMFT, docc 
displays a downturn at a critical value Uc, indicating a 
transition from a metallic to an insulating state. An in- 
spection of the spectral weight A(k, 0) at the chemical 
potential reveals that the gap opens first in the region 
around k = (tt, 0) (,Parcollet et al. . 2003^ . These results 
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Figure 25 Double occupancy as a function of U in the frus- 
trated 2D Hubbard model with t' = t for different tempera- 
tures (from top to bottom) T/t = 1/5, 4/30. 1/10, 1/11 calcu - 
lated with CDMFT/QMC. Taken from lParcollet et aZj (120031) . 

were substantiated by a CPT study of the frustrated 



2D Hubbard model l|Senechal and Tremblavl |2003|) for 

t' — —OAt. Although not the focus of this study, the 
results show further evidence of a Mott transition at a 
finite value of U in the filling dependence of the chemical 
potential. 

The existence of the Mott transition in the frustrated 
Hubbard model and its absence in the unfrustrated model 
seems to indicate that antiferromagnetic correlations play 
a key role in the opening of a Mott gap at weak coupling. 
Since the opening of the gap occurs in the paramagnetic 
solution, it cannot be attributed to the existence of an- 
tiferromagnetic ordering. Thus the conclusion reached 
from these cluster studies is that a symbiosis of local mo- 
ment formation and short-ranged antiferromagnetic cor- 
relations cause the gap to o pen at finite temperatures 
l|Moukouri and Jarrel]Ll200l|) . 



2. Antiferromagnetism and precursors 

If the simulations are performed without enforcing the 
spin symmetry or frustrating the lattice, the system is 
given the possibility to transform to a state with antifer- 
romagnetic order. Since the system is two-dimensional, 
we know from the Mermin- Wagner theorem that the true 
Neel temperature is necessa r ily zer o. As found in in- 
finite dimensions ll Jarrell et all \199^ . the DMFT pre- 
dicts a finite temperature transition even in 2D due to 
its mean-field character. As discussed in Sec. II. Bl clus- 
ter approaches restore non-local fluctuations and thus are 
expected to drive the Neel temperature systematically to 
zero as the cluster size increases. 

As discussed in Sees. III.FI and III. El phase transitions 
can be identified from the disordered (here: paramag- 
netic) state by calculating the corresponding susceptibil- 
ity, or from the ordered state by computing the order 
parameter. The calculation of order parameters is exem- 
plified in Fig. 1^ where we plot the DCA/NCA result for 
the sub-lattice magnetization m — \/N'Y^^ ^e^^'^^ariia 
(see Eq. (IllOll ) as a function of the reduced temperature 
t = T/Tn in the 2D half-filled Hubbard model for the 
cluster sizes Nc = 1 and Nc — ^ when U = At. The 
Nc — 4: Neel temperature Tn — 0.208t is reduced from 
the Nc — 1 result Tn = 0.304t and the order parame- 
ter is strongly suppressed. As expected, non-local spin 
fiuctuations suppress antiferromagnetism. 

Figure 123 shows the DCA/QMC result for the temper- 
ature dependence of the inverse antiferromagnetic sus- 
ceptibility 1/xaf B,t U = 6t for various cluster sizes Nc 
in the paramagnetic state. At high temperatures the 
susceptibility is independent of Nc due to the lack of 
non-local fluctuations. In contrast to ESS, the low tem- 
perature susceptibility diverges at T = Tn indicating the 
instability to the antiferromagnetic state. When Nc ^ 1 
the susceptibility diverges with a critical exponent 7 « 1 
as expected for a mean-field theory. Consistent with the 
NCA results the susceptibility diverges at lower temper- 
atures when Nc > 1 with larger exponents indicative of 
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Figure 26 Sublattice magnetization as a function of temper- 
ature in the half-filled 2D Hubbard model calculated with 
DCA/NCA for cluster sizes Nc = 1,4 when C/ = 4t. Ins et: 
Neel-temperature versus doping. Taken from lMaieil i2001f) . 



fluctuation effects. However as discussed in Sec. lI.BI these 
critical exponents reflect the behavior at intermediate 
temperatures. Very close to the transition, there must 
a region of mean-field behavior. However, this region is 
very difficult to resolve with DCA/QMC, due to numer- 
ical noise, which is especially large near the transition. 



tuations which are responsible for suppressing the an- 
tiferromagnetism result in precursors of the antiferro- 
magnetic phase transition. The onset of antiferromag- 
netic correlations on short time- and length-scales may 
be signaled by a pseudogap in the DOS as a precur- 
sor to the antiferromagn e tic ga p. This was predicted 
by iKampf and Schrieffed l)l99nl) using a phenomeno- 
logical ansatz for the weak coupling Hubbard model 
based on the presence of strong antiferromagnetic spin- 
fluctuations. On a microscopic level, this question has 
b een addressed by FSS/Q MC in the 2D Hubbard model 
in lVekic and White! ill QQ.ll) and lCreffield oil (119951) and 
by ap prox imate many-body tech niques in iDeisz et all 
l | l 99 6 t ) and lMoukouri et all l|l999|) . But the results have 
been inconclusive as to the existence of a pseudogap at 
low temperatures, due to the limitations of these tech- 
niques. 

Within quantum cluster appro aches the p seudogap 
phenomenon was first studied by iMaier etal. (2000^ 
within the DCA/NCA formalism. In contrast to A^c = 1 
where a Kondo-Iike quasiparticle peak emerges at the 
chemical potential as the temperature is decreased (rem- 
inisce nt of the D — oo DMFT result l)Ceorges et all . 
Il996j) '). a pseudogap was found when non-local correla- 
tions were included in Nc = 4 simulations. For larger 
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Figure 27 Inverse antiferromagnetic susceptibility versus tem- 
perature in the half-filled 2D Hubbard model calculated with 
DCA/QMC for various cluster sizes Nc when U — 6t. The 
Unes are fits to the function (T — Tn)''. Inset: Corresponding 
Neel-tempcraturcs as a function of the cluster size. Taken 
from JarrcU et al. ( 2001b). 
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Figure 28 Spectral function j4(k,a;) and the real 5ReE(k,tj) 
and imaginary SjmE(k, uj) parts of the self-energy for vari- 
ous temperatures at k = (tt, 0) in the 2D half-filled Hubbard 
model for U = 5.2t calculated with DCA/QMC for a 64-site 
cluster {Nc = 64). Taken from .Huscroft e t al. (2001). 



As shown in the inset, the transition temperatures fall 
very slowly with the cluster size Nc, but the iVc = 4 
falls off the curve. As detailed in Sec. III. ('.41 fluctua- 
tion effects in clusters with linear size L = 2 are over- 
proportionally enhanced since its coordination number 
is reduced compared to the original system. Hence the 
Nc = 4: result does not fall on the curve, similar to the 
behavior seen in DCA studies of the FKM l)Hettler et all . 
12000') . 

The question arises of whether the same non-local fiuc- 



cluster sizes, the emergence of t he pseudogap i n the DOS 
was explored b y DCA/QMC in.Huscroft all ^QQl) and 
DCA/FLEX in lArvanpour et all (|2003a|) . Figure dis- 
plays the DCA/QMC results for the spectral function 
v4(k,a;) and the self-energy S(k,tt') at the Fermi wave- 
vector k — (tt, 0) in the paramagnetic state for a 64-site 
cluster {Nc — 64) at various temperatures. With decreas- 
ing temperature a pseudogap develops in A(k, uj) at the 
Fermi wave- vector k = (tt, 0). Simultaneously the slope 
of 3fJeE(k, 0) becomes positive at k = (7r,0) signahng 
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the appearance of two new solutions in the quasiparticle 
equation w — ek + M — 5ReI](k, u) — 0. In addition to the 
strongly damped solution at w = which is also present 
in the non-interacting system these two new quasiparticle 
solutions appear on both sides of a; = 0. A consequence 
of the antiferromagnetic order on short time- and length- 
scales, they can be viewed as precursors of the doubling 
of the unit cell in the antiferromagnetic state. The pseu- 
dogap is generated by the local minimum in 9mS(k,w) 
which signals the breakdown of Fermi-liquid behavior. 

By studying the syste m on a triangular lattice, 
llmai and Kawakamil l)2002() investigated the effects of 
frustration on the pseudogap in the half-filled 2D Hub- 
bard model using the DCA/NCA and DC A/FLEX ap- 
proaches. Figure EHl schematically illustrates the trian- 
gular lattice and the choice of cluster wave- vectors in the 
corresponding Brillouin zone for Nc = 4. 




Figure 29 Illustration of the triangular lattice with hopping 
amplitudes t and t' (a) and DCA coarse-graining cells (b) in 
the first Brillouin zone (dashed line) of the triangular lattice 
when Nc = 4. Cluster K-point s are i ndicated by the dots. 
Taken from llmai and Kawakaniil i2002l) . 

For t' = this setup corresponds to the unfrustrated 
system and the effects of frustration can be systemat- 
ically studied as t' is increased to its maximal value 
t' — t. Figure 1201 reproduces the results for the DOS 
and coarse-grained spectra A(K.,lu) for different values 
of the frustration t' . As the geometrical frustration in- 
creases from t' = to t' = t, antiferromagnetic spin fluc- 
tuations are suppressed. Consequntly the pseudogap in 
the unfrustrated system diminishes and a quasiparticle 
peak develops at the chemical potential. The change in 
the DOS mainly originates in the region in momentum 
space around K = (7r,7r/V3) where the Fermi surface is 
located. These results are thus consistent with an anti- 
ferromagnetic spin fluctuation driven pseudogap. 



3. Pseudogap at finite doping 

The properties of the Hubbard model away from half- 
filling are of great interest especially in the context of 
HTSC. Contrary to Fermi-liquid theory, low-energy spin- 
excitations in HTSC are suppressed at low temperatures 
as evidenced by Knight-shift experiments. Concomi- 
tantly, the Fermi surface is gapped along certain direc- 
tions in the Brillouin zone as indicated in ARPES experi- 
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Figure 30 The DOS (a) and coarse-grained single-particle 
spectral functions A{K.,lj) in the 2D half-filled frustrated 
Hubbard model for K = (0,0) (b), K = (7r,7r/y3) (c) and 
K = (0,27r/V3) (d) for various values of the frustration t' 
when U = 6t and T = O.Qt. Taken from .Imai and Kawakarnl 
(I2002D . 



ments. This pseudogap phenomenon^ has proven a great 
challenge for theories of strongly correlated systems. 

DMFT has provided great insight into the evolution 
of spectra in doped Mott-insulators. Exact results based 
on the self-consistent mapping onto an Anderson impu- 
rity model show that the system is a Fermi-liquid in 
the metallic state in the absence of symmetry break- 
ing below a coheren ce temperature reminiscen t of the 
Kondo-temperature l|Georges and Kotliailll99^ . Hence 
the spin-susceptibility becomes finite at low temperatures 
in contrast to the ex perimental results in und erdoped 
cuprates. Furthermore lMiiller-HartmannI l)l989a|) showed 
that because the self-energy is momentum independent, 
volume and shape of the Fermi surface are identical to 
the non-interacting Fermi-surface. Thus DMFT does not 
include the effects that lead to the emergence of a pseu- 
dogap in the spin- and quasiparticle spectrum and cluster 
extensions are necessary. 

Within quantum cluster approaches the pseudogap 
phenomenon in the doped 2D Hubbard model was stud- 
ied wi th DCA/NCA in IMaier et al\ feOOObI) and [Maie3 
1I2OOID. with DCA/OMC in I.Tarrell et al\ ll2n01ah md 
iMaier et al\ (|2002a), wit h the two-si t e compo site opera- 
tor method inTstanescu and Phillios (2003a"b') and with 
the CPT in lSenechal and Trcmblav 12003) . 

Figure EH shows the DCA/NCA result for the low- 
energy DOS in the 2D Hubbard model at 5% doping 
{5 — 0.05) around the chemical potential {lo = 0) cal- 



For a review on the pseudogap phenomenon see 
iTimusk and Statd 119991) . 
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culated on a 4-site cluster {N^ = 4). At high tempera- 
ture T > 0.3t no pseudogap is seen in the DOS. As the 
temperature is lowered the DOS distorts at the chemi- 
cal potential and a pseudogap emerges at lu = when 
T < 0.3t. 



T/t=0.400 

T/t=0.270 

T/t=0.068 




Figure 31 DOS for various temperatures in the 2D Hubbard 
model at 5% doping when U — 12t near the chemical potential 
[uj = 0) calcul ated wi t h DC A/NCA for a 4-site cluster, Nc = 
4. Taken from lMaieJ (1200 ll) . 



These perturbative res ults^ were confirm ed in 
DCA/QMC simulations by iJarrell et all l|2001a|) which 
we reproduce in Fig. At low temperatures a pseudo- 
gap is observed in the DOS at dopings S < 0.2. This de- 
pression of quasiparticle states at the chemical potential 
is accompanied by a downturn of the uniform magnetic 
susceptibility shown in the inset. For low to intermediate 
doping it develops a maximum defining a crossover tem- 
perature T* . Below T* quasiparticle and low-energy spin 
excitations are suppressed by non-local correlations sim- 
ilar to what is observed in the experiment. At the same 
time, the charge susceptibility (not shown) displays qual- 
itatively different behavior forming a strong low-energy 
peak at low temperatures fsee .Maier et al. . 2004. 2 002al 
and Fig. I38|l which signals the emergence of coherent 
charge excitations below T*. The TVc = 4 DC A results 
are thus consistent with a spin-charge separated picture 
as in Anderson's RVB theory. This is not surprising since, 
as we discussed in Sec. lIV.D.2l fluctuations are enhanced 
in the Nc — 4: cluster due to the "too small" coordination 
number. It is known that small coordinatio n numbers 
favor the spin-charge separated RVB state l)AndersonL 




Figure 32 DOS for various dopings 5 in the 2D Hubbard 
model at r = 0.092t and U ^ 8t calculated with DCA/QMC 
for a 4-site cluster, Nc — 4. Inset: Uniform spin- susceptibility 
as a f unction of temperature. Taken from IJarrell et al\ 
ll2001af) . 



I1987D over the Neel state. Indeed the RVB state was 
shown to be the ground state of a 2 x 2 Hcisenberg model 
with periodic bounda ry conditions with a large gap to 
the first excited state (jDaeotto and Moreol Il988|) . 
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The NCA is perturbative in the coupling to the host, not in the 
interaction, as detailed in Sec. 1111. B. 31 



Figure 33 The spectral function A(k, ui) near the chemical 
potential in the 2D Hubbard model at T = 0.088t and U = 
8t along high-symmetry directions in the first Brillouin-zone 
between F = (0,0), X — (7r,0) and M — (7r,7r) calculated 
with DCA/QMC for a 4-site cluster Nc = 4. 

The corresponding 4-site DCA/QMC result for the mo- 
mentum resolved spectral function A(k, lu) in the pseu- 
dogap regime is displayed in Fig. 1331 for energies near 
the chemical potential between points of high symmetry 
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r = (0,0), X = (7r,0) and M = (tt, tt) in the first Bril- 
louin zone. The overall dispersion of the band crossing 
the chemical potential uj — follows that of the non- 
interacting system et- While coherent quasiparticles ex- 
ist along r M, the pseudogap is seen near X — (tt, 0) at 
the chemical potential uj = 0. The anisotropy of the pseu- 
dogap is thus consistent with that observed in ARPES 
measurements on underdoped hole-doped cuprates. 

Qualitatively similar results for the emergence of the 
quasiparticle spectrum in the doped 2D Hubbard model 
were obtained by Stanescu and Phillios (2003b) using 
the two-site com posite operator approach disc ussed in 
Sec. HLHI (see also IStanescu and Phillip l2003a|) . Figure 
1341 illustrates their results for the doping dependence of 
the chemical potential /i, the imaginary part of the self- 
energy SmS and the U dependence of the low-energy 
DOS. At half-filling the chemical potential has a dis- 
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Figure 34 The doping dependence of the chemical potential in 
the 2D Hubbard model calculated with the two-site composite 
operator method for T = 0.15t (dashed line) and T = 0.07t 
(solid line). Right inset: DOS for various values of U at 5% 
doping (n = 0.95). Left inset: Imaginary part of the self- 
energy evaluated at a Fermi momentum (0.3,2.10) for n = 
0.97. (0.3,1.84) for n = 0.80 and (0 3,1.06) for n = 0.3. Taken 
from IStanescu and Phillip^ ll2003bl) . 



continuity indicating the abs ence of mid-gap state s. In 
agreement w ith DCA/QM C (jJarrell et ad. l2001albl ) and 



DCA/NCA l)Maieil l200lD results, |StoE| is large in the 
underdoped pseudogap regime (n = 0.97) and acquires 
Fermi-liquid behavior at larger doping n < 0.80 indicated 
by the parabolic minimum at the chemical potential. As 
illustrated in the inset, the depth of the pseudogap de- 
creases as U increases suggesting a pseudogap scale com- 
p atible with t'^/U. 

ISenechal and Tremblavl l)2003l) recently investigated 
the difference in pseudogap behavior between electron 
and hole doped HTSC using CPT for the 2V) t ~ t' - t" 
Hubbard model. As illustrated in the Fermi-surface plots 
in Fig. 123 their results a± U — 'Ht demonstrate that the 
pseudogap in hole-doped systems (right side) occurs near 
X = (tt, 0) at optimal doping consistent with the results 



discussed above. In electron-doped systems (left side) 
however, the pseudogap appears at the crossing points 
of the Fermi surface with the antiferromagnetic Brillouin 
zone boundary at moderate interaction U — At. When U 
is large however (not shown) , the Fermi surface only sur- 
vives in the neighborhood of (tt, 0) and (0, tt). As shown 
in the lower panels of Fig.|23the pseudogap in both cases 
is generated by a large scattering rate |3r7iS](k, 0)| at 
the chemical potential. A unified picture of the spectral 



o 
< 



U=4 
















(0,0) 












Figure 35 Intensity plot of the spectral function A(k, 0) (top) 
and imaginary part of the self-energy £smE(k, 0) (bottom) of 
the 2D Hubbard model calculated with CPT on a 3 x 4-site 
cluster. The left side shows the results in the first quadrant 
of the Brillouin zone for the 17% electron-doped system at 
U — it and the rig ht side for the 17% hole-doped system at 
U = 8t. Taken from ISenechal and Tremblavl ||2003D . 

properties of the electron- and hole-doped cuprates thus 
emerges from these results if the interaction strength U 
is allowed to be doping dependent. To reproduce the 
experimental observations in optimally doped cuprates, 
large values of U seem necessary in hole-doped {U « 8t) 
systems, while smaller va lues of U describe the electron - 
doped systems {U < 6t) l)Senechal and Tremblavl l2003|) . 



4. Superconductivity 

FSS QMC simulations for the doped 2D Hubbard 
model in the intermediate coupling regime U ^ W sup- 
port the idea of a spin-fiuctuation driven interaction 
mediating d-wa ve superconductivity (for a review see 
iScalapinol Il999l) . However the fermion sign problem lim- 
its these calculations to temperatures too high to study 
a possible transition. These calculations are also re- 
stricted to relatively small system sizes, making state- 
ments for the thermodynamic limit problematic, and in- 
hibiting studies of the low-energy physics. These short- 
comings do not apply to embedded cluster theories which 
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are built for the thermodynamic limit. Cluster sizes 
larger than one are necessary however, to describe a pos- 
sible transition to a state with a non-local (d-wave) order 
parameter as discussed in Sec. III. El 

It is well known from weak coupling FSS FLEX re- 
and phenomen ological theories 



suits llBickers et pi.Lll98fl| )^ 
fMo nthoux et a/.lll99ll: IS 



TScalapinol Il999l) that antiferro- 
magnetic spin-fluctuations mediate pairing with d-wave 
symmetry and cause a pseudogap in underdoped systems. 
In optimally doped cuprates, these spin- fluctuations are 
known to be short-ranged, extending over a few lat- 
tice spacings. Hence quantum cluster approaches should 
provide an adequate methodology to study supercon- 
ductivity in these systems. Pairing in t he 2D Hub- 
bard model was studied using DC A/NCAbviM aier et a,l\ 
l2000a). and with DCA /QMC bv l.Tarrell etai (,.2001a«l¥ 
and lMaier et all ^00^. The possible coexistence of su- 
perconductivity with antifcrromagnetic order was inves- 
tigated by Lichtenstein and Kat snclson ( 2000). 

The results of 4-site (A^c — 4) DCA simulations for 
the doped 2D Hubbard model show an instability to a 
superconducting phase with a d^2_y2--w&Ye order param- 
eter at low enough temperatures. As a typical example 
of this transition. Fig. 1361 shows the DCA/NCA result 
for the DOS and the coarse-grained anomalous Green 
function Gi2{'K,uj) = F(K,w) defined in Eq. at 
different cluster K-points near the chemical potential in 
the superconducting state. Gi2(K,a;) vanishes at the 
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Figure 36 DOS near the chemical potential (a) and coarse- 
grained anomalous Green function Gi2(K.,lj) = _F(K,(.<j) for 
different DCA cells (b-c) in the superconducting state of the 
2D Hubbard model at 19% doping, T = 0M7t, U = 12t 
calc ulated with DCA/NC A for a 4-site cluster, Ac = 4. Taken 
from lMaier et al\ ll2000aD . 

cluster K = (0,0) and (tTjTt) but is finite at (tt, 0) and 
(0, tt) with opposite signs. Since the K-dependence of the 
coarse-grained order parameter Ak is given by the K- 
dependence of the coarse-grained anomalous Green func- 
tion (see Eq. H99|l '). this result is consistent with a d^2_y2- 
symmetry of the order parameter. The finite pair ampli- 



tude is also reflected in the DOS depicted in the upper 
left part where the lower subband of the full spectrum 
is shown. It displays the superconducting pseudogap at 
z ero frequency as expected for a d-wave order parameter. 

iJarrell et aU (t2001a.b.) used DCA/QMC to search for 
many different types of superconductivity, including s- 
, extended s-, p- and d-wave, of both even and odd 
frequency. Of these, only the odd-frequency s-wave 
and even-frequency d-wave pair-field susceptibilities were 
strongly enhanced, and only the d-wave susceptibility di- 
verged. This is illustrated in Fig. 1371 where the pair-field 
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Figure 37 Pair-field susceptibilities versus temperature in the 
even frequency s-wave, extended s-wave (xs), d-dave and odd- 
frequency s-wave channels in the 2D Hubbard model at 5% 
doping, U = 8t calculated with DCA/QMC for a 4-site clus- 
ter, Nc = 4. Inset: Inverse d-wave pair-field susceptibility ver- 
sus temperature for different dopings and cluster sizes. The 
line is a fit to b(T - T^V with Tc = 0.084t and 7 = 0.72. 
Taken from iJarrefi et aZ.H2001al'l . 

susceptibilities are plotted versus temperature at 5% dop- 
ing. The behavior of the inverse d-wave pair-field suscep- 
tibility as a function of temperature for Nc — \ and 4 and 
5 = Q and 0.5 is shown in the inset. For Nc = 1 there 
is no tendency towards pairing. As detailed in Sec. III. El 
the DMFT is not able to describe pairing with symme- 
tries lower than the lattice symmetry (i.e., p-, d-wave, 
etc.). For Nc ~ ^ and (5 = the inverse susceptibil- 
ity rises abruptly as the temperature is lowered and the 
Mott gap opens in the DOS. The Mott gap becomes more 
pronounced as Nc increases (see Sec. lIV.LrTll . so that for 
larger clusters the gap prevents superconductivity even 
for U < W. If charge excitations are gapped, then pair- 
ing is suppressed. At half-filling, for U ^ 8t the gap 
is of order U, and thus much larger than the magnetic 
exchange energy J ~ At'^/U — 0.5t. Hence the opening 
of the Mott gap suppresses any magnetically mediated 
pairing. Away from half-filling the width of the pseudo- 
gap in the charge excitation spectrum is much smaller, 
of the order of J (see Sec. IIV.D.3|) . so that magnetically 
mediated pairing is possible. 
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More insight in the nature of pairing was gained from 
further D CA/QMC studies of the 2D Hubbard model 
ijMaier et a l. 2004). Figure EHl shows the DCA/QMC re- 



energies, Tr(tG) and Tr(SG) respectively, of the super- 
conducting (SC) and normal state (NS) solution as a 
function of temperature at low doping 5 — 0.05 (left) and 
optimal doping 5 = 0.20 (right). For both doping levels, 
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Figure 38 DOS (left), local dynamic spin susceptibility (cen- 
ter) and local dynamic charge susceptibility (right) in the 2D 
Hubbard model at 5% doping at different temperatures above 
and below the critical temperature Tc = 0.0218 = 0.087t cal- 
culat ed with DCA/QM C for a 4-site cluster, Nc = 4. Taken 
from lMaier et al\ ll2004f) . 
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suit for the evolution of the DOS, the local dynamic spin- 
and local dynamic charge-susceptibility^ as the tempera- 
ture decreases below the critical temperature Tc. As dis- 
cussed in Sec. IIV.D.3I the normal state low-temperature 
DOS and spin-susceptibility display a pseudogap, i.e. a 
depression of low-energy quasiparticle and spin excita- 
tions. Both quantities evolve smoothly across the super- 
conducting transition with the pseudogap changing to 
a superconducting gap^ below Tc. However since the 
charge susceptibility is peaked at zero frequency even 
slightly above Tc, it changes abruptly upon pairing to 
show the same behavior as the spin-susceptibility, includ- 
ing the superconducting gap at low frequencies. Remark- 
ably, well below Tc all quantities display narrow peaks at 
LD « O.leV, delimiting the superconducting gap. This re- 
sult clearly indicates the formation of quasiparticles be- 
low Tc. The absence of quasiparticles in the normal state 
however undermines the very foundation of the BCS the- 
ory of conventional superconductors where pairing is a 
result of a Fermi surface instability that rel ies on the 
existe nce of quasiparticles in a Fermi-liquid ijSchriefferl 
I1993D . 

DCA/QMC results for the condensation energy further 
establish the unconventional c haracter of superco nduc- 
tivity in the 2D Hubbard model ijMaier et a/.ll2004|) . Fig- 
ure 123 presents the kinetic (top) and potential (bottom) 



Note that in the DCA, local quantities are identical in the lattice 
and on the cluster and thus are easily calculated. 
Note that due to the finite resolution in momentum space, the 
DCA underestimates low-energy spectral weight in superconduc- 
tors where the gap has nodes on the Fermi surface. As a result, 
a fully developed gap is found at low temperatures instead of a 
DOS that vanishes linearly in frequency as expected for a d-wave 
superconductor. 



Figure 39 Kinetic (top) and potential (bottom) energies ver- 
sus temperature in the normal (NS) and superconducting 
(SC) states in the 2D Hubbard model at 5% (left) and 20% 
(right) doping for U — 8t. Tc is indicated by the vertical 
dotted hues. Taken from lMaier ei al\ (12004^ . 

the kinetic energy in the superconducting state is reduced 
compared to the normal state, while the potential ener- 
gies are almost identical. This result is in agreement with 
recent optical experiments which show that the supercon- 
ducting transition in the cupr ates is due to a lowerin g of 
the electronic kinetic energy l)Molegraaf et oi.L|2fl02ft . It 
further supports the evidence that pairing in the Hub- 
bard model is fundamentally different from BCS pairing 
which occurs through a reduction of the electronic po- 
tential energy accompanied by a slight increase in the 
kinetic energy. 

The possibility of coexisting d-wave superconducting 
and antiferromagne tic order in the 2D Hubbard model 
was investigated by iLichtenstein and KatsnelsonI (^2000) 
using a 4-site cluster approach similar to the DCA/QMC 
method (see Fig. EHIa)). In this approach, an 8 x 8 
matrix representation of the Green function is required 
to account for both the antiferromagnetic order param- 
eter (c||Cj|) and the superconducting order parameter 
(cijCji). Figure HUT c) reproduces the results for the two 
order parameters as a function of doping at fixed tem- 
perature in the weak coupling regime {U = 4.8t). The 
authors find that the antiferromagnetic order parameter 
coexists with the d-wave superconducting order param- 
eter over a wide range of doping. Consistent with the 
DCA/QMC results, the antiferromagnetic order parame- 
ter is maximal at zero doping where the superconducting 
order parameter vanishes due to the opening of the gap. 
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Figure 40 (a) Schematic representation of an antiferromag- 
netic d-wave 2x2 periodically repeated cluster, (b) Generic 
phase diagram of HTSC. (c) Magnetic (M) and d-wave super- 
conducting (F) order parameters versus hole-doping in the 
2D Hubbard model at l3t = 15, t' = -0.15t, U = 4.8t calcu- 
lated with a 4-site cluster approach simila r to th e DCA/QMC. 
Taken from iLichtenstein and KatsnelsonI l(2C 



5. Phase diagram 




Figure 41 Temperature-doping phase diagram of the 2D Hub- 
bard model when U — 8t calculated with DCA/QMC for a 
4-site cluster, Nc — 4. The error bars on T* result from 
the difficulty in locating the maximum in the uniform spin 
susceptibility. Regions of antiferromagnetism, d-wave super- 
conducting and pseudogap behavior are seen. 



The results reviewed in the preceding sections illus- 
trate that quantum cluster approaches applied to the 2D 
Hubbard model are able to capture the complex behav- 
ior observed in HTSC. The qualitative agreement with 
experiments is summarized in the Nc = 4 DCA/QMC 
temperature-doping (T-S) phase diagram of the 2D Hub- 
bard model in the intermediate coupling regime U — W 
shown in Fig. 1411 The phase boundaries were determined 
by the instabilities of the paramagnetic phase as indi- 
cated by the divergence of the corresponding susceptibil- 
ities. Therefore these results do not allow any conclu- 
sions about a possible coexistence of the antiferromag- 
netic and d-w ave superconducting phases for S < 0.5. 
The results bv ILichtenstein and KatsnelsonI l)200(]() how- 
ever suggest this coexistence at least for weak coupling 
(see Fig. note)). 

The pseudogap crossover temperature T* determined 
by the peak in the uniform spin susceptibility (see 
Fig. I32|l serves as a boundary separating the observed 
Fermi liquid and non- Fermi liquid behavior. For x < T* 
the self-energy shows non-Fermi liquid character for the 
parts on the Fermi surface near k = (tt, 0) (see Fig. I34II . 
Quasiparticle and low-energy spin excitations are sup- 
pressed as indicated by the pseudogap in the DOS and 
the spin-susceptibility (see Fig. I38|l . At S > 0.2 Fermi- 
liquid behavior is recovered (see Fig. l34|l . At low temper- 
atures, the systems is antiferromagnetic near half-filling. 
The d-wave superconducting phase at finite doping has 
its maximum transition temperature at (5 ~ 0.05. 

As indicated in Sec. IIV.D.3I the = 4 DC A cluster 
favors the spin-charge separated RVB state. The A^'c = 4 
results may thus be interpreted within the RVB picture 
iAndcrso il 1198/1) : The pairing of spins in singlets be- 
low the crossover temperature T* results in the suppres- 
sion of low-energy spin excitations and consequently in 



a pseudogap in the density of states. Charge excitations 
are quasi-free as indicated by the zero-frequency peak in 
the charge susceptibility (see Fig. IS5|). Well below the 
transition spin and charge degrees of freedom recombine, 
forming electrons which pair. Frustrated kinetic energy 
is recovered as indicated by the reduction of the kinetic 
energy as the system goes superconducting (see Fig. I39|l . 

Although the properties of HTSC are well described 
by the Nc — 4 results, it is important to ask the ques- 
tion whether the phase diagram and in particular the 
observed RVB nature of the results are stable when the 
cluster size is increased. While the pseudogap temper- 
ature T* may be expected to persist with increasing 
Nc, the transition temperatures to the antiferromagnetic 
and superconducting state are expected to fall to zero, 
as discussed in Sec. II. Bl (see also Fig. I?7jl . This has 
been confirnied thr ough simulations of larger clusters 
ijjarrell et all . l2004|) . A finite inter-planar coupling be- 
tween an infinite set of Hubbard planes may then be used 
to stabilize the transit ions at larger Nc, at least to the an- 
tiferromagnetic state l|.Tarrell et oiJ . l2001lJ) . The stability 
of the RVB behavior and even the superconductivity is 
less clear and more complex models may be necessary for 
them to persist in larger clusters. 

6. Studies of related models 

a. Stripes in the t-J-model. Several numerical studies 
indicate that there is a tendency for doped holes to 
form stripes separated by antiferromagnetic d omains in 
stron gly correlated systems (for a review see iDagottol 
Il994l) . CPT is the quantum cluster method of choice 
to study the large unit cells in the inhomogeneous stripe 
phase, which complicate if not preclude the application of 
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the numerically more expensive embedded cluster tech- 
niques (CDMFT, DCA). A thorough study of stripes 
in the large U limit of the Hubbard model, the t-J 
model, was conducted bv^ Zacher et al\ |2002a"b') within 
the CPT. To implement the CPT for this problem, the 
authors divided the lattice into alternating clusters of 
metallic stripes and antiferromagnetic domains. The 
inter-cluster hopping linking these clusters was treated 
perturbatively within the CPT. The enforced stripe pat- 
tern in this implementation prohibits to explore the sta- 
bility of stripes, but allows to investigate the effects of the 
stripe pattern on the single-particle excitations. In sys- 
tems with less than 12% doping the technique was shown 
to reproduce salient ARPES features in selected HTSC 
if a site-centered 3-f 1 stripe pattern, i.e. half-filled anti- 
ferromagnetic three-leg ladders separated by doped one- 
leg chains, was chosen (jZacher et a/.l l2002b() . At higher 
dopings the comparison with ARPES indicates that the 
weight of bond-centered stripes with a 2-f 2 pattern in- 
creases in which excess holes prolif erate out of the stripe s 
into the antiferromagnetic domain l)Zacher et a/J . l2002a(l . 



b. Spectral properties of the three-band Hubbard model. In 
the context of HTSC, the single-band Hubbard model can 
be viewed as a low-energy approximation of the more 
complex and more realistic three-band Hubbard model 
l|Emervl[l987tlEmerv and Reiteill 1988(1 . The three-band 
Hubbard model takes into account the Px and py oxygen 
orbitals in addition to the Cu d degrees of freedom in 
the superconducting Cu0 2 planes. The CPT s tudy of 
its spectral properties by 'Da hnken et al\ l|2002f) shows 
very good agreement with ARPES data on HTSC at half- 
filling as well as in the doped system including a hole- 
like Fermi surface at high doping which splits into an 
electron- and hole-like branch when a bilayer hopping t± 
is included. 



c. Cluster simulations of the periodic Anderson model. The 
periodic Anderson model (PAM) is widely considered to 
be a paradigm for a variety of rare-earth and actinide 
compounds including the heavy Fermion systems. It is 
composed of a strongly hybridizing band of d electrons 
and a weakly hybridizing band of correlated /-electrons 
described by the Hamiltonian 

H = -tJ2idld,, + h.c.) + Y,^d{nf^+0 

ij,cr i 
i i 

DMFT simulations of the PAM reveal an antiferromag- 
netic insulating phase at half filling of both the / and 
d bands. The gap is set by the Kondo coherence scale 
To which is strongly enhanced compared to the single- 
impurity model (SIAM) scale. When the c?-band is doped 
away from half filling while the /-band remains roughly 



half filled, the system becomes more metallic, and the 
Kondo scale is strongly suppressed compared to that of 
the SIAM. In both the insulator and metal, the tempera- 
ture dependence of the impurity susceptibility and many 
other low temper ature quantities deviate strongly from 
that of the SIAM l|Tahvildar-Zadeh et ad ll997|). 

Non-local corrections were studied bv IShimizul l|2002|) 
who used the DCA together with the NCA as a clus- 
ter solver to study the single-particle spectra. He finds 
large deviations from the DMFT result due to the ef- 
fects of RKKY exchange. At half filling, where the 
RKKY exchange is strong and antiferromagnetic, he finds 
a large gap of the order of the RKKY exchange energy. 
When the filling of the conduction band is small and the 
RKKY exchange is weaker and ferromagnetic, the coher- 
ence peak is restored. 



V. CONCLUSIONS AND PERSPECTIVES 

In this review we tried to convey the message that 
quantum cluster approaches provide powerful theoreti- 
cal tools to study the rich phenomenology in systems 
dominated by strong electronic interactions, such as most 
notably transition metal oxides, heavy Fermion and one- 
dimensional systems including superconducting and mag- 
netic compounds. Quantum cluster approaches are non- 
perturbative in nature, their quality can be systemati- 
cally improved by increasing the cluster size, and they 
provide complementary information to finite size simula- 
tions. By mapping the lattice problem to a finite size 
cluster they describe short-ranged correlations within 
the cluster accurately while approximating longer-ranged 
physics on the mean-field level. Of the various attempts 
to add non-local corrections to local approximations such 
as DMFT, we focused in this review on three established 
quantum cluster approaches which we believe to play the 
major role in the description of many-particle systems. 

The cluster perturbation theory provides a very eco- 
nomic way to calculate the single-particle dynamics by 
continuing the results of an isolated finite size cluster to 
the thermodynamic limit. For a cluster consisting of a 
single site only this method is identical to the Hubbard-I 
approximation, while it recovers the exact result in the 
infinite size cluster limit. When combined with the self- 
energy functional approach, the CPT can also be used to 
study instabilities to symmetry broken phases. 

Both embedded cluster techniques, the dynamical clus- 
ter approximation and the cellular dynamical mean field 
theory are superior to the CPT in that they map the 
lattice to an embedded cluster instead of the CPT's iso- 
lated cluster. This leads to a self-consistent theory with 
a single-particle coupling between the cluster and the 
host. As a result both DCA and CDMFT naturally allow 
for the study of phase-transitions and they provide ther- 
modynamically consistent results on the one- and two- 
particle level. 

CDMFT can be viewed as a direct generalization of 
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DMFT to a cluster in real space. The mapping between 
lattice and cluster problems is identical to that of the 
long established MCPA for disordered systems. It leads 
to a cluster with open boundary conditions which vi- 
olates translational symmetries. In contrast the DCA 
cluster is defined in cluster reciprocal space. Hence the 
DCA cluster has periodic boundary conditions and there- 
fore preserves the translational symmetries of the lattice. 
This difference in boundary conditions translates directly 
to different asymptotic behaviors for large linear cluster 
sizes Lcj smd the decision which method to use for a given 
problem strongly depends on the quantities of interest: 
Local quantities, such as the local density of states when 
defined on central cluster sites converge faster in the 
CDMFT since they do not directly couple to the mean- 
field. Due to the large mean-field coupling of the surface 
sites however, the CDMFT converges slowly, with correc- 
tions of order 0{1/Lc), for quantities extended over the 
cluster. The DCA converges more quickly, as 0(1/L^), 
due to the periodic boundary conditions on the cluster. 

Quantum cluster approaches reduce the complexity of 
the infinite lattice problem by mapping it to a cluster 
with fewer degrees of freedom. The numerous methods 
employed to solve the DMFT equations are in princi- 
ple available to study the effective cluster model. How- 
ever as the complexity of this task rapidly increases 
with cluster size, potential cluster solvers are faced with 
severe size limitations. As the most promising tech- 
niques we reviewed perturbative approaches including 
the fluctuation-exchange approximation and the non- 
crossing approximation as well as non-perturbative tech- 
niques including quantum Monte Carlo and numerical 
renormalization group. 

We discussed the application of quantum cluster ap- 
proaches to a wide range of problems in condensed matter 
physics. The information gained from these studies has 
led to significant progress in the field of strongly corre- 
lated electron systems. Even studies using small cluster 
sizes opened up new insight in problems such as one- 
dimensional systems, the Mott-Hubbard transition in two 
dimensions and high-temperature superconductivity. 

To gain conclusive evidence however, larger cluster size 
studies are inevitable to verify the information obtained 
from small clusters. This task is severely hampered by 
the rapidly increasing complexity with cluster size. Fu- 
ture progress therefore is closely linked to improvements 
in the efficiency and flexibility of the techniques used 
to solve the effective cluster problem. Within quantum 
cluster approaches we explore a coarse-graining approx- 
imation in fc-space. To further reduce the complexity, 
the same idea could be transfered to the frequency do- 
main. By coarse-graining the frequency dependence of 
irreducible quantities, correlations on long time scales are 
neglected, while the short time scale behavior is described 
accurately. An important aspect in this context is again 
the causality question. First test results are encourag- 
ing; they show that coarse-graining in Matsubara space 
leads to acausalities, while coarse-graining the real fre- 



quency axis does not face this problem l)Arvanpour et all . 
I20n8al) . 

Another route to defeat the cluster size problem is 
to develop hybrid algorithms that treat different length- 
scales in the problem with di fferent accur acy. As a 
promising step in this direction. iHaeue et al\ ((2003) have 
developed a hybrid technique which maps the infinite lat- 
tice onto two embedded clusters of different size, thus 
dividing the problem into three length-scales. Short- 
ranged correlations described by the small cluster are 
treated accurately within QMC, correlations of interme- 
diate length-scale are treated perturbatively in the large 
cluster using FLEX, and the long-ranged physics beyond 
the size of the larger cluster is again approximated on 
the mean-field level. 

To improve comparisons with experiments and to 
achieve predictive capability, the inclusion of the specifics 
of the actual materials is required. Along the lines of the 
LDA-f DMFT approach, one can use electronic structure 
calculations to parameterize the models studied by quan- 
tum cluster approaches. First steps in this direction have 
been made by [Potervaev et al. (2003). A more integrated 
approach to the ab-initio description of strongly corre- 
lated systems by combining the ideas of density func- 
tional theory and quantum cluster approaches remains 
an important and challenging task. 
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